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ABSTRACT 


This  is  the  first  volume  of  the  final  report  in  a  study  of  reentry  effects  on  signal  transmission 
from  slender  cone  bodies  at  ICBM  velocities.  The  report  gives  results  of  state  of  knowledge 
surveys  of  fluid  mechanics  methods  and  electromagnetic  phenomena.  Results  of  illustrative 
calculations  are  also  given  in  many  of  these  areas.  Antenna  breakdown  is  found  to  be  the 
most  important  problem  in  reentry  transmission  effects,  as  well  as  the  least  well  understood. 
Recommendations  for  needed  research  are  given,  including  hig  .  gas  temperature  and 
convection  experiments  and  theoretical  analyses  of  temperature,  convection,  and  flow  field 
gradients  effects.  High  altitude  flow  field  development  and  studies  of  thermal  nonequilibrium 
effects  are  recommended  also. 
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SECTION  1 


INTRODUCTION 


This  is  the  final  report  of  a  study  of  reentry  effects  on  ECM  antenna  systems,  in  which  data 
have  been  collected  with  the  aim  of  aiding  the  antenna  designer  predict  and,  if  possible, 
avoid  such  reentry  effects.  Much  of  these  data  are  qualitative  and  some  of  it  based  on  poorly 
understood  phenomena  in  the  present  state  of  knowledge  of  the  limiting  technologies.  Thus 
there  is  the  need  for  critical  evaluation  and  the  formulation  of  recommendations  for  further 
study,  as  have  been  done. 

1.1  OBJECTIVE  OF  STUDY 

The  objective  of  this  study  is  the  preparation  of  a  document  to  serve  as  a  compendium  of 
design  data  and  principles  for  .’ntennas  compatible  with  slender  body  vehicles  which  reenter 
at  ICBM  velocity.  This  study  gives  consideration  to  flow  field  analysis,  antennas,  linear 
plasma  effects,  antenna  breakdown,  nonlinear  attenuation,  alleviation  techniques,  and  the 
preliminary  design  of  recommended  flight  experiments.  In  most  of  these  areas  the  emphasis 
is  placed  on  a  survey  of  present  knowledge  and  methods  of  prediction  and  control.  Limited 
illustrative  calculations  have  been  carried  out  in  many  cases  The  problem  of  reentry 
effects  on  antenna  breakdown  is  stressed. 

1.2  B/CKGROUND 

The  above  objective  is  stated  in  terms  of  a  technology  problem,  the  solution  of  which  has 
application  to  the  design  of  antennas  for  reentry  vehicles.  The  basic  problem  is  the  predic¬ 
tion  and. control  of  the  effects  of  the  reentry  induced  environment  oh  the  antenna  performance. 
Attention  is  limited  to  the  reentry  flow  field  of  a  slender  cone  vehicle  reentering  the  atmos¬ 
phere  with  ICBM  velocity.  The  aspects  of  the  flow  field  which  are  important  to  the  antenna 
performance  are  the  profiles  of  gas  species  concentrations,  including  excited  states,  tempera¬ 
tures,  velocities,  and  electron  density  in  the  body  flow  field  (wake  flow,  except  possibly  base 
flow  and  near  wake  in  rare  cases,  is  unimportant).  An  important  part  of  the  problem  of  flow 
field  characterization  is  the  variation  with  altitude  during  reentry  of  each  of  the  above 
quantities. 


The  presence  of  appreciable  levels  of  ionization  in  the  reentry  flow  field  can  lead  to  various 
linear  plasma  effects.  Linear  effects,  which  are  defined  as  being  independent  of  the  antenna 
input  power,  can  be  expected  at  low  input  signal  levels.  They  include  but  are  not  limited  to 
attenuation,  antenna  input  impedance  changes,  and  antenna  pattern  distortion.  These  effects 
depend  on  the  electron  density,  gas  density,  and  temperature  profiles  over  the  body. 


When  the  input  signal  level  to  the  antenna  cannot  be  assumed  to  be  arbitrarily  small,  then  it 
is  possible  for  nonlinear  plasma  effects  and  breakdown  to  be  important.  Nonlinear  plasma 
effects  can  take  many  more  forms  than  linear  effects,  in  that  the  signals  which  are  input  to 
the  antenna  can  change  the  electromagnetic  characteristics  of  the  flow  field  in  many  different 
ways  which  depend  on  the  nature  of  the  signals  as  well  as  on  the  flow  field  properties.  How¬ 
ever  the  same  basic  flow  field  properties  are  important  in  nonlinear  as  in  linear  effects, 
provided  breakdown  has  not  yet  been  reached. 


Antenna  breakdown  is  reached  when  the  input  signal  level  becomes  so  strong  that  the  antenna 
field  ionizes  the  air  around  the  antenna  to  the  extent  that  signals  can  not  propagate  through 
to  the  receiver.  This  phenomenon  takes  place  in  a  static  cold  air  environment  at  a  given 
altitude  even  without  the  effects  of  the  reentry  environment;  however,  the  presence  of 
hypersonic  flow  field  conditions  can  change  significantly  the  input  signal  level  at  which  it 
takes  place.  Since  cold  air  static  breakdown  can  easily  be  measured  in  a  laboratory  vacuum 
chamber  for  a  given  antenna,  the  important  thing  to  be  determined  about  the  effects  of  the 
reentry  environment  is  just  what  change  in  the  breakdown  signal  level  versus  altitude  they 
will  cause.  Here  the  antenna  field  is  a  cause  of  electron  sources  and  sinks  as  a  function 
of  the  gas  species  concentrations,  excitation  levels,  and  velocity;  thus  these  parameters 
as  well  as  flow  field  ionization  must  be  established  to  determine  the  reentry  effects.  In 
addition,  the  altitude  region  in  which  breakdown  is  important  is  so  high  (especially  at  the 
lower  end  of  the  frequency  spectrum  of  reentry  signal  transmission  systems)  as  to  require 
study  of  the  limits  of  the  continuum  regime  of  fluid  mechanics. 
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1.3  FORMAT  OF  REPORT 

After  a  short  summary  of  the  findings  of  this  study,  this  report  gives  each  technical  prob¬ 
lem  area  in  the  form  of  a  survey,  followed  by  a  few  illustrative  calculations.  The  flow 
field  section  treats  a  limited  number  of  calculations  of  importance  in  antenna  breakdown 
for  a  typical  vehicle  configuration  and  velocity-altitude  trajectory.  The  main  emphasis  on 
antennas  is  the  estimation  of  diffusion  length  for  use  in  breakdown  calculations.  In  the 
section  on  linear  plasma  effects,  the  electron  density  levels  which  would  be  detrimental 
to  the  proper  functioning  of  antennas  on  conical  vehicles  are  indicated  so  that  the  vehicle 
designer  can  determine  to  what  extent  he  should  attempt  to  reduce  the  flow  field  plasma 
level.  A  complete  survey  of  the  theory  of  classical  breakdown  and  of  reentry  effects  is 
given,  followed  by  illustrative  calculations  of  normalized  antenna  breakdown  power  in  a 
typical  reentry  environment.  This  section  also  presents  a  convenient  closed  form  expression 
for  breakdown  power  of  an  antenna.  The  section  on  nonlinear  effects  is  concerned  mainly 
with  a  survey  intended  to  assess  the  relative  importance  of  these  effects  without  going  into 
details  of  the  theory,  since  this  is  a  tremendously  large  and  complicated  subject.  Alleviation 
techniques  are  surveyed  with  the  purpose  of  indicating  which  might  be  useful.  Recommenda¬ 
tions  are  given  on  the  importance  of  various  aspects  of  the  problem  and  on  further  research 
which  is  needed  to  provide  definitive  design  data  and  systems  limitations.  Classified  aspects 
are  given  in  the  second  volume,  including  analyses  of  previous  flight  experiments  and 
recommendations  for  future  ones. 


SECTION  2 
SUMMARY 


The  electron  density  levels  at  which  linear  and  nonlinear  plasma  effects  on  signal  transmission 
from  a  slender  cone  reentry  body  become  important  are  evaluated.  Flow  field  technology 
and  illustrative  calculations  show  that  these  plasma  effects  can  be  avoided  on  a  properly 
designed  vehicle.  A  survey  of  reentry  antenna  breakdown  brings  out  many  important  factors. 
The  state  of  knowledge  of  antenna  breakdown  effects  is  such  as  to  leave  several  important 
technology  problems  still  to  be  solved.  There  is  a  need  for  more  experimental  and  theoretical 
research  in  these  areas.  However,  a  closed  form  expression  for  reentry  antenna  break¬ 
down  power  is  given.  This  expression  is  sufficiently  general  to  fit  the  available  applicable 
experimental  data,  and  it  is  probably  accurate  enough  to  give  crude  but  conservative  extimaces 
of  reentry  effects. 
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SECTION  3 

FLOW  FIELD  ANALYSIS 


3.1  SURVEY  OF  FLOW  FIELD  TECHNOLOGY 

A  literature  survey  of  theorectical  techniques  which  are  capable  of  describing  the  flow  field 
environment  about  a  body  during  atmospheric  reentry  is  presented  in  the  form  of  a  categorized 
review  and  reference  list.  The  survey  was  limited  to  within  the  United  States,  and  the  sources 
of  information  include  the  prominent  journals  and  texts.  In  addition,  thirty-seven  organiza¬ 
tions  and  individuals  known  to  be  involved  in  this  type  of  work  were  solicited  for  inputs.  How¬ 
ever,  because  only  seventeen  of  these  requests  were  acknowledged,  the  survey,  although 
extensive,  cannot  be  considered  complete.  It  is  restricted  in  scope  to  axisymmetric  or 
three-dimensional  bodies  within  the  continuum  regime  of  reentry  and  to  flows  which  are 
steady,  compressible, and  attached  (as  distinguished  from  separated  boundary  layer,  shear 
layer,  near  or  far  wake  flows).  The  material  is  divided  into  three  main  categories:  the 
inviscid  shock  layer,  the  boundary  layer,  and  the  fully  viscous  and  merged  layers.  These 
main  categories  are  sub-divided  under  the  headings  of  thermodynamic  equilibrium/ideal  gas 
and  nonequilibrium  gas,  and  the  inviscid  shock  layer  is  farther  divided  into  flows  which  are 
locally  subsonic  or  supersonic.  For  the  consideration  of  the  reader  who  is  not  completely 
familiar  with  the  reentry  flow  field  problem,  an  introduction  is  presented  which  includes  a 
discussion  of  the  problem  and  definitions  for  many  of  the  specialized  terms. 

3.1.1  INTRODUCTION 

In  the  reentry  body  flow  field  problem  the  word  "reentry"  is  quite  generally  (and  herein)  taken 
to  imply  that  the  velocity  of  a  body  as  it  begins  its  descent  into  a  planet's  atmosphere  is  less 
than  the  velocity  required  to  escape  the  planet's  gravitational  attraction.  The  escape  velocity 
from  earth  is  approximately  37, 000  feet/second,  and  typical  reentry  velocities  will  range 
from  15, 000  to  25,000  feet/second.  This  aspect  is  important  in  that,  for  reentry  velocities, 
the  energy  losses  within  the  flow  field  due  to  radiation  are  of  a  secondary  nature  and,  as 
pointed  out  by  Grusyczynsld  and  Warren1,  can  be  neglected;  however,  for  entry  velocities, 
radiation  effects  must  be  taken  into  acoount. 
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Another  aspect  of  the  problem  which  should  be  clarified  is  that  the  words  "body  flow  field" 
denote  the  flow  field  about  the  forebody  of  the  reentry  body  as  distinguished  from  the  flow 
field  in  the  base  region  or  the  near  or  far  wake  regions. 

The  flow  field  environment  which  developes  about  a  vehicle  during  atmospheric  reentry  has 

2 

been  divided  into  seven  regimes  by  Probstein  in  an  analysis  based  on  the  stagnation  region 
of  an  approximately  spherical  body.  With  A  denoting  the  ambient  mean  free  path,  R^  the 
body  radius,  and  e  the  ratio  of  the  free  stream  density  to  the  average  density  in  the  stagna¬ 
tion  region,  the  four  regimes  within  the  continuum  are  defined  by  Probstein  as  follows: 

Fully  Merged  Layer  Regime 

R. /10  <  A  SR 

D  oo  D 

Incipient  Merged  Layer  Regime 
e  ^  R./10  <  A  SR./10 

b  oo  D 

Viscous  Layer  Regime 


R./10  <  A  s  c*  R. /10 

b  oo  D 


Boundary  Layer  -  Inviscid  Layer  Regime 


A  s  e  R, /10 

oo  D 
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For  typical  earth  reentry,  the  average  density  ratio  e  will  be  about  1/10.  For  a  nose  radius 
of  1/2  inch  and  taking  A  from  the  1952  U.S.  Standard  Atmosphere  Tables,  the  above  four 

oo 

regimes  in  terms  of  altitude  (kilofeet)  are: 


235  <  Fully  Merged  Layer  £  285 
210  <  Incipient  Merged  Layer  £  235 
175  <  Viscous  Layer  £210 
Boundary  Layer  -  Tnviscid  Layer  £  175 

Under  the  same  conditions  except  for  a  nose  radius  of  one  foot,  the  above  regimes  will  shift 
upwards  by  about  70  kilofeet. 


In  both  of  the  merged  layer  regimes,  the  shock  wave*  which  envelops  the  reentry  body  is  not 

a  discontinuity  but  as  Probstein  points  out  is  relatively  thick,  ranging  from  a  value  of  about 

R  /3  a*  the  upper  limit  of  the  fully  merged  layer  regime  to  about  R.  /90  at  the  lower  limit  of 
b  D 

the  incipient  merged  layer  regime.  Throughout  both  of  these  regimes  transport  effects 

(thermal  conductivity,  viscosity  and  difusion)  are  of  primary  importance,  and  must  be 

included  in  the  analysis  of  both  the  shockwave  and  the  intervening  shock  layer  between  the 

back  face  of  the  shock  wave  and  the  body  surface.  For  both  of  the  merged  layer  regimes, 

Probstein  recommends  the  solution  of  the  Navier-Stokes  equations  (along  with  the  continuity, 

energy,  and  state  equations)  including  transport  effects,  with  the  boundary  conditions  deter- 

199 

mined  by  the  body  and  free  stream.  However,  Cheng  has  pointed  out  that  when  the  shock 
wave  thickness  can  be  neglected  compared  to  the  body  radius  (which  on  a  1/10  basis  includes 
all  of  the  incipient  merged  layer  regime  and  about  70  percent  of  the  fully  merged  layer  regime), 
the  shockwave  thickness  effects  are  secondary  compared  to  the  substantial  changes**  in  the 


*  The  shockwave  is  sometimes  referred  to  at  the  shock  transition  zone,  particularly  when 
it  is  not  of  infinitesimal  thickness. 

**  For  flow  across  an  oblique  shock  wave  of  infinitesimal  thickness  bounding  a  Fow  field 
region  in  which  transport  effects  (viscosity,  conductivity  and  diffusion)  can  be  neglected, 
i.e. ,  for  inviscid  flow,  the  velocity  component  tangent  to  the  shockwave  and  the  stagnation 
enthalpy  (total  energy)  remain  unchanged  across  the  shockwave.  The  value  of  this  velocity 
component  will  depend  on  the  shock  wave  angle,  whereas  the  value  of  the  stagnation  enthalpy 
does  not.  The  equations  governing  this  type  of  flow  can  be  found  in  Shapiro3  and  are  gen¬ 
erally  referred  to  as  the  Ranldne-Hugoniot  relations. 


3-3 


tangential  velocity  component  and  the  stagnation  enthalpy  across  the  shock  wave  caused  by 
the  transport  processes.  Cheng,  therefore,  maintains  that  the  shockwave  can  be  treated  as 
a  discontinuity  whenever  its  thickness  is  negligible  compared  to  the  body  radius  and  the  flow 
across  it  adequately  described  by  the  Rankine-Hugoniot  relations  after  they  have  been  modified 
to  account  for  transport  effects.  Thus,  over  a  substantial  part  of  the  merged  layer  regimes, 
the  modified  Rankine-Hugoniot  relations  can  be  used  to  provide  outer  boundary  conditions  for 
the  shock  layer  solution. 

In  the  viscous  layer  regime,  transport  processes  are  of  primary  importance  within  the  entire 
shock  layer,  but  the  shock v-ave  is  thin  enough  to  be  considered  a  discontinuity  and  the  Rankine 
Hugoniot  relations  provide  the  properties  behind  it  to  be  used  as  boundary  conditions  for  the 
shock  layer  solution. 

In  the  boundary  layer  -  inviscid  layer  regime,  the  shockwave  thickness  is  on  the  order  of 
several  ambient  mean  free  paths  and  is,  therefore,  a  real  physical  discontinuity.  Transport 
effects  a-'e  important  in  the  region  of  the  shock  layer  adjacent  to  the  body  surface  (the 
boundary  layer),  whereas  they  are  ngeligible  in  the  shock  layer  region  adjacent  to  the  shock 
wave  (the  inviscid  layer*).  In  this  case  the  reentry  body  flow  field  problem  cp  ">s  solved 
separately  in  each  of  these  regions  by  matching  the  flow  properties  at  the  interface,  i. e. ,  the 
boundary  layer  edge.  This  matching  process  is  aone  as  follows.  An  inviscid  layer  solution 
is  obtained  as  if  there  existed  no  boundary  layer  on  the  body,  i.  e. ,  the  body  surface  is  treated 
as  the  zero  streamline.  The  distribution  of  flow  properties  along  the  Inviscid  zero  stream¬ 
line  are  then  used  to  prrvide  boundary  layer  edge  conditions  and  a  boundary  layer  solution  is 
obtained  which  determine  the  actual  mass  flow  within  the  boundary  layer.  By  matching  this 
mass  flow  to  that  contained  within  the  appropriate  streamline  in  the  inviscid  solution,  a  new 
distribution  of  properties  for  the  boundary  layer  edge  conditions  is  determined.  If  this 
distribution  of  properties  differs  significantly  from  the  first  distribution,  the  boundary  layer 
must  be  recalculated  using  the  new  edge  conditions  until  satisfactory  agreement  in  the  edge 

*  Because  the  free  stream  is  isoenergetic  (constant  stagnation  enthalpy)  and  the  flow  across 
the  shock  wava  is  adiabatic,  the  inviscid  layer  is  also  isoenergetic. 
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conditions  from  one  iteration  to  the  next  _s  obtained.  (When  this  iteration  process  is  required 
to  match  edge  conditions,  it  is  referred  to  as  the  vorticity  interaction  problem.)  After  the 
edge  conditions  are  matched  satisfactorily,  the  boundary  layer  solution  is  used  to  calculate 
the  boundary  layer  displacement  thickness,  which  is  the  distance  the  inviscid  layer  is 
displaced  away  from  the  body  surface  by  the  boundary  layer.  If  the  displacement  thickness 
can  be  neglected  compared  to  the  shock  layer  thickness,  then  the  matched  inviscid  layer  - 
boundary  layer  solution  has  been  obtained.  However,  if  it  cannot  be  neglected  (which  is 
referred  to  as  the  displacement  thickness  problem),  then  the  displacement  thickness  must  be 
added  to  the  body  surface  and  the  inviscid  shock  layer  recalculated  treating  the  displaced 
"body"  surface  as  the  zero  streamline.  Starting  with  the  boundary  layer  in  which  the  edge 
conditions  were  matched  for  the  original  body,  the  matching  process  descr'Uxl  above  must 
be  repeated  for  the  displaced  "body. " 

Due  to  the  nature  of  the  boundary  layer  problem,  the  coordinate  bystem  selected  for  the 
solution  is  always  located  on  the  body  surface.  In  the  axlsymmetric  case  (as  a  simple 
example),  the  boundary  layer  equations  in  the  surface  coordinate  system  will  involve  both  the 
transverse  radius  of  curvature  (the  radial,  or  normal,  distance  from  &  local  point  in  the 
boundary  layer  to  the  axis  of  symmetry)  and  the  longitudinal  body  radius  of  curvature  (the 
body  surface  radius  of  curvature  in  a  meridinal  plane).  If  the  boundary  layer  thickness  is 
small  compared  to  the  transverse  radius  of  curvature,  the  transverse  radius  of  curvature 
can  be  approximated  by  the  body  radius,  i.e. ,  the  radial  distance  from  the  body  surface  to 
the  axis  of  symmetry. 

If  in  addition  the  transverse  curvature  effect  is  neglected,  the  governing  equations*  constitute 
what  is  known  as  first  -  order  boundary  layer  theory.  If  either  transverse  or  longitudinal 
curvatuit-  effects  are  taken  into  account,  it  is  referred  to  as  second  -  order  boundary  layer 
theory. 

*  Within  the  boundary  layer,  neglecting  the  longitudinal  curvature  effect  is  equivalent 
to  neglecting  the  transverse  pressure  gradient  as  shown,  for  example,  in  Cheng199. 
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To  make  the  boundary  layer  problem  more  tractable,  the  longitudinal  and  traverse  coordin¬ 
ates  are  usually  recast  to  simplify  the  form  of  the  governing  equations.  For  compressible 
flows  the  Howarth  -  Dorodnitsyn  transformation  is  often  applied,  which  reduces  the  equations 
to  the  incompressible  form  by  setting  the  transformed  transverse  coordindate  proportional  to 
the  integral  of  the  gas  density  c'rer  the  physical  transverse  coordinate.  In  general,  the 
boundary  layer  edge  and  wall  (body  surface)  conditions  will  vary  with  the  longitudinal  coordin¬ 
ate.  In  some  instances  it  is  possible  to  recase  the  longitudinal  coordinate  and  the  dependent 
variables  so  that  the  transformed  boundary  conditions  are  constant  and  the  wall  and  edge 
conditions  appear  in  the  governing  conditions  grouped  as  terms  which  are  constant  or  which 
vary  slowly  with  the  transformed  longitudinal  coordinate.  In  the  case  when  these  terms  are 
constant,  it  is  called  a  similarity  transformation,  and  the  problem  reduces  to  a  two  point 
boundary  value  problem  because  the  dependent  variables  are  functions  only  of  the  transverse 
coordinate.  In  the  cne  when  these  terms  vary  slowly  with  the  transformed  longitudinal 
coordinate  (so  that  they  may  be  considered  constant  at  any  fixed  value  of  the  coordinate), 
the  boundary  layer  problem  again  reduces  to  a  two  point  boundary  value  problem,  and  the 
transformation  is  referred  to  as  locally  similar. 

For  the  inviscid  shock  layer  solution,  the  origin  of  the  coordinate  system  is  usually  located 
at  the  apex  or  nose  of  the  shockwave  or  nose  of  the  shockwave  or  body,  but  is  not  necessarily 
made  to  coincide  with  the  shock  wave  or  body  surface.  If  the  reentry  body  has  a  sharp  nose, 
the  shock  wave  will  be  attached  at  the  apex  and  the  flow  within  the  shock  layer  will  be  locally 
supersonic.  If  the  reentry  body  has  a  blunt  nose,  the  flow  within  the  shock  layer  in  the  nose 
region  is  locally  subsonic  and  becomes  locally  supersonic  as  it  expands  around  the  body. 

This  fact  complicates  the  blunt  body  problem  because  the  governing  partial  differential 
equations  are  elliptic  in  the  subsonic  region,  parabolic  on  the  sonic  surface  and  hyperbolic 
in  the  supersonic  region.  Where  the  flow  is  locally  supersonic,  the  well  established  method 
of  characteristics  can  be  used.  Where  the  flow  is  loyally  subsonic,  the  inviscid  shock  layer 
solution  must  be  obtained  by  other  methods  such  as  finite  difference  methods,  the  method  of 
integral  relations,  or  Taylor  series  or  other  expansions.  These  methods  are  classified  as 
either  inverse  or  direct.  In  the  inverse  methods,  the  coordinate  system  is  made  to  coincide 
with  the  shockwave,  the  contour  of  which  is  assumed.  The  blunt  body  problem  is  then  solved 
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as  an  initial  value  problem,  proceeding  away  from  the  specified  shock  wave  contour  until  the 
body  surface  is  determined  by  mass  flow  conservation  requirements.  By  judiciously  selecting 
the  shockwave  contour,  the  resulting  body  shape  can  be  made  to  closely  approximate  the 
desired  one.  However,  because  the  governing  equations  are  elliptic,  the  initial  value  problem 
is  not  well  posed  and  numerical  techniques  for  its  solution  are  inherently  unstable.  Although 
this  presented  serious  difficulties  in  some  of  the  first  applications  of  the  method,  it  has  to  a 
large  extent  been  overcome  by  better  integration  techniques,  coordinate  transformations  and 
other  improvements.  In  the  direct  method  the  coordinate  system  is  located  on  the  body,  and 
the  inviscid  blunt  body  problem  is  solved  as  a  boundary  value  problem  whether  the  posi¬ 
tion  of  one  of  the  boundaries  is  unknown.  The  location  and  shape  of  the  shock  wave  is ,  there¬ 
fore,  determined  in  the  process. 

A  relatively  new  and  very  promising  method  of  solving  the  steady  reentry  body  flow  field 

4 

problem  is  the  asymptotic  time  method  proposed  by  Crocco  in  which  the  steady  flow  field 
is  recognized  to  be  the  state  attained  through  any  transient  process  after  infinite  time.  The 
inclusion  of  the  transient  terms  in  the  governing  equations  makes  them  hyperbolic  and  there¬ 
fore  amenable  to  solution  as  an  initial  value  rather  than  as  a  boundary  value  problem.  In 

theory  the  method  is  capable  of  solving  the  complete  viscous  shock  layer  problem  as  has 

5  6 

been  demonstrated  by  Scala  and  Gordon. 

At  reentry  velocities,  the  thermal  environment  within  the  shock  layer  is  so  high  that  the 
internal  energy  modes  of  the  gas  molecules  (rotational,  vibrational, and  electronic  states) 
become  excited  and  chemical  reactions,  including  dissociation  and  ionization,  occur.  In  this 
case  the  shock  layer  properties  cannot  be  accurately  prescribed  by  treating  the  gas  as  ideal, 
i.  e. ,  as  both  thermally  perfect  (for  which  the  pressure  is  proportional  to  the  product  of  the 
ga«  density  and  temperature)  and  calorically  perfect  (for  which  the  specific  heat  capacities 
are  constant).  At  the  lower  altitudes  (below  approximately  100  kilofeet  for  typical  earth 
reentry),  this  problem  is  rather  easily  resolved  because  the  gas  is  in  local  thermodynamic 
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equilibrium.*  At  the  higher  altitudes  of  the  continuum  regime,  however,  this  problem  is 
considerably  more  difficult  because  the  gas  is  not  in  local  thermal  or  chemical  equilibrium 
throughout  a  substantial  part  of  the  flow  field.  In  this  case, to  accurately  describe  the  flow 
field  properties,  the  finite  rate  processes  for  both  thermal  nonequilibrium  and  chemical 
nonequilibrium  and  their  coupled  effects  must  be  taken  into  account. 


3.1.2  REVIEW 

In  the  following  review  the  attempt  has  been  made  to  identify  each  reference  with  enough 
clarity  to  enable  the  reader  to  determine  whether  or  not  the  referenced  work  is  pertinent  to 
his  particular  interest.  In  most  cases  the  method  of  solution,  salient  characteristics,  and 

results  are  indicated.  In  certain  cases  where  (in  die  author's  opinion)  a  substantial  contri- 

* 

button  has  been  made  which  might  not  be  extensively  known,  the  work  is  discussed  in  more 
detail.  It  is  by  this  aspect  alone  that  the  review  is  intended  to  be  a  critical  one. 


In  the  interest  of  brevity  the  following  abbreviations  have  been  employed: 


MOC 

A/2D 

3D 

E/IG 

Le 

Pr 


method  of  characteristics 
axisymmetric  or  two-dimensional 
three-dimensional 
equilibrium  or  ideal  gas 
Lewis  number 
Prandtl  number 


*  The  word  "local"  is  used  in  the  macroscopic  sense  to  denote  a  fluid  element  which  con¬ 
tains  enough  gas  particles  to  constitute  an  ensemble,  i.e. ,  its  average  properties  are  the 
most  probable  ones  determined  by  statistical  analysis  of  all  possible  quantum  states  of  the 
gas.  Thermodynamic  equilibrium  is  defined  to  include  mechanical,  thermal  and  chemical 
equilibrium.  For  mechanical  equilibrium,  which  is  defined  in  the  D'Alembert  sense  treating 
the  fluid  mass  acceleration  as  an  inertial  body  force,  the  sum  of  forces  acting  on  a  fluid 
element  is  zero.  In  thermal  equilibrium  the  energy  of  each  gas  particle  is  partitioned 
equally  per  degree  of  freedom  (equilabrated)  among  the  energy  modes.  In  chemical  equil¬ 
ibrium  the  relative  concentrations  of  the  component  species  of  the  gas  are  established  by  the 
equilibrium  of  all  possible  chemical  reactions.  The  concept  of  local  thermodynamic  (or 
mechanical,  thermal  or  chemical)  equilibrium  in  a  flowing  gas  is  applicable  only  when  the 
effects  of  property  gradients  throughout  a  fluid  element  (ensemble)  can  be  neglected.  It  is 
not  applicable  for  example  in  strong  shock  wave  flow  for  which  some  of  the  property  gradients 
tire  appreciable. 


3, 1,2.1  Iiiviscid  Shock  Laver 


3. 1,2. 1. 1  Equilibrium  or  Ideal  Gas 

3. 1.2.  l.l.l  Loyally  Supersonic  Flow.  The  conical  flow  problem  is  treated  in  detail  by 
Shapiro  .  For  an  ideal  diatomic  gas  tabulated  results  are  presented  by  Kopal8  and  Wang 
et  al.  and  graphical  results  by  Dailey  and  Wood19  and  the  Ames  Research  Staff11.  The 
method  was  extended  to  equilibrium  air  by  Johnson  whose  results  show  that  at  hypersonic 
speeds  the  conical  solution  is  substantially  affected  by  equilibrium  gas  properties  being 
characterized  by  a  thinner  shock  layer,  lower  pressure,  significantly  lower  temperature, 
higher  velocity, and  significantly  higher  density  than  for  an  ideal  diatomic  gas. 


The  A/2D  MOC  solution  for  rotational  flow  of  an  ideal  gas  is  treated  at  length  by  Ferri13 
and  Isenberg  and  Lin  .  The  method  has  been  extended  to  equilibrium  air  by  Gravalos  et 
al.  and  many  others  (e.  g. ,  Powers15,  Dresser16,  Inouye  et  al. 59).  Edsall17  improved 
the  standard  method  by  interpolating  for  the  entropy  of  a  streamline  in  a  table  constructed 
from  the  starting  data  and  calculated  shock  points,  thereby  eliminating  the  necessity  of  using 
the  projected  velocity  vector  to  linearly  interpolate  for  the  entropy  at  a  grid  point.  Gravalos 
et  al.  have  expanded  the  method  to  provide  for  bodies  with  compression  corners  which 
result  in  secondary  shock  waves.  Field  and  surface  properties  are  furnished  for  a  sphere- 
cone-cone  at  Mach  10  for  both  equilibrium  air  and  ideal  gas.  Davis19  has  developed  an 
A/2D,  E/IG  solution  which  provides  for  expansion  comers,  compression  corners, and  shock 
wave  intersections. 

Flow  over  an  axisymmetric  body  at  small  angle  of  attack  can  be  determined  by  a  limited 
application  of  the  3D  MOC.  The  method  is  described  in  detail  by  Ferri13,  20  and  Whitham21 
and  consists  of  solving  the  3D  characteristics  equations  by  a  perturbation  technique,  retain¬ 
ing  only  linear  terms  in  angle  of  attack  and  obtaining  the  necessary  zeroth  order  quantities 
from  the  axisymmetric  MOC  solution.  Ideal  gas  results  for  cones  are  presented  by  Stone22 
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■  i 


23  23  24  25 

and  Rakich  and  for  sphere  cones  by  Rakich  ,  Ferri  and  by  Bronp  and  Edeifelt  who 

26 

also  extended  the  method  to  equilibrium  air.  Rakich  broadened  his  previous  work  to 
equilibrium  air  and  presents  a  comparison  of  equilibrium  air  and  ideal  gas  results  for 
cones,  sphere  cones  and  ogives. 


27  28 

The  3D  MOC  solution  for  steady  rotational  flow  is  treated  by  Holt  and  Coburn  . 

Moretti  et  aL^  have  developed  and  applied  the  method  to  a  variety  of  3D  and  axlsymmetrlc 

bodies  at  angle  of  attack  for  both  equilibrium  air  and  ideal  gas.  Pridmore  Brown  and 
30 

Franks  have  extended  the  method  to  provide  for  secondary  shock  waves;  however,  no 
applications  are  presented.  C.  W.  Chu  has  developed  the  method  with  an  improved 
integration  technique.  Instead  of  satisfying  the  compatability  relation  exactly  alon^three 
arbitrarily  selected  bicharacteristics  from  the  base  of  the  Mach  cone  to  the  vertex,  Chu 
approximately  satisfies  it  for  all  (infinitely  many)  of  the  bicharacteristics  on  the  Mach  cone. 
He  does  this  by  m  inimizing  the  root  mean  square  of  the  residual  function  determined  from 
the  compatability  relation  which  requires  the  evaluation  of  three  contour  integrals  a^-ound 
the  Mach  cone  base.  3trom  -  has  extended  the  method  to  provide  *"*r  a  gas  which  is 
frozen  along  streamlines  but  may  havo  a  different  composition  on  each  streamline. 


Thommen  and  D'Attore33  have  shown  that  the  finite  difference  scheme  of  Lax  and  Wendoroff 
can  be  applied  successfully  to  the  3D  supersonic  inviscid  flow  problem. 


3. 1. 2. 1. 1.2  Locally  Subsonic  Flow.  The  blunt  body  problem  was  first  reviewed  by  Van 

Dyke34,  who  also  developed  an  inverse  numerical  solution,  and  later  extensively  by  Hays 

and  Prcbsiein  and  Swigert36,  45 .  The  work  which  was  covered  in  these  reviews  will  not  be 

35 

repeated  here,  and  the  reader  is  referred  to  Hays  and  Probstein  for  reviews  on  references 

36  45 

34  and  37  to  44,  and  to  Swigert  ’  for  reviews  on  references  45  to  53. 


Maslen54  and  Cheng  and  Gaitatzcs  55  have  developed  inverse  blunt  body  solutions  using  the 

thin  shock  layer  approximation  which,  considering  their  simplicity,  show  excellent  agreement 

56 

with  results  of  exact  solutions.  Van  Tuyl  has  developed  an  approximate  blunt  body  solution 

using  Pade  fractions  (ratios  of  polynomials)  to  approximate  the  terms  in  the  double  power 

series  expansion  method  and  shows  good  agreement  with  exact  solutions  for  spherical  noses. 

57 

Kao  has  developed  a  direct  series  truncation  blunt  body  solution  for  the  zero  angle  case 

34 

which  is  analogous  to  Swigert's  method  for  angle  of  attack. 

58 

Lomax  and  Inouye  have  solved  the  inverse  blunt  body  problem  by  numerical  integration  in 

physical  coordinates  using  a  predictor  -  co  -rector  and  a  smoothing  technique.  Extensive 

tables  for  spheres  in  eqi’.i librium  air  which  include  surface  and  field  properties  are  presented 

at  altitudes  between  J.OO  to  300  kft  and  free  stream  velocities  between  10  and  45  kft  /sec 

57 

(a  best  buy).  Inouye  et  al*  have  coupled  the  solution  of  Lomax  and  Inouye  with  the  method 
of  characteristics.  Ideal  gas  results  for  the  surface  and  shock  of  a  sphere-cone,  an  ellipsoid 

ft 

cylinder,and  a  sphere-cone-cylinder-flare  are  presented.  Dresser  and  Anderson  have 

developed  a  numerical  inverse  axisymmetric  solution  and  furnish  results  for  a  aphere-cone 
6 1 

cylinder.  Joss  has  broadened  the  numerical  Inverse  method  to  include  axisymmetric 

62 

bodies  at  angle  of  attack  for  itleal  gas.  Webb  et  al.  have  produced  a  finite  difference  inverse 

61 

solution  similar  to  Joss  and  furnish  resuits  for  the  Apollo  body. 


The  axisymmetric  direct  blunt  body  problem  using  Dorodnitsyn's  method  of  integral  relations 

has  been  solved  by  Traugott  3’  64  for  an  ideal  gas  and  by  Feldman65  and  Kuby  et  al.64  for 

47  67 

equilibrium  air  and  extended  to  angle  of  attack  by  Waldman  ,  Vaglio-Laurin  and  Leigh  and 
68 

Rothrock  .  However,  Leigh  and  Rothrock  found  that  when  the  entropy  differential  equation 
is  properly  treated  the  method  of  results  in  an  overdetermined  system  of  equations  evidently 
caused  by  the  boundary  conditions  of  the  shock.  Gravalos  et  al.  have  developed  an 
axisymmetric  direct  finite  difference  method  for  the  blunt  body  problem  in  which  the  body 
pressure  distribution  and  shock  wave  are  initially  prescribed  and  the  governing  equations  in 
intrinsic  coordinates  are  solved.  Based  on  the  results  of  the  solution, the  shock  wave  and 
surface  pressure  are  modified  and  the  solution  is  repeated  until  the  boundary  conditions  are 
satisfied  to  within  a  specified  tolerance.  Waiter  and  Anderson  have  also  produced  this 
type  of  solution. 


The  direct  blunt  body  problem  has  also  been  solved  using  the  asymptotic  time  approach. 

48  49 

Bohuchesky  and  Rubin  and  Bohachesky  and  Mates  have  developed  this  type  of  method 

in  which  the  governing  equations  are  solved  by  finite  difference , with  the  free  stream  and 

body  as  initial  conditions  allowing  the  shock  wave  to  develop  and  present  ideal  gas  results 

71  72 

for  the  Apollo  body  at  angle  of  attack.  Morettl  et  al.  *  have  achieved  a  significant 
improvement  in  the  finite  difference  asymptotic  time  method  by  treating  the  shock  wave  as 

a  moving  discontinuity,  and  furnish  ideal  gas  results  for  different  configurations.  Webb  and 

74  75  76 

Dresser  have  produced  a  solution  similar  to  that  of  Moretti's.  Sauerwein  ’  and 

77 

Boerleke  and  Brong  have  applied  the  asymptotic  time  approach  to  the  blunt  body  problem 
using  the  unsteady  three-dimensional  method  of  characteristics.  However,  this  method 
appears  to  be  considerably  less  efficient  and  accurate  than  that  of  Moretti's  method. 

3.1. 2.1. 2  Nonequilibrium  Gas 

3. 1. 2. 1. 2. 1  Locally  Supersonic  Flow.  The  two-dimensional  axisymmetric  method  of 

78  79 

characteristics  was  first  extended  to  a  nonequilibrium  gas  by  Wood  and  Kirkwood  ’  and 
30 

B.  T.  Chu  who  showed  that,  for  gas  in  thermal  or  chemical  nonequilibrium,  the 

characteristics  waves  are  propagated  at  the  local  frozen  acoustic  speed.  Sidney  et  al 
Hi.  82  83  84 

*  *  *  have  applied  the  nonequilibrium  method  of  characteristics  to  wedges  and 

cones  for  the  vibrational  relaxation  of  an  ideal  diatomic  gas  using  the  harmonic  oscillator 
as  a  model.  Their  results  show  that  vibrational  nonequilibrium  causes  curvature  in  the 
shock  wave  and  the  surface  pressure  overexpands  from  the  frozen  value  at  the  cone  tip  (which 
is  larger  than  the  equilibrium  value  because  of  the  larger  shock  angle  at  the  tip)  to  below 
the  equilibrium  value  and  asymptotically  approaches  from  below, whereas  the  gas  temperature 

gr  gg 

asymptotically  approaches  the  equilibrium  value  from  above.  Spurk  et  al.  ’  have 

developed  a  thermal  equilibrium  chemical  nonequilibrium  characteristics  solution  for  a  5 

species  -  10  reactions  air  system  and  furnish  results  for  the  distribution  of  surface 

87  88  hQ 

properties  and  the  shock  wave  shape  for  both  cones  and  wedges.  Gravalos  ’  ’  et  al. 

have  developed  a  thermal  equilibrium  chemical  uonequilibrium  characteristics  solution  for 
a  7  species  -  7  reactions  air  system  which  provides  for  an  expansion  corner  on  the 
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body.  Results  aro  presented  (in  reference  87)  for  a  45°  cone  with  a  39°  expansion  corner. 

Wood  et  al.9°  have  developed  a  characteristics  solution  with  both  vibrational  nonequilibrium 
and  chemical  nonequilibrium  for  a  7  species  -  7  reaction  air  system.  The  harmonic 
oscillator  model  Vi  used  for  the  vibrational  relaxation  of  molecular  oxygen  and  nitrogen  and 
the  effect  of  vibration  on  their  nonequilibrium  dissociation  rates  are  taken  into  account. 

Thermal  equilibrium  results  are  presented  for  a  sphere  -  cone,  and  both  vibrational  and 
chemical  nonequilibrium  results  for  an  ogive.  Their  sphere-cone  results  also  show  that  the 
nonequilibrium  streai.itube  gives  excellent  results  for  the  electron  density  distributions. 

Kliegel  et  al.91  have  developed  a  characteristics  solution  for  chemical  nonequilibrium 
(6  species-6  reactions  air  system)  with  the  provision  for  treating  vibrational  relaxation 
(without  vibrational  -  dissociation  coupling)  using  the  anharmonic  oscillator  model.  Thermal 
equilibrium  results  are  presented  for  the  surface  profiles  of  a  cone-cyunder.  South 
and  South  and  Newman94  have  applied  Dorodnitsyn's  method  of  integral  relations  to  flow 
over  pointed  bodies  (locally  supersonic  flow)  of  a  vibrationally  relaxing  ideal  diatomic  gas 
and  found  that  the  method  did  not  give  the  proper  asymptotic  behavior  because  of  the  improper 
weighting  of  the  frozen  flow  boundary  condition  at  the  shock  wave.  This  difficulty  was  overcome 
(in  reference  94)  by  weighting  the  shock  wave  boundary  condition  and  their  results  for  a^two 
strip  approximation  agree  well  with  the  characteristics  solution  of  Sedney  and  Gerber  . 

3. 1. 2. 1. 2. 2  Locally  Subsonic  Flow.  The  blunt  body  problem  for  a  nonequilibrium  gas  was 

first  solved  by  Lick99,96  using  a  numerical  inverse  method  for  the  subsonic  region  and  the 

method  of  characteristics  for  the  supersonic  region.  His  results  are  for  thermal  equilibrium 

and  although  necessarily  based  on  poor  reaction  rates  show  that  nonequilibrium  effects  do  not 
°  97 

significantly  affect  the  static  pressure  in  blunt  bocfy  flows.  Hall  et  al.  have  also  used  a 

numerical  inverse  method  to  solve  the  blunt  body  problem  for  a  thermal  equilibrium  chemical 

nonequilibrium  7  species  -  7  reactions  air  system.  Results  are  furnished  for  body  surface  and 

field  properties  for  a  spherical  nose  and  binary  scaling  is  demonstrated.  Additional  results 

98,99,  1 00 

whicn  have  been  obtained  from  this  solution  are  presented  by  Gibson  and  Marrone 

mi  i  c\9  .  „  103,104,105 

and  by  Wurster  and  Marrone101’  .  Based  on  the  work  of  Treanor  and  Marrone 

Marrone106  expanded  the  solution  of  Hall  et  al.97  to  include  vibrational  nonequilibrium  which 
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treats  both  coupled  vibration  -  dissociation  and  dissociation  -  vibration  effects.  The 

vibrational  energy  levels  are  prescribed  from  spectroscopic  data  and  "preferred" 

dissociation  from  the  upper  vibrational  levels  is  allowed  by  an  arbitrary  exponential 

probability.  A  Boltzman  distribution  is  assumed  for  the  vibrational  level  populations. 

Vibrational  chemical  nonequilibrium  results  are  presented  for  the  hypersonic  flow  of  oxygen 

107 

over  a  spherical  nose.  Lee  and  Chu  have  produced  a  numerical  inverse  blunt  body  solution 

for  thermal  equilibrium  chemical  nonequilibrium  (5  species  -  6  reactions  air  system)  and 

97  108  109 

furnish  a  detailed  comparison  of  their  results  with  those  of  Hall  et  al.  Shih  et  al.  ’  , 

have  applied  Dorodnitsyn's  method  of  integral  relations  to  the  direct  blunt  body  problem  for 
thermal  equilibrium  chemical  nonequilibrium  (5  species  -  6  reactions  air  system).  Results 
using  the  one  strip  approximation  are  furnished  and  frozen  and  equilibrium  air  results 
compared  for  a  sphere. 

The  nonequilibrium  stream  tube  method  is  an  approximation  technique  whict  can  be  applied  to 

either  locally  subsonic  or  supersonic  flows.  Thermal  equilibrium  chemicf  1  nonequilibrium 

streamtube  solutions  for  air  have  been  developed  by  Bloom  et  al.  *10,  112 ,  Lin  and  Teare 

Eschenroeder  et  al. l14, 115, 116  and  Lord!  et  al. 117  and  superimposed  on  Ideal  gas  or 

118  119 

equilibrium  air  flow  field  solutions  by  McMenamin  and  O'Brien  and  Karydas  .  Chen  and 
120 

Eschenroeder  have  developed  a  streamtube  solution  which  provides  for  vibrational 

121 

relaxation  using  the  harmonic  oscillator  model.  Dresser  et  al.  have  produced  a  streamtube 
solution  which  allows  for  both  vibrational  and  electronic  relaxation. 


3. 1.2.2 


3. 1. 2. 2. 1  Equilibrium/Ideal  Gas 


3. 1. 2. 2. 1. 1  First  Order,  The  similar  boundary  layer  solution  for  this  category  is  treated 

122 

by  Hayes  and  Probstein  and  the  stagnation  point  boundary  layer  is  an  important  case  of  its 
123 

application.  Libby  has  developed  a  low  temperature  air  solution  with  air  injection  or  suction. 
He  solves  the  two  point  boundary  value  problem  after  a  transformation  to  the  Crocco  variables 
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numerically  by  the  "shooting  through17  or  initial  value  method  determining  the  necessary  initial 

slopes  at  the  edge  by  means  of  an  asymptotic  expansion  which  is  valid  at  the  edge.  Tabulated 

and  graphical  results  are  presented  for  a  wide  range  of  wall  conditions  and  a  comparison  with 

an  approximate  integral  method  is  made,  Weston1  ^  has  produced  a  high  temperature 

equilibrium  air  solution  with  air  injection  or  suction.  He  uses  the  shooting  through  method  to 

oltain  the  numericxd  solution  by  assuming  velocity  and  enthalpy  slopes  at  the  wall  and  iterating 

until  the  boundary  conditions  at  the  edge  are  satisfied.  Tabulated  profiles  for  an  axisymmetric 

stagnation  point  are  provided  for  a  variety  of  wall  conditions,  Axisymmetric  stagnation  point 

125 

boundary  layers  for  air  have  also  been  developed  by  Herring  (without  mass  transfer) 

126  127  128 

and  by  Hoshezaki  and  Smith  and  Howe  and  Mersman  (with  mass  transfer).  Strom 

has  developed  a  three-dimensional  stagnation  point  finite  difference  solution  with  mass  transfer. 

Axisymmetric  similar  solutions  for  other  than  the  stagnation  point  have  been  produced  by 

Waiter  and  Anderson129  and  Grabow130.  Brant  and  Burke  have  developed  an  axisymmetric 

locally  similar  solution  for  an  arbitrary  equilibrium  gas  which  may  contain  many  species. 

Locally  similar  three-dimensional  (with  linearized  secondary  flow)  boundary  layers  have  been 

132  133 

developed  by  Beckwith  and  Beckwith  and  Cohen  for  equilibrium  air  or  ideal  gas  with 

no  mass  transfer,  by  Polak  and  Li4^4  for  an  ideal  gas  with  no  mass  transfer,  and  by 
135,  136 


Hang  et  al. 


with  mass  addition  for  frozen  chemistry. 


137 

Integral  solutions  for  the  axisymmetric  boundary  layer  have  been  developed  by  Moran  and 
138  139 

Carlson  for  a  binary  gas  and  by  Liu  and  Kuby  for  the  combustion  of  a  graphite  surface 


with  frozen  chemistry  in  the  boundary  layer.  Kang' 
integral  solution  for  ideal  gas  with  mass  transfer. 


140 


has  produced  a  three-dimensional 


Li141  has  solved  the  axisymmetric  boundary  layer  in  the  stagnation  region  by  a  series 

expansion  method  in  which, after  a  transformation  to  the  similarity  variables,  the  velocity, 

enthalpy  and  stream  function  are  expanded  in  a  power  series  for  which  the  coefficients  are 

142  143  144 

functions  only  of  the  transverse  coordinate.  Parr  ,  Krause  ’  and  Waiter  and 
145 

Anderson  have  developed  A/2D,  E/IG  boundary  layer  solutions  based  on  an  implicit  finite 

146 

difference  method  develo  ?ed  by  Flugge-Lotz  and  Blottner  .  In  this  method  the  momentum 
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equation  is  linearized  and  the  boundary  layer  is  solved  step  wise  in  the  longitudinal  coordinate 

147 

by  inverting  a  tri-diagonal  matrix  at  each  step.  Flannelop  has  expanded  the  method  to 

148  149 

three-dimensional  boundary  layers.  Raetz  and  Der  and  Raetz  have  also  developed 
numerical  solutions  for  the  three-dimensional  boundary  layer. 


The  turbulent  boundary  layer  problem  has  been  treated  by  Crocco  .  Denison  has 
produced  results  for  turbulent  boundary  layer  for  equilibrium  air  on  an  ablating  graphite 
surface  using  the  Crocco  solution  (total  enthalpy  and  species  cross  fractions  linear  function 
of  velocity)  with  the  Van  Driest  velocity  profile. 


3. 1. 2..2. 1. 2  Second  Order.  The  second-order  ax‘ symmetric/two-dimensional  boundary 

152  153  154 

layer  problem  has  been  treated  and  reviewed  by  Van  Dyke  *  *  .  Lewis  and 

co-workers155, 156  and  Adams157  have  applied  the  second  order  theory  of  Van  Dyke  using  the 

158 

numerical  method  developed  by  Davis  and  Flugge-Lotz  for  a  perfect  gas.  Smith  and 

co-workers159’160,161,162  have  solved  the  second  order  (transverse  curvature)  boundary 

layer  problem  as  a  sequence  of  two  point  boundary  value  problems  by  converting  the  governing 

equation  to  ordinary  equations  through  the  process  of  replacing  the  longitudinal  derivatives  by 

their  finite  difference  approximations.  The  two  point  boundary  value  problem  is  then  solved 

at  each  longitudinal  station  by  the  shooting  through  (or  initial  value)  technique.  This  method 

163 

has  been  applied  to  equilibrium  air  flow  on  a  sphere- cone  by  Lewis  and  Whitfield  and 
164 

Mayne  et  al. 1  for  a  binary  gas  on  a  sharp  cone.  Using  the  implicit  finite  difference  technique 

1  Jk  /J  ICC 

of  Flugge-Lotz  and  Blottner*  ,  Levine100  has  developed  an  A/2D,  E/IG  second  order  boundary 
layer  solution  which  includes  transverse  curvature  and  mass  addition. 

166  167 

Similar  second-order  solutions  have  been  developed  for  an  ideal  gas  by  Maslen  ,  Yasuhara 

1 68 

and  Probstein  and  Elliot  and  for  equilibrium  air  or  a  mixture  of  ideal  gases  (binary  diffusion) 
169 

by  Sagendorph  .  Maslen  includes  both  transverse  and  longitudinal  curvature,  whereas  the 
others  include  transverse  curvature. 
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3. 1.  2.  2.  2  Nonequilibrium  Gas 


3. 1.  2. 2. 2. 1  First  Order  Theory.  The  similar  nonequilibrium  gas  4./2D  boundary  layer 

170 

problem  is  treated  by  Dorrance  .  Finite  difference  solutions  have  been  developed  and 
applied  at  the  stagnation  point  qy  Fay  and  Riddell1  1  and  Fay  and  Kaye172.  Fay  -nd  Riddell 
used  the  shooting  through  or  initial  value  technique  and  present  catalytic  and  noncatalytic 
wall  results  for  air  with  P=0.  71,  a  Sutherland  viscosity  law  and  Le  =  1,  l.  4  and  2.  Fay  and 
Kaye  used  the  implicit  method  for  nitrogen  dissociation  treating  the  chemical  production 
term  as  a  parameter.  0*her  similar  solutions  have  been  developed  by  Moore  and  Pallone173 
for  air,  Libby  and  Pierucci  74  for  air  plus  hydrogen, and  by  Libby  and  Liu1'5.  Locally 
similar  solutions  for  trace  contaminants  in  air  have  been  obtained  by  Lenard176, 177. 


178 

Pallone  et  al.  have  applied  Doroditsyn's  method  of  integral  relations  to  the  nonequilibrium 
air  boundary  layer  problem.  Property  profiles  at  three  stations  on  an  8°  cone  for  three  free 
stream  conditions  and  a  comparison  of  the  effects  of  diffusion  is  made.  Smith  and  Jaffe179 
have  extended  their  previous  finite  difference  method  to  the  nonequilibrium  gas  bou.idaiy 
layer.  Results  are  presented  for  a  binary  dissociating  gas.  Blottner180,  181»  182  has 
developed  an  implicit  finite  difference  solution  for  the  nonequilibrium  gas  boundary  layer* 
based  on  the  work  of  Flugge-Lotz  and  Blottner.  His  solution  it,  capable  of  treating  catalytic 
(recombined)  or  noncatalytic  wall,  mass  transfer,  and  multicomponent  diffusion  for  arbitrary 
(25  species  -  40  reactions)  chemical  system.  Air  results  are  presented  for  a  sphere-cone, 
cone  ,and  hemisphere-cylinder  (binary  giis).  Additional  applications  of  Blottner's  solution  have 
been  presented  by  Lew  .  Moore  and  Lee  hrve  developed  a  solution  analogous  tr 
Blottner's. 

Approximate  solutions  for  the  nonequilibrium  gas  boundary  layer  have  been  produced  bv  Fox185 
186  J 
and  Haydav  .  Fox  solves  the  problem  in  terms  of  the  eigenfunctions  of  the  basic  boundary 

layer  operator  for  a  ideal  dissociating  gas.  A  particular  solution  is  required  for  the  species 

equation  and,  hence,  the  method  is  not  applicable  to  multicomponent-multireaction  system. 

Hayday  uses  an  asymptotic  series  expansion  method  in  which  the  coefficients  are  evaluated  by 

approximately  integrating  the  governing  equations.  Results  are  present  only  for  an  equilibrium 
air  stagnation  point. 


187  188 

Binary  scaling  laws  have  been  developed  by  Lee  and  Levlnsky  and  Levinsky  for  cones 

1  RQ 

in  pure  air  and  will  be  extended  to  cones  at  angle  of  attack  by  Lee  and  Baker  . 


3. 1.  2.  2.  2,  2  Second  Order  Theory.  No  publications  available  in  this  category. 


3. 1.  2.  2  *  is cous  Layer  -  Merged  Layer 


3. 1. 2.  3. 1  Equilibrium  or  Ideal  Gas 

190 

The  fully  viscous  layer  and  merged  layer  problems  are  treated  by  Hayes  and  Probsteln 
191 

Chen  et  al.  have  solved  the  viscous  Hyer  stagnation  region  problem  for  an  ideal  gas  with 

mass  transfer.  Modified  Rankine-Hugonoit  relations  and  a  linear  viscosity  are  employed. 

Results  are  presented  for  the  axisymmetrk  stagnation  point.  Goldberg  and  Goldberg  and 
193 

Scala  have  developed  a  series  expansion  solution  for  the  viscous  layer  stagnation  region 
for  equilibrium  air  with  mass  transfer.  A  Sutherland  viscosity  law  is  used  and  Lewis  and 
Prandtly  number  are  constant.  Results  are  presented  for  a  variety  of  wall  conditions  and 
Reynolds  numbers. 


195 

Kao  has  done  a  comparison  of  in  vise  id,  viscous  layer  and  third-order  boundary  layer 

theory  for  the  stagnation  region.  Results  are  presented  for  Reynolds  numbers  (based  on  shock 

1 95 

viscosity)  of  10,100  and  1, 000.  Probetein  and  Kemp  have  used  the  constant  density 
assumption  to  obtain  the  viscous  layer  stagnation  region  solution  for  an  ideal  gas.  Lenar d^ 
has  developed  a  merged  layer  stagnation  point  solution  for  an  ideal  gas.  Wei197  has  used 
the  method  of  inner  and  outer  expansions  to  obtain  the  viscous  layer  solution  on  a  slender 
body  for  an  ItKal  gas.  Waldron  has  used  a  perturbation  technique  to  obtain  the  viscous  layer 
solution  for  an  ideal  gas  on  a  slender  body.  Heat  transfer  results  are  presented. 


3. 1. 2.  3. 2  Nonequilibrium  Gas 

1 99 

The  viscous  layer  merged  layer  problem  for  a  nonequilibrium  gas  is  treated  by  Cheng 
Cheng  also  developed  a  numerical  solution  for  the  stagnation  region  and  presents  results  for 
a  binary  gas. 
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Shih  and  Krupp  have  solved  the  viscous  layer  stagnation  region  problem  for  an  ideal 
dissociating  gas  using  the  series  expansion  method  and  numerically  integrating  the  ordinary 
differential  equations. 

201 

Chung  has  developed  a  viscous  layer  stagnation  region  solution  for  a  binary  gas  and  Chung 
,  202  . 

et  al.  have  developed  a  merged  layer  st«^nation  region  solution  for  a  binary  gas.  Lee  and 
Zierten  have  developed  a  merged  layer  stagnation  region  solution  for  air  by  superimposing 
on  frozen  flow  field  as  calculated  by  Chung’s  method. 
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Table  3-3.  Summary  of  Geometric  Variables  for  Selected  Body  Surface  Normals 
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In  Figures  3-1  through  3-10  the  contrast  between  equilibrium  and  nonequilibrium  effects  at 
an  altitude  of  150  Idlofeet  is  demonstrated  by  plotting  tne  streamwise  variation  of  the  flew 
properties  against  the  distance  along  the  body  surface  for  three  streamlines  wlv  h  intersect 
the  bow  shock  at  angles  of  approximately  50,  70iand  89  degrees.  In  sequence,  these  plots 
show  this  variation  for  the  pressure,  temperature,  density,  velocity  component  parallel 
to  the  surface  (all  normalized  by  the  corresponding  free  stream  quantities),  and  the  number 
.entities  of  molecular  oxygen,  molecular  nitrogen,  atomic  nitrogen,  atomic  oxygen,  nitric 
oxide- and  electrons.  Before  proceeding  to  a  brief  discussion  of  some  of  the  more  interesting 
aspects  of  these  results,  it  should  be  pointed  out  that  in  Figure  3-1,  due  to  the  positive 
normal  pressure  gradient  in  the  nose  region,  the  higher  shock  angle  streamlines  initially 
expand  more  rapidly;  however,  in  the  over-expansion  region  this  trend  is  reversed  prior  to 
converging  at  the  cone  surface  pressure  value. 

For  the  nonequilibrium  results  shown  in  Figures  3-1  through  3-10,  it  may  be  noted  that  ihe 
flow  on  the  89  degree  streamline  is  almost  in  equilibrium  at  the  first  normal.  This  is  caused 
by  its  passage  through  the  near  stagnation  region.  At  this  point  oxygen  is  almost  entirely 
dissociated.  As  the  flow  continues  to  expand  over  the  nose,  atomic  oxygen  freezes  in  the 
nonequilibrium  case,  whereas  for  equilibrium  it  begins  to  recombine.  This  single  feature, 
more  than  any  other,  determines  the  nonequilibrium  effects  which  are  indicated  in  these 
results. 

In  Figures  3-11  through  3-25,  nose  bluntness  and  altitude  effects  on  the  nonequilibrium  flow 
properties  are  shown  by  plotting  the  flow  properties  against  shock  layer  thickness,  with 
altitude  as  a  parameter,  for  three  selected  normals  located  along  the  body  surface  at 
equal  to  0.49S2,  61.  25, and  245.3.  Sequentially  in  groups  of  three,  the  normal  profiles  are 
giver  for  the  pressure,  temperature,  density, and  velocity  parallel  to  the  surface  (all  normalized 
by  the  corresponding  free  stream  quantity)  and  the  electron  number  density.  Profiles  of  the 
remaining  species  number  densities  are  presented  in  Table  3-4,  3-5,  and  3-6  for  the  different 
altitudes  at  the  selected  normals.  In  Figure  3-26  the  integral  of  the  electron  density  across 
the  shock  layer  is  plotted  against  the  normalized  coordinate  along  the  body  surface  with 
altitude  as  a  parameter.  Some  of  the  more  noteworthy  aspects  of  the  results  are  discussed 
in  the  following  paragraphs.  3-53 
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Figure  3-1.  Stream  wise  Variation  of  Normalized  Pressure  vs.  Normalized  Coordinate  Along  Bod 
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Figure  3-3.  Streamwise  Variation  of  Normalized  Density  vs.  Normalized  Coordinate 

Along  Body  Surface 


mmmm  mmfeBiiiMiiiliHiElEini 

SSSSSISU!99£»-"-»W!S«|| 


mm*  Mtfm 


M*^|ilM||aHapptllll 

ii  mm  niiii  aili- -  i  iii 


HMHIil 

IBBIJBBI 


:::::::: 


ALTITUDE  =>  150  KILOFEET 

Rw  -  1.  27  CM 

N  0 
v„  =  6,  60O\  CM/SEC 


INDEX 

SYMBOL  SMOCK  ANGLE 


NORMALIZED  COORDINATE  ALONG  BODY  SURFACE, 


=“  "’EEiiii:5S^jTlMfiBiFijeigj 
i.  .'ti'iiw  !ii.ii*mnmtu«iiin»iii!|rw!"n 


mam 


lliwi- 


milWN  nrp  <MW\ 

r  y-  - 

S^Efer^iinj 

nr  n  ’«fcl 

prlm^S^  mihiW 

feaaJSt* 


ilk  fill  ■  -  -J  m  ’I « 


.H*wy- -u-  i+rnimv  SHJW 

Hiwprfiidiifiimwfl-fciiiii 

emhbsm 

MSSSg 


:LnlHIVW  ■  *(H  it  fr  II  »sn 

[•bp*****!**  M  i 


altitude  -  i&o  kilofeft 

■  1.  27  CM 


INDEX 

SVMBOL  SHOCK  ANGLE 
NONEQ  iDtO) 


«- . . . 

;  ii. 

•|.::  ;!;:;IM  1 115! H*P lUi^r  =5  V 

HMlipBili 

:uiuiiu  mi, ,  uniutfuKiui  rjmi  — ■ 

•  ,  .  '.t 

mm  IBMfrl»Sft«agliMMP 

mBiigiigiin 

iisr"  Si ".  :Ai\  !asB3tinPj 

wan  n  mi  Ltfiii'i  .in  MitttuiiL'himSim-'iU  lOrntwni 

hi  ii  i  v '  mi  r  pnwgBTi 

5*  •4<m  i»‘!wBSiai 


. . . 

HIl'IBlflillH! 


IHigiiSll! 


i^illk  fJl.litlH!  .1  i  SStfilfl!  IlH  liilllllllHlIWHl 

ill!  Ilii Oil  n  ■  ■  MM  M  IlfilfLu!  h= 

J13!«l!!ig!S';i=HEi5ia!!  Hi! 

N(iill|l-'ll|Li  fLidJ  ii  1 

[niii  jlpilft  3  lfN!llll!ii:Hn||ijH3ijpiMi 

jjjj;  uiijHjKIE-'li r  1!H!  ||IPh:'^IMIMlS3Uff 4|1j 


..  ..  'WiWk  h  ~ar||>1  >>j|i  ||IM'|  .  I  «'J  k-Pi 'llfl'  ilt -  :  Jll‘^1  ‘/-Hlflj  l  llllllil 

i  life .  *  fc  i=!i  iHBttiWfc !  Mi  HtnUHl* IKWiWlt  £ M  dblliiiih 
ipii  ini  -  srK'  «nnmratth.:,  -iiT;||fcimMIWwuiai  r<;==  nr  nii> 


ra«ll|lp|B4i . . 

. . .  ■  :,i  jfJL  It  VUIPJW'illiiSiP*., 

inrMI:  r-  i  B  k  Enr  ±c  =i  -lS±i  *  'mil  lirji::.  -  ii  : 


n\nHi 


.It,...*'*  41  „  i^./iitflifin  -?  nvui'lftlMlItHulHimiK..:'  rn.  "  —  ■  ■ 

V? Ia>. HU  :Hj  i:(f St  ill; ■  ^MNSSlil  I^IOl  t ■  J 1 » a^ia  'III Halt. >> 

huifl'T  :"if . nn  viMirri  m  t,j“ 

limuffi'.;*:.  *e  nr;»di:  iei  'KdiiiiiiiikiuntiiiiiiU'i-i  nwifpHftMfyiffniffliBiiB 

fegfcS imz Hi  ~  II'  II  l|i|  ■IIIHIMI'  II  ~  JfBTi~~ .  ' : 


|J|_3U. _a  =  ; i!4--Ht=  dHUHL  -IWUnlli  '■r'|'jr‘i|  IM 

lili:  Mr  ■H-leiEEsF7lHfl>='HL'fr!iii  '■ttM'Sh  **** 

!  "  -  -I!  _  !L  V Ltl  I  1_.  i  !■  ._"ia!iul_  I-,  I  ■  1  J  ■  .......  ~E  ■  .  -  i 


tfttUhnfcctLrtas^  >.<  .AL::-\  flllttlH!ii5^»R^Miiufl:ai 

:  iLfciHfisr  EiiiSEKHiiiFsiiiiillsiiii  -!?i  1\\B  i  *  HJF  Ifilliiirfl 

vm 

.  rtMMFWHlMliHh;1-"  . . Hit  UHMlEtl 

■  ■=j.NHi=i=  IN . .  -  Wfiaii— FT «■=  ■  feSfei 


.vHi;j|!L.|l!Bd  idltti  *■■■  IlMl;  IMUHHlIlHIlik^W ill 

===esHai'^_i:::ULJia=»a;  «ww4l=^:;2L .  nrric 


-s.H  itfMHKF  iH^T  :  J«'  '  '  ■  “Cl . 

til  -  s  ii = :  ^isyunsiiiEEiiS ■taiisiiiBEp ■  ifni^ii  i tta&miimiffl fli= 


SKiBircriliBirai 

■  ibSiat;  ill-  -I1  -l-_  -i 1  'r: 


iig^jgg 


ALTITUDE  - 150  KILOFEET 


INDEX 

SYMBOL  SHOCK  ANGLE 
EQ  NONEQ  (DEG] 


NUMBER  DENSITY  OF  ATOMIC  OXYGEN 


[gKiKagssssagsssBgasai 


::::a::si!:::;B::ESK:Rs::::::r.:siss!&sas!BSiK:::sB:: 


IMII  IMtl  ■■  ■■Mil 

siSHsitSSiisiss:::::: 


ssssi:s:si8ss:::ss=Bs:ss£::i£!i:::!ii3as:£!s!;!ii:!!sssi 


i»»sn:::BssssssBSBSKS! 


■■Mllll 


(i  mu  inn  inn  n  ■■  »ii  rnwrnmm 


■MuSuSSImMHiiiSUSSiiS! 


«»*■' »aJ 
leaat  nM 

■■■■  ■■■■■  ■■■■■■■in  i 
■■■■■■■■111111111111 
■■■■■■■■hhiiiiiii 
■■■■■■■Miiiiimi 


IIIIL.  Kill  (/*■■■■■■■  ■■■■■■■■■■■■■■■■■■«■  ■■■niiiiiMB  ■■  5 f ^  SI  - 

miih  . . . . . 

liiiimi  ■■■  ii ■■■■■■■  an iiiiiiiiii|HaiiiiliiHii ■■■■■■■■  ■■■■■■■■mi 

■  mil 


,  ,.VI^^^^n^mHlSpw8i9RSni8iSMHHiiSgjiSSSi|mSmn 

I SSSS!  SS  SSI  Zm  "  .'ok  SKm»CSS  S  8S8S  SSSSS  «■!•<  SSSSS  SSSSa  StSS!  SS  SS  SS  SS  "SSSSSSSSmS* 

I  miMUaiHHa  b  ■■■■■■■■  ■■■UMIII ■  ■■■■■■■■■■••■■■£ 

I  ■■■nil ■!!■■■■  I  IWIItM ■■■■■■  IMIimilUMIIIHI  mil M ■■■■■« ■■■ 

E■■■lM>|■|■■■■»  ^^■^'■■■■■aMaBRRaiaiain9Blf !!"!!!!!! SEES 5! £!SS59 
I  uni iiiiihm  ■■  'ik  wkiii  ■■■■  mil  urn  ■■■■■inn  mu  ■■■■ 11,1  aP""9 

I  llllllllll■■■■■■  % 

t  «^imiiiHnmmjiiiii^ 


lillllHHIIlUX  MlBHIMIIIIIIllJl 


ALTITUDE  =150  KILOFEET 
R  =  1.  27  CM 
W  INDEX 

SYMBOL  SHOCK  ANGLE 


■■  ■■ 


la  at  ■  •■•aaaaii ■ 

mb  ■■■■■■■  Mid 
!■■■■■■■■  miH 


_ !■■■■■■■■  mini  hi 


I  nil  nil. ■■■■■■■■  ■■■■■■■■■  U*  U  llli  ‘III  II  !■■>  '■  ■■■■■■■■■! Illlil**  ■■mill  Il|lip4 

iiirliii’SiEfiiiiliiuiHiitUiiii  ■■■MIL.  ■««  filial  fek?«iiuMia':.iiiii  III . . 

I  mil IIIIIMmSb ■■  II  ■■■■■■■■lllk1llll».'*pMIIIII  lllln^  II  !!■■■■■  HIM IIIII11II  IIM  illlil 

1  llfll  1111  ■■!■■■! Hill ^VII  IIII iMBII ■*  '  ■■■■■  IIUIH|B  1111! ■ 

la  aaaaa  BB  BB  MB  ■  !  MMB  BB  B  B  BB  P  BIBB!  IMBi  .  ‘BBS  ■  ■■■■!:■  .  ■«  Bl  BB  ■■  M  BBI B  a.  '  ^BBB  BUT  *.•■■■  BPP**  ■■■■’  ^■■BHI  ■MBIIIHMiimIIHII 

ImSiSSSi  BBaSSSiiSSaBBiBBiimmimBk "*■■■■■■  mu  a. ~  iubbibbibm.  -  ml  m*' .*■■■»•  .«*■  bbivbi i.  ■■■BiuBBiinuiHin 

;  m!SS!SI{iiiiiiiiliiiiHi>iiiiiiiHiimiiik.,'iiinMmi«  -nmn  ■■»■  .!■■  b*,,:mmii>iiii  ■■■■■nir 

I . mi . 


M _  mmb  . . .  m  fcilPPBi 

nEiiin.  5SSiSI«!SSmSij5«MMjjj!jj5iBiiiiiii 


■  ■■■■■•■ill 

I  ■■■■■mm 

■  ■■■■iiiiiii 
I  ■■■■■■uni 


I  . .  ■ 

iiaBBiiniBMaBBaaiiiiiiiiiiii  ■ 
iiBBaBaaaiBBBBBaaaaiiiaBaaUiia 


HiSS  SS  SSSSS  !mirSSS!njBB55rB5BBM»lMiiSSSSrMMii»iiiiS«SEBBri5Eii5iiiSIiiiIJ 

|s:ss::ss====sss::::s;ssi;::;ss:sss::si;:2s;sssrs:::s=ssssssss«s::::K;ss::sssss:Ri;sssss;ss=ssss;s:::::s:; 

[n  mi  minlH  im  ■■  ■■  ai  a  ■■■  tn  iin  im 
^kii  ■■■■■■■■■  ini  ■■  si  ai  c  ■■■mi  ■■;  up  mi  ini 

■hiinuimiiiiiiiMHiiiiHMaiMMiiiiiiiiii 

inimtik'.  mtt  liii’ttiiiiin  ■■  ■■  it  ■■  w  b#  iiiii  iiniHlii-  n  ft  !■'■■■■  ■■■it  imt  nil? 

liiiliiiiiiSSSlMMBiiiiiiiiiiiji . ■■iiiiiiiiiiiiiBaaniiiBBgBBiaiiiiiiiiiniiiiiiMiiiiiiiiBmniiBBBa»aaiiiiii 


lI!(.M  —  £wiaaa»aiiBE5SSraS??ir?r»«iitSSM«aSa^aT^ 

IMiaaB  mm  a  m  aanHaaMB  mmb  a  •■••a«*i«liaa««"»»Mi  ■ 

■Mb  MB  mm  mm  mb  hhmmb  Mb  mmb  I1BI  mMb  Ba  a  BiBi  m 

■  Bill  MMB  MMB  MB  MMMMB  BMMBB  ■■aillBIII  M MM ■  B  Bill  B 

^BBMMMlia  BMalBBaillB  . . . 

aiMBM  bB  BB  aagMBB  jMaai  a  a  BBa  a  a  Bat  BBB  BBBBBBJ 


ss  ss  si  ssssSbUss  ssss?  isns  I 

mm  ■■■■«■■■■■■■■■■  ■IIHHIII I 
MMUAMBMamHiiaaii  ■■■■•  I 
HMMini  I 


hS 


¥ 


nil:  Urn 


IBWtJPiftllWiSii' 

(RHUIkMijW 


U*  frV 


|  'IT  ill i  L:;!?/  ,Ky  isr  tt  ih"i!:  wm  -j;-i 

I  i|h  ju  triNi  vmm . . 

\m  '  4v:^eU!  ilg=sll!ilE;  ^BEimillliilEI39Si^lflSiyEi  slil:  .Sr!rf:! 

4il|  *L+MKHlL'.l!!'tr  ■  s  ^  -ftfjilHIill  «M*fMjhlMIN!j  UHHIlftlh.. 


mi 


Mfe^aMnwoniNiviu 


.r.w  ttll*  jMtWWu  HI.  "  % 

bimiWMi  «k^ 


■nisiHeii*  ’miiiswMel 


1 

i,jf.  . -jishiMHUa  :ni , 

III  HiUKfflHIffiHryPii, 


^liiiiHi  biinifegftf  p  iiiii^HiKLisiiiK 

t«M<  ^j^wwrna^SS!" 


iBsuiHin  iunjn!HMmABHniwwin»i,iit-.t4i*b:.:"' 

innM'nttHwiiiivi  *«Elii 

iK6>ifraimittiM®sw»Bf4  -'-mm 

,jh,»  i.^iniui. wim iuinriaa»HKiiuiiy^»M>~.  win  mu 

warn  in—B  m iwimvpi*^  &  iiti1! 

»s  mammmm  HMfeaiihl|>!l! 


®mjs3sas$m*8ia 


.:  : ; : : . : — l; :  -si**--  - 

r  wm-iSt* . . 

!1  iUmMHf-XH’iM 


I'nii*  Him  ■1*1(0  jii  ii  nmwr-BMiw  i  .Hu-iittiikiuirdttiwiH.Aiir"  sh-j?  '■■  ►=  r  ■  non 

L  m  (WM*  nea  ih'  t#>i  ■•- 1«  ujhwi  wjiKin*Miki»iiiiiMw»MMn«lit«tW 

•,  <•  **  :i:ii  asm . . 

.. ms  :::  ■  t  itfiijiij: :ii  iiw;aw»ipiBii..  ^iintitHiiiMtaftfiiimiii'aiittiwai 


rittiMflittuaiBsii 


iiiiiiiiinii! 


”  (Wll»l«l  :ipr  Will'  ■BiKfc  '■  -  «i  ■WpllL'iW!  ■  iiii Hi*  linflim.  ■  ...'nismn1 v%r i 

mmmsm  «HHflBiwnlU .:  WBiwiiiMtiflmiimmiii  mmmtwmeszszi 


>-  iiteiMn'it.siW  -■stf.  iii- 


-  :'i.4i  jKffliij:  irnm* 


HSlEMimii 

.U- ,  K>4i  ■IIUIUb.1  ..»■■■■! 

:™  1!ES5|51 


SYMBOL  SHOCK  ANGLE 
NGNEQ  PEG) 


ALTITUDE  *  150  KILGFEET 


LIS  CM 


Coordinate  Along  Body  Surface 


NORMALIZED  COORDINATE  UX3NG  BODY  SURFACE,  »^/R 


Figure  3-10.  Streamwise  Variation  of  Electron  Number  Density  vs.  Normalized 

Coordinate  Along  Body  Surface 
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Figure  3-11.  Normalized  Pressure  vs.  Shock  Layer  Thickness;  s./r  =  4. 962^; 

Nonequilibrium  N 


Figure  3-12.  Normalized  Pressure  vs.  Shock  Layer  Thickness;  s  /R  =  6.125^- 

Non  equilibrium  b  N  * 
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Figure  3-13.  Normalized  Pressure  vs.  Shock  Layer  Thickness;  s./R. 

Nonequilibrium 


Normalized  Temperature  vs.  Shock  Layer  Thickness 

Norequi  librium 
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SHOCK  LAYER  THICKNESS,  Y/Y^ 

Figure  3-16.  Normalized  Temperature  vs.  Shock  Layer  Thickness;  s^/R 

Nouequllibrium 
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Figure  3-18.  Normalized  Density  vs.  Shock  Layer  Thickness;  s^/R^  =  6.125^ 


3-71 


1 


Figure  3-20.  Normalized  Velocity  Parallel  to  Surfaoe  vs.  Shock  Layer  Thickness • 

sb/R=4.982V 

*- 
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Figure  3-2i.  Normalized  Velocity  parallel  to  Surface  vs.  Shock  Layer  Thickness;  e, /R^  =  6.125 


ELECTRON  NUMBER  DENSITY  N  (no./cc) 


Figure  3-23.  Electron  Number  Denelty  vs.  Shock  Layer  Thickness;  =  4. 962s\1 
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Normalized  Coordinate  Along  Body  Surface;  Nonequilibrium 


Figure  3-26 


Table  3-4.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles 


Species  Number  Density  (particles  per  cc) 

Altitude 

(kllofeet) 

Y/Y 

sn 

°2 

N2 

N 

0 

NO 

50 

0 

2. 723X1 

7.128. 

1.000N 

5.70X6 

6.520 

9.689 

6. 662  x8 

2.2.8^\9 

2.375 

2.  651 

2.509 

6.  785N'}8 
6.405 

5.228 

0 

1.296s}9 

1.329 

1.418 

0 

I.I2X8 

1.203 

1.476 

0 

100 

7. 742K2 
2.429V 
4.838 
7.005  0 
1.000  N 

8. 29(^J4 

8. 742 
9.556 
1.374s*5 
3.100 

6. 399s}7 

2. 056^  8 

2.  080 
2.1584 

2. 279 

2.384 

2.410 

,.mK" 

1.707 

1.660 

1.503 

1.195 

0 

1.515s}8 

1.526 

1.553 

1.  566 
1.510 

0 

3.810s}6 

3.  867 
4.161 
5.293 
8.962 

0 

150 

8. 06X2 

2.393V 

4.780 

7.044 

1.000X 

7.407s}2 

1.126X4 

2.655 
7.244  ,c 
1.910s}5 
6.717X6 

2. 097s}7 

2.  000 

2.127 

2.187 

2.238 

2.530 

2.  013\17 
1.838 

1.484 

7.l08k6 

0 

!.617\7 

1.562 

1.452 

1.313 

1.113 

0 

3. 75^s}5 
4.945 
8.684  - 
1.693s}6 
3.174 

0 

200 

7.98K? 

2. 706 X1 
4.782 

7.  098 

1.  000N 

s.sieX}3 

1.226V 

2*77ll5 

l.osr}5 

4.637 

1.027s}6 

3.07i^}6 

3.040 

3.068 

3.220 

3.463 

3.867 

2.194s}6 

1.707 

1.145  1K 

6. 54  ^}5 

2.973 

0 

2.  08(^6 
1.828 

1.541 

1. 324  „ 

8.  002s}5 

0 

1.624s}5 

2.  878 

4.567 

4.656 

2.588 

0 

250 

4. 89«X2 

1. 975\ 
4.612 

7.594 

1.000 

1.519s}4 

3.244 

6.774 

1.384 

3.731s}5 

3.845 

4.101 

4.412 

4.693 

5. 211 

8.  299s}4 

5.944 

3.517 

US°kl3 

6.  529s} J 

0 

1.522s}5 
1.252  ,  . 

7. 516s}4 
3.348 

1.062 

0 

5.1oX4 

4.  070 

2*  23\13 

8. 84<XJ 
2.459 

0 

Table  3-5.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles, 

V  rn  *  s-125\l _ 


Species  Number  Density  (particles  per  cc) 


Altitude 

(kllofeet) 

Y/Y 

sn 

°2 

N2 

N 

O 

NO 

50 

1.6lH6 

1.647^? 

2. 962s\6 

7.  956^ 

7.437s\6 

4.2lK® 

2.259 

1.699 

2.214 

8.009 

8. 426  - 

1.247\2 

5.303 

1.844 

9. 882H5 

7.917 

i.uK7 

2.416 

1.366 

2.018 

3*  10\.  14 

6.872 

1.523 

4.088 

2.873 

2.469 

7’62fl3 

6.238 

1.535 

5.726 

5.636 

3.169 

8*651  1A 

3.738  ,e 

2.473 

7.875 

1.094^8 

4.138 

1.614\18 

4.946s?8 

5.603^ 

i.  i45y 

1.847 

6.  955 

2. 339X1 

1.39$S 

1. 362\ 

1.  000\ 

3.413 

1.285^ 

0 

0 

0 

100 

1 

0 

1.10H2 

1. 573s'\7 

8.295\16 

1.  0G2\17 

7. 680V1 3 

3.  832 

1.206 

1.636 

8.057 

1.082 

7.306 

8.272 

1.351 

1.739 

7. 733 

1.128 

6.396 

1.279X 

1.706 

1.853 

6.865 

1.166 

6.684 

1.904 

2.582 

2.  048 

5.433 

1.231 

7.097 

2.439 

4.354 

2.  231 

3*844  1 . 

1.286 

9-189 

4.013 

5.815 

2.  558 

3.992V* 

1.184 

8. 13^!® 

5.907 

4.500X6 

2.  898 

4.402X8 

4.!84\46 

3-  00l\ 

8.275 

9.934 

3.763 

2.417^° 

6.729V" 

6.632^1 

1. 166\ 

1.615NJ7 

6.082  1  0 

5.249\  iK 

1.  733VSg 

1.936\19 

2.  214 

3.493 

1.315^8 

1.7518\15 

1.  005V 

1.49K9 

1.000 

3.330 

1.254 

0 

0 

0 

Table  3-5.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles, 

s  /R  =  6. 125V  (Cont) 


Species  Number  Density  (particles  per  cc) 


Altitude  Y/Y 
(ldlofeet) 


9.506^1 
3.69<K 
7. 235  _ 

1.176V® 
1.821 
2. 381 
3.932 
5.796 
7.596  , 
1.07  sV 
1.000  V 


4.829\?2 
3-251 
1.569  V3 
8.078  .. 
2. 2M\!4 
7.121 
1.126\15 
1.568 
2.695 
5. 256 


2.78K15 

2. 655 

2.559 

2. 543 

2.557 

2.674 

2.809 

3.288 

4.305 

5.908  .. 

1.015^6 

1.980 


1.780V5 

6. 763\ 
3.753 

7]  925^  3 
2.821  19 
1.83oV. 
2.864V° 
5.450 V 

2. 168V 
0 


1.95K15 
1.728 
1.523 
1.393 
1.239 
1.043 
8.013^4 
2.365  10 
2. 317V1 
4.  089\A 
3.47l\6 
0 


2.7!2\“ 

2.11l\, 

1.886\3 

6.783 

1.982^4 

3.214 

3.279 

la64\13 

1.435V3 

3.617V1 

3. 877V 

0 


8. 088\"t 

3. 543V 
9. 200  2 
1. 668 V 
2. 511 
3.639 
4. 668 
7.038 

9.495 
1.411X1 

2.495 

1.  000\ 


5.39^J2 

1.155V3 

2.575 

4.199 

5.486 

6.556 

7.713 

9.  14 

1.286V 

1.631 
2. 745 
6.982 
7.302 


S.eiT^3  1.223S\4  3.509X13 


2. 719 

2.711 

2. 769 

2.869 

3.010 

3.234 

3  971 

4.902 

6.M9 

1.034V 

2.622 

2. 750 


3.993 
2. 228 
1.209  19 
6.996  V 
4.101 
2.060 

9-455>io 

9.  069\ 
2.478\® 

3.  099 V 
6.678V15 


Tok13 

5.188 

3.460 

2.309 

1.394 

1.086  10 

2.382V12 

2.902V 

0 


3.411 

2.764 

1.959 

9.12lV 
5.583 
4.720 
1.163  ,, 
1.852V11 
1.805  V 
2.969V1 
0 


Table  3-6.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profiles, 

sb/%  - 2- 4  A2 


100 


150 


- 

Species  Number  Density  (particles  per  cc) 

Altitude 

(ldlofeet) 

Y/y 

sn 

°2 

N2 

N 

0 

NO 

50 

s.es^C4 

1.  077V1" 

2.092 

3. 575 
5.073 

7,011V2 
1.  013'. 

1.  000 V 

1.577V® 

2.205 

5.  234 
1.367V 

3. 003 

6.048  1Q 

1.100V8 

1.847 

3.467 

1.624V 

1.675 

1.818 

2.  030 

2, 327 

2.982 

4.162 

6.955  .. 

i.soeV 

2.  85  8V® 
2.132 
9.445\15 
3.154 

8.  382%. 

n  207^ 

2.34l\ 

0 

7. 853N!7 
7.910 

7.812 

6.929 

4.986 

2.577 

4.990^ 

1.40oV 

0 

7. 20JV® 

i.nU1 

1.532 

1.864 

1.583 

5. 729v%5 
1.363N 

0 

!.mV 

3.755 
8.123  _ 
1.253\3 
1.852 
2.349 
3.710 
5.265 
7.216  « 

l.OOlV 
1,871 
1.000V 


6.892^ 
2.177X4 
4.986 
8*  056  . 

1.272\3 
1. 713 
2. 973 
4. 366 
6.160 
8.868  Q 
1. 000  \ 


3\2 


7.503 
7.409 
7.253 
6.793 
6.877 

7,5<\l3 
1. 049S._ 

r.S\“ 

1.005V7 

1.625 

3.449 

3.357 

2.125^° 

1.960 

1.914 

2.544 

7* 945  . g 

1*030^13 

3.146V** 

1.866^® 

6*  355  1  £ 

1.063\16 

1.749 

3.405 
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1.927^® 
1.929 
1.934 
1.960 
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2.  814 
4.009 
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0 
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2.  271 \ 

8.  016^® 
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8.570^® 
8.621 
8.801 
9. 160 
9.618  . 
1.044V 
1.126 

4a99\J® 

7.574\® 

0 

1.48K® 
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1.49lNg3 

1.969V 


4.821sv4 
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4.600 

4.224 

3. 984 

3.626 

3.696 

8.444 

2.334^® 

7.368%. 

1.968^? 

1.109V8 

0 

4.386S\1 
5.525  ^2 
1.  031V 
3.102 
9.555 

3. 290%* 
2.487V 
1.166  to 
1.399V3 

2.189  V 
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Table  3-6.  Altitude  Dependence  of  the  Species  Number  Density  Normal  Profile 

Sh/RN  '  2-  433\2  <C0»t) 


Species  Number  Density  (particles  per  cc) 


(kilofeet) 
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CD 

• 

f-H 

2.  868 

2.  094 

8.794\4 

3. 758  1 

3. 255 

6.715  ,c 

3.355 
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0 

0 

0 
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234 
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4.892 
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2.202 

1.476 
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0 
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It  is  worth  noting  in  Figures  3-11,  and  3-13  that  the  pressure  profiles  obtained  from  the 
perfect  gas  or  equilibrium  flow  fields  actually  follow  the  free  stream  Mach  number  variation 
rather  than  the  altitude  variation.  Also,  for  these  pressure  profiles,  which  are  derived  from 
the  perfect  gas  or  equilibrium  streamwise  pressure  distributions  and  reproduced  during  the 
construction  of  the  approximate  nonequibbrium  flow  field,  the  pressure  behind  the  shock  is 
replaced  by  the  value  from  the  thermal  equilibrium  frozen  shock  flow  solution.  This  change 
at  Y/Yg^  =  1  is  particularly  noticeable  on  the  first  normal,  where  the  result  is  a  lower  value 
for  the  equilibrium  cases  and  a  higher  value  for  the  perfect  gas  case. 

Some  attention  should  also  be  directed  to  the  development  of  an  off- wall  peak  in  the  temperature 
profiles  back  on  the  conical  surface  (Figures  3-15  and  3-16)  as  the  altitude  is  increased. 

This  trend  is  unquestionably  one  of  the  more  marked  effects  of  nonequilibrium  flow,  inasmuch 
as  the  temperature  for  equilibrium  flow  decreases  monotonically  from  the  wall  to  the  shock 
(except  in  the  region  where  the  streamlines  have  traversed  that  portion  of  the  shock  where 
the  minimum  shock  angle  occurs  as  a  result  of  the  overexpansion). 


3.  2. 2  VIBRATIONAL  NONEQUILIBRIUM  CALCULATIONS 

The  relative  importance  of  vibrational  nonequilibrium  on  the  shock  layer  flow  properties  has 
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been  investigated  by  applying  the  following  superpositional  method.  Landau  and  Teller 
have  shown  that  the  Master  Equation  for  a  system  of  harmonic  oscillators  reduces  to  the 
familiar  vibrational  relaxation  equation: 

dey  1 

*  -  y  <£v~-  €v>  <3-3> 


where  y  denotes  the  relaxation  time  of  the  system;  the  equilibrium  vibrational  energy; 
and  the  vibrational  energy  at  time  t.  For  the  harmonic  oscillator 


Q/T 

S  =  RV  (e  V  -1) 
e  /r 

e  =  R0  /  (e  v  -1) 
v»  v  ' 


(3-4) 

(3-5) 
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in  which  R  is  the  universal  gas  constant;  6^  is  the  characteristic  vibrational  temperature;  T, 
the  kinetic  or  gas  temperature;  and  T^,  the  vibrational  temperature.  For  applications  in 
which  the  relaxation  time  and  the  equilibrium  vibrational  energy  are  implicit  functions  of 
time,  the  general  solution  for  equation  1  is; 


t  e  t*  i  t 

‘  _£!_  f  dt**  1  dt* 

y  J  ~Z  dt*  +  const  exp-  f  — 

exp  o  T  J  J  7 


(3-6) 


r  t  e  -  e 

=  /  — , 


exp  f 


c  dt* 

exp  -  I  -  +  € 

y  vo 


(3-7) 


in  which  t*  and  t**  are  dummy  variables  for  time.  The  second  equation  is  obtained  from 

applying  the  boundary  conditions  at  t  =  0,  e  =  e  . 

v  vo 

For  steady  flow,  the  distance  along  a  streamline,  s,  and  the  ordered  velocity,  V,  are  related 
by 


so  that  the  integrals  in  Equation  3-7  may  be  transformed  into; 


(3-8) 


r  s  €  -  € 

f  V00  VO 

ev =  j  “vv  “P 


j  ds*j 


/ds* 

W  +  €vo  <3‘9> 

o 


Using  Equations  3-4  and  3-5  the  vibrational  energy  may  be  rewritten  as; 


[  /  vy  (exp(0v/T)-l  '  exp  <9/T)-l) 


s* 

exp  / 


/ds*  1 

-  +  - — - 

Vy  exp  (S  /T  )-l 

o  V  vo 


(3-10) 
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and  the  local  vibrational  temperature  along  the  streamline 


is  given  by 


T  =9  /in  (1  +  R  0  /  e  ) 

v  v  '  V  V  (3-11) 

By  superimposing  the  solution  of  Equations  3-10  and  3-11  over  the  thermal  equilibrium 
calculations  for  the  invisdd  flow  field  at  250  Idlofeet,  the  relative  importance  of  the  vibrational 
nonequilibrium  problem  can  be  assessed.  The  highest  altitude  flow  field  was  selected  for 
investigation  because,  since  vibrational  relaxation  times  vary  inversely  with  pressure, 
conditions  for  vibrational  nonequilibrium  are  most  favorable  at  higher  altitudes.  Based  on 
this  simplified,  order-of-magdtude  analysis  it  is  concluded  ...at  vibrational  nonequilibrium 
can  cause  substantial  changes  m  the  flow  properties,  particularly  the  gas  temperature  and 
density,  and  consequently  must  be  taken  into  account  to  accurately  pr*  Hct  the  shock  layer 

properties  at  higher  altitudes.  The  analysis  and  the  results  which  form  tie  basis  for  this 
conclusion  arc  discussed  below. 

The  superpositional  method  used  in  this  investigation  consists  of  solving  the  uncoiled 
streamwise  vibrational  relaxation  equation  using  superposition  with  the  thermal  equilibrium 
chemical  nonequilibrium  invisdd  flow  field  to  provide  the  values  of  the  gas  pressure  and 
temperature  required  to  calculate  the  local  relaxation  time  and  vibrational  equilibrium 
energy.  Only  the  primary  molecular  constituents  N2  and  Og  were  considered  and  the 
relaxation  times  for  each  of  these  species  were  obtained  from  Millikan  and  White  208 ; 

e  (129T"1/3- 22.3) 

y  -  i 

°2  P  (3-12) 

(220T_1/3  -  24.  8) 

N2  P  (3-13) 

where  the  relaxation  times  >fc2  andy^  are  in  seconds,  p  is  the  gas  pressure  in  atmospheres 
and  T  is  the  gas  temperature  in  degrees  Kelvin.  The  harmonic  oscillator  model  is  used  to 
define  the  vibrational  temperature  corresponding  to  the  vibrational  energy  obtained  from  the 
solution  of  the  relaxation  equation  and,  therefore,  is  given  by 
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VV  ln«*rovAV  <s-i4> 

In  the  above  expression  T  and  e  are  the  vibrational  temperature  and  energy,  R  is  the 
universal  gas  constant  and  6  is  the  characteristic  vibrational  temperature  which  has  a  value 
of  2239  K  for  oxygen  and  3395°K  for  nitrogen. 

Some  typical  results  of  the  analysis  are  shown  in  Figures  3-27,  3-28, and  3-29  for  streamlines 

which  intersect  the  bow  shock  at  shock  angles  of  89,  70, and  46  degrees,  respectively.  In 

these  figures  the  streamwise  variation  of  the  vibrational  temperatures  of  O  and  N  as  well 

2  2 

as  the  gas  temperatures  obtained  from  solutions  for  thermal  equilibrium  chemical  nor 
equilibrium  and  for  perfect  gas  with  y  -  7/5  (undissociated  air  with  no  vibrational  excitation) 
are  given  as  functions  of  the  normalized  di’.ance  along  the  body  surface.  In  all  cases  the 
vibrational  temperatures  of  both  oxygen  and  nitrogen  are  characterized  by  an  initial  raoid 
rise  from  the  free  stream  value  at  the  shock.  However,  as  the  streamline  continues  to 
expand  around  the  body,  the  drop  in  the  gas  pressure  and  temperature  causes  a  rapid  increase 
in  the  local  relaxation  time  with  a  corresponding  rapid  leveling  off  ("freezing")  of  the 
vibrational  temperature.  If  at  this  point  the  vibrational  temper  ature  is  less  than  the 
equilibrium  temperature  it  will  continue  to  rise  gradually  until  the  two  are  equal, which  for 
this  application,  particularly  for  nitrogen,  requires  a  considerable  distance  along  the  body. 

If  the  vibration  temperature  at  the  leveling  off  point  is  greater  than  the  equilibrium  temperature 
it  will  decrease  gradually  until  the  two  are  equal.  It  can  be  seen  that  this  latter  situation  is 
the  case  for  the  89  degree  shock  angle  streamline,  whereas  the  former  situation  applies  for 
the  other  two  streamlines.  Also  for  the  89  degree  streamline  ft  may  be  noticed  that  the 
difference  between  the  no  vibrational  excitation  temperature  and  the  thermal  equilibrium 
temperature  back  on  the  conical  portion  of  the  body  is  about  2500°K  or  about  45  percent  of  the 
equilibrium  temperature.  In  this  case,  however,  this  temperature  difference  should  not  be 
interpreted  as  indicative  of  the  extent  of  vibrational  nonequilibrium  effects  because  both  the 
oxygen  and  nitrogen  vibrational  temperatures  are  above  the  thermal  equilibrium  temperature. 

As  a  result  of  this,  if,  in  the  solution  of  the  coupled  problem,  gas  temperature  changes  from 
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Figure  3-27.  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temoeraturc  vs. 

Normalized  Coordinate  Along  Body  Surface 
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Figure  3-28.  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temperature  vs 

Normalized  Coordinate  Along  Body  Surface 
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Figure  3-29.  Vibrational  Nonequilibrium;  Streamwise  Variation  of  Temperature  vs. 

Normalized  Coordinate  Along  Body  Surface 
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thermal  equilibrium  do  not  significantly  affect  the  chemical  activity,  the  gas  temperature 
from  the  coupled  solution  would  be  below  the  equilibrium  temperature  shown  in  Figure  3-27. 

To  evaluate  the  extent  of  vibrational  nonequilibrium  effects  in  this  case  would  require  a 
comparison  with  the  gas  temperature  from  a  perfect  gas  solution  withy  =  9/7  (full  vibrational 
excitation).  In  the  other  two  cases,  however,  where  the  vibrational  temperatures  are  below 
the  thermal  equilibrium  temperature,  the  gas  temperature  from  the  coupled  solution  would  be 
between  they  =  7/5  curve  and  thermal  equilibrium  curve  and,  consequently,  the  difference 
between  these  two  curves  would  indicate  the  extent  of  the  vibrational  nonequilibrium  effects. 

3.2.3  STAGNATION  POINT  MATCHED  INVISCID-BOUNDARY  LAYER  CALCULATIONS 
Profiles  of  the  flow  properties  along  the  stagnation  streamline  within  the  shock  layer  have 
been  obtained  from  matched  invisdd  -  boundary  layer  chemical  nonequilibrium  pure  air 
solutions  for  altitudes  of  200  and  150  kilofeet.  These  profiles  are  required  as  initial  conditions 
for  the  nonequilibrium  boundary  layers  which  are  presented  in  Paragraph  3.  2.  4;  they  are, 

however,  of  considerable  interest  in  themselves  and  as  such  are  presented  separately  from 
the  boundary  layer  solutions. 


In  shock  layer  applications,  the  stagnation  line  matched  inviscid-boundaiy  layer  problem  is 
coupled,  i.  e. ,  the  inviscid  solution  must  be  used  to  provide  edge  conditions  for  the  boundary 
layer  solution  which  in  turn  must  be  known  to  locate  the  bow  shock  as  influenced  by  boundary 
layer  displacement.  The  matched  solution  must,  therefore,  be  obtained  by  iteration.  The 
procedure  used  herein  is  to  first  assume  the  boundary  layer  profiles  necessary  to  locate  the 
edge  at  which  the  inviscid  solution  is  tenninated  thereby  supplying  the  necessary  edge  conditions, 
The  boundary  layer  solution  with  these  edge  conditions  is  obtained,  these  profiles  are  then 
used  to  relocate  the  edge,  and  the  process  is  repeated  until  the  change  in  the  edge  is  less  than 
a  specified  small  tolerance  after  which  the  solutions  are  considered  matched. 


3-92 


Before  proceeding  to  a  description  of  both  the  invisdd  and  boundary  layer  solutions  it  is 
appropriate  to  point  out  that  both  solutions  have  identical  thermodynamic  and  nonequilibrium 
chemical  reaction  systems  and  both  assume  the  condition  of  thermal  equilibrium. 

3. 2.  3.1  The  Invlsc  d  Layer  Solution 

In  the  invlscid  region  of  the  shock  layer,  the  solution  is  provided  by  the  streamtube  method 
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of  McMenamin  and  O'Brien  which  was  described  in  detail  in  Paragraph  3. 2. 1. 2.  The 

governing  variable  used  with  the  streamtube  method  is  the  streamwise  velocity  which  is 

2q9 

obtained  from  a  correlation  of  data  from  Burke,  Curtis,  and  Boyer  which  are  the  result 

of  the  inverse  shock  nonequilibrium  invlscid  shock  layer  solution  of  Cornell  Aeronautical 
97 

Laboratory  .  The  correlation  which  reproduces  the  data  of  reference  to  within  +  8  percent 
is  given  by: 

v  =  vs  (!  -  s/  A)b  (3-15) 

in  which  v  is  the  streamwise  velocity,  v  is  the  velocity  behind  the  shock  wave,  A  is  the  shock 

s 

detachment  distance  and  the  exponent  is  given  by: 


b  =  1 

for 

b  =  0. 000596  v 

for 

s 


v  £  1670  ft/sec 
s 

vg  2: 1670  ft/sec 


The  nonequilibrium  shock  detachment  distance  required  for  the  invlscid  shock  ip.yer 

solution  is  obtained  in  the  following  maimer.  By  examining  a  multitude  of  shock  detachment 

distances  from  equilibrium  air  and  perfect  gas  (1  £  y  £  5/3)  flo\  ?  field  solutions  for  spherical 
210 

noses,  Storer  found  that  all  cases  could  be  correlated  with  excellent  agreement  by  the 
formula 


(3-16) 
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In  which  is  the  sphere  radius,  p^  is  the  free  stream  density,  and  pg  is  the  density  behind 

a  normal  shock  wave  corresponding  to  either  the  equilibrium  air  or  perfect  gas  shock 

solution  lb  the  case  may  be.  insofar  as  the  gas  density  in  the  stagnation  region  for 

either  equilibrium  air  or  a  perfect  gas  is  nearly  constant,  it  might  be  expected  that  the  shock 

detachment  distance  could  be  essentially  specified  by  the  density  behind  the  normal  shock. 

This  cannot  be  true  for  the  nonequillbrium  flow,  however,  because  the  density  varies  from 

the  frozen  value  ( y  -  7/5)  at  the  shock  wave  to  the  equilibrium  value  at  the  stagnation  point. 

Nevertheless,  because  Storer’s  formula  applies  equally  well  to  both  of  the  extremes  of 

nonequilibrium  flow,  it  might  be  expected  to  apply  also  in  the  nonequilibrium  case  if  the 

density  behind  the  shock  were  replaced  by  a  suitable  average  density  along  the  stagnation 

streamline.  Indeed,  this  has  been  found  to  be  the  case  and  good  agreement  with  the  non- 
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equilibrium  shock  detachment  data  of  reference  has  been  obtained  whenog  in  Storer's 
formula  is  replaced  by 


^FROZ 


A 

L 


FROZ 
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p  ds 


(3-17) 


in  which  A  is  the  frozen  shock  detachment  distance,  and  p  Is  the  density  from  the  non- 
FROZ 

equilibrium  stream  tube  solution  which  is  obtained  using  the  velocity  correlation  with  A  =  A^^. 
3. 2.  3. 2  Boundary  Layer  Solution 

In  the  boundary  layer  region  of  the  shock  layer  the  solution  is  provided  by  the  implicit  finite 
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difference  method  of  Blottner  .  At  the  stagnation  point  the  boundary  layer  equations 
degenerate  to  a  set  of  ordinary  nonlinear  differential  equations  with  two-point  boundary 
conditions.  After  a  transformation  to  the  Howarth-Dorodnitsyn  coordinates  system, 

Blottner  solves  the  linearized  set  of  governing  equations  by  an  implicit  finite  difference 
scheme  which  requires  the  successive  inversions  of  a  tridiagonal  matrix.  The  governing 
equations  which  must  be  solved  are  of  the  form 
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in  which  t|  is  the  transverse  Howarth-Dorodnitsyn  coordinate  and  H  represents  the  longitudinal 

velocity  ratio  u/ug,  the  temperature  ratio  T/T^,  the  species  mass  fractions  or  the 

element  mass  fractions  (pending  on  whether  the  momentum,  energy,  species  continuity, 

or  element  conservation  equation  is  being  considered.  Because  of  the  linearization  process 

the  coefficients  a  ,  a  ,  a  and  a  depend  on  the  unknown  variables  H.  Thus,  to  start  the 
1  2  3  4 

solution,  initial  profiles  of  the  dependent  variables  must  be  provided  from  which  to  calculate 
the  coefficients.  The  set  of  governing  finite  difference  equations  forms  a  tridiagonal  matrix 
which  yields  the  solution  of  the  dependent  variables  after  inversion.  Based  on  these  values 
of  the  dependent  variables; the  coefficients  are  recalculated  and  the  process  repeated  until 
the  change  in  the  dependent  variables  for  successive  iterations  is  less  than  a  specified  small 
tolerance. 

Calculations  for  a  0. 5  inch  nose  radius  body  at  velocities  of  22. 0  and  21. 95  kilofeet/sec  at 
altitudes  of  200  and  150  kilofeet  respectively  have  been  done.  The  results  of  the  matched 
inviscid-boundary  layer  stagnation  line  solutions  are  presented  in  Figures  3-30  through  3-36 
where  successfully  in  groups  of  two,  the  invisdd  and  boundary  layer  solutions  are  compared 
at  altitudes  of  200  and  150  kilofeet  for  the  gas  temperature,  density,  electron  density,  and 
pressure.  In  both  cases  the  boundary  layer  edge  was  arbitrarily  defined  as  the  point  where 
the  derivative  of  the  longitudinal  velocity  ratios  equaled  1/10  i.  e. ,  at  the  point  where 

du/ 

— =  1/10  (3-19) 

dtl 

In  both  cases  the  boundary  layer  displacement  thickness  is  negative  due  to  the  cold  wall 
condition  and,  it  is  so  small  as  to  be  unnoticeable  in  the  results.  Also,  in  both  cases  the 
catalytic  wall  condition  (thermodynamic  equilibrium)  is  imposed. 

It  should  be  noted  that  the  rapid  decrease  in  temperature  and  corresponding  increase  in  density 
for  the  invisdd  solution  near  the  wall  is  caused  by  the  equilibrium  condition  which  is 
intrinsically  required  by  reacting  invisdd  flows  at  a  stagnation  point. 
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Figure  3-31.  Matched  Invisdd  -  Boundary  Layer  Solutions;  Gas  Temperature  vs.  Normalized 

Distance  from  Shock  Wave 
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Figure  3-32. 
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Figure  3-33.  Matched  Inviscid  -  Boundary  Layer  Solutions;  Gas  Density  vs.  Nonnaliaed 

Distance  from  Shock  Wave 
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Figure  3-34. 


Matched  Invlscid  -  Boundary  Layer  Solutions;  Electron  Density  vs.  Normalized 
Distance  from  Shock  Wave 
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Figure  3-35.  Matched  Inviscid  -  Boundary  Layer  Solutions;  Electron  Density  vs.  Normalized 

Distance  from  Shock  Wave 
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Figure  3-36.  Matched  Inviscid  -  Boundary  Layer  Solutions;  Gas  Pressure  vs.  Normalized 

Distance  from  Shock  Wave 


3.2.4  INVISCID  LAYER-BOUNDARY  LAYER  CALCULATIONS 

Profiles  of  the  nonequilibrium  flow  properties  (temperature,  velocity,  density,  and  electron 
density)  at  body  stations  corresponding  to  s^/R^^  =  0. 4962,  61. 25,  and  243. 3  are  shown  in 
Figures  3-37  to  3-48.  In  each  of  these  figures  the  dashed  lines  represent  the  results  of  a 
pointed  cone  teflon-air  boundary  layer  calculation  with  simulated  nose  bluntness  and  a 
noncatalytic  wall.  Nose  bluntness  effects  on  electron  production  were  simulated  in  the  pointed 
cone  calculation  by  adding  electrons  at  the  boundary  layer  edge  in  the  tip  region  so  that  the 
integrated  electron  mass  flow  at  a  station  farther  back  on  the  body  matches  that  of  a 
corresponding  blunt  cone.  This  procedure  was  necessary  for  the  teflon  calculations  because 
the  teflon  chemistry  has  not  yet  been  incorporated  into  the  blunt  body  boundary  layer.  The 
integrated  electron  mass  flows  were  obtained  from  the  blunt  body  pure  air  calculation  at  150 
kilofeet  and  scaled  for  the  two  higher  altitudes.  The  matching  station  was  taken  at  2-1/2  nose 
radii  back  on  the  conical  section.  From  the  matching  station  to  the  end  of  the  bodyt  teflon 
was  injected  at  a  rate  inversely  proportional  to  the  square  root  of  the  distance  along  the 
surface,  and  the  wall  temperature  was  equal  to  1800°R.  Also,  at  the  present  time,  the 
noncatalytic  wall  option  is  the  only  one  available  with  mass  adoition  other  than  air. 

The  broken  dashed  lines  in  each  of  the  figures  represent  the  results  of  pure  air  blunt  body 

boundary  layer  calculation  with  no  mass  transfer  and  a  catalytic  wall  at  a  temperature  of 
o  2 

1800  R.  Based  on  Frobstein's  criterion,  the  solution  at  altitudes  of  150  and  200  kilofeet  are 
in  the  boundary  layer  -  inviscid  layer  and  the  viscous  layer  regimes,  respectively,  and  for 
these  two  altitudes  the  inviscid  layer  solutions  are  shown  as  solid  lines.  However,  at  an 
altitude  250  kilofeet,  the  flow  is  in  the  merged  layer  regime  and  in  this  instance  the  results 
were  obtained  using  the  boundary  layer  solutions  with  the  shock  layer  increased  to  induce  the 
shock  transition  zone  and  the  edge  conditions  obtained  from  the  adiabatic  Rankine-Hugoniot 
relations.  The  results  at  this  altitude  are,  therefore,  only  approximate  because  the  effects 
of  transverse  curvature,  normal  pressure  gradient,  and  temperature  and  velocity  jump  at  the 
edge  have  not  been  properly  accounted  for. 
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At  the  two  lower  altitudes  the  first  iteration  in  matching  mass  flows  of  the  inviscid  and 
boundary  layer  solutions,  to  determine  the  proper  boundary  layer  edge  conditions,  has  been 
done.  It  is  significant  to  note  that,  although  the  two  solutions  have  not  been  matched  exactly, 
a  smooth  joining  could  be  accomplished  by  additional  iterations  without  apparent  difficulty. 
This  plainly  indicates  that  in  vorticity  interaction  problems,  the  only  real  difficulty  lie?  in 
determining  at  what  finite  distance  the  boundary  layer  edge  is  located.  For  the  two  cases 
presented  here,  tho  boundary  layer  edge  is  arbitrarily  prescribed  at  the  same  value  of  the 
transformed  normal  coordinate  determined  as  the  edge  in  the  stagnation  point  solutions  of 
Paragraph  3.2. 3.  It  is  also  interesting  to  note  that  back  on  the  cone  surface  the  entropy  layer 
is  contained  within  the  boundary  layer. 
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Figure  3-40,  Electron  Density  vs.  Shock  Layer  Thickness 
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Figure  3-41,  Normalized  Temperature  vs.  Shock  Layer  Thickness 
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Figure  3-42.  Normalized  Velocity  Parallel  to  Surface  vs.  Shock  Layer  Thickness 
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Figure  3-43,  Normalized  Density  vs.  Shock  Layer  Thickness 


Figure  3-45.  Normalized  Temperature  vs.  Shock  Layer  Thickness 
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Figure  3-47.  Normalized  Density  vs.  Shock  Layer  Thickness 
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SECTION  4 
ANTENNAS 


The  importance  of  the  antenna  in  a  reentry  transmission  system  goes  beyond  the  usual 
factors  associated  with  electrical  performance,  radiation  pattern,  cuid  efficiency.  The 
antenna  near-field  distribution  has  an  intimate  effect  on  the  form  and  magnitude  of  the 
influence  of  the  reentry  induced  environment  on  the  antenna  performance.  That  this  is  true 
is  fairly  obvious  for  both  linear  and  nonlinear  plasma  effects,  since  the  reentry  plasma  is 
formed  in  the  immediate  vicinity  of  the  antenna.  But  it  is  especially  true  in  the  case  of 
antenna  breakdown  effects. 

The  near-field  distribution  affects  the  breakdown  characteristics  of  an  antenna  in  various 
ways.  The  effect  on  the  characteristic  diffusion  length  is  the  most  important,  even  in  the 
static  environment  inside  a  vacuum  chamber.  In  the  reentry  environment  the  near-field 
distribution  is  partially  responsible  for  determining  which  flow  field  region  will  breakdown, 

(i.  e.  the  region  having  the  highest  field  intensity).  Also  the  near-field  distribution  affects 
the  antenna  impedance  and  thus  it  controls  the  relationship  between  the  input  power  at 
breakdown,  which  is  significant  to  the  transmitting  system,  and  the  breakdown  field  strength, 
which  is  predicted  by  the  theory. 

This  section  is  devoted  to  a  survey  of  the  problem  of  antenna  near-field  calculation  and  the 
derivation  of  simple  estimates  of  the  ratio  of  the  characteristic  diffusion  length  to  the 
wavelength  of  the  RF  signal. 

4.1  SURVEY  OF  HELD  DISTRIBUTIONS 

The  distribution  of  radiated  fields  of  an  antenna  can  in  principle  be  determined  from  a 
solution  of  Maxwell's  equations  with  appropriate  boundary  conditions.  If  the  distribution  of 
currents  and  charges  on  a  given  antenna  is  known ,  then  the  field  determination  reduces  to 
the  evaluation  of  definite  integrals  involving  these  currents  and  charges.  The  difficult  part 
of  the  problem  lies  in  determining  the  distributions  of  the  source  currents  and  charges. 
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These  source  distributions  are  often  assumed  in  order  to  make  the  problem  tractable,  and 
in  the  calculation  of  far-field  patterns  and  total  radiated  power  this  assumption  is  not  very 
critical.  But  in  the  determination  of  the  near  field  and  the  antenna  impedance,  the  assump¬ 
tion  of  a  source  distribution  may  lead  to  serious  errors. 

The  current  distribution  on  an  antenna  which  consists  of  good  conductors  is  often  solved  for 
by  integral  equation  techniques.  Such  integral  equations  are  based  on  the  requirement  that 
the  tangential  electric  field  at  the  surface  of  the  antenna  must  vanish,  whei°  the  field  at  a 
given  point  is  determined  by  an  integral  over  the  current  distribution  at  all  other  points. 

011  212 

King"  and  Hallen  have  made  extensive  calculations  of  currents  and  fields  on  finite 
cylindrical  antennas,  and  Harrington  has  studied  thin  wire  antennas  in  detail.  While 
many  variations  on  the  techniques  of  solution  of  these  equations,  such  as  variational  and 
perturbation  methods,  have  been  studied;  it  is  possible  to  place  them  all  on  a  common 
ground  in  terms  of  numerical  methods  of  solution  of  integral  equations.  Unfortunately  the 
numerical  analysis  nature  of  the  theory  makes  it  d'Tlcult  to  generalize,  except  wnen  the 
form  of  the  antenna  i »  restricted  to  a  very  simple  type. 

The  electric  field  E ,  caused  by  a  source  current  density  distribution  J,  is  given  by  (MKS 
units) 


i.M  w  /  ikl  r~-~r '  I 

E(r5  J  T<~>  C,  .  t  ,  <•»’  <4-D 

I r  -  r'  I 

where  Is  the  permeability  of  free  space,  u>  is  the  radian  frequency,  i  is  the  wave  number, 
r1  is  the  radius  vector  in  the  source  coordinates,-?' is  in  the  field  coordinates,  and  dv'  is 
the  source  volume  element.  This  equation  is  useful  in  deriving  an  integral  equation  for  the 
unknown  current  distribution,  as  well  as  in  calculating  the  field  distribution  for  an  assumed 
scarce  current  distribution. 
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4.2  ILLUSTRATIVE  CALCULATIONS 


The  main  concern  in  this  section  is  the  derivation  of  simple  estimates  of  the  effective 
diffusion  length  for  the  prediction  of  antenna  breakdown.  As  discussed  in  Section  6.i,  the 
diffusion  length  in  classical  breakdown  theory  comes  from  the  solution  of  the  equation. 

7 2  *  -*-(»VD)  *  =  0  (4-2) 

where  v  is  the  high  frequency  ionization  coefficient,  and  D  is  the  diffusion  coefficient.  In 

the  uniform  field  theory  of  cavity  breakdown,  the  boundary  conditions  on  the  unknown 

function  $  determine  the  solution  to  this  equation.  If  the  function  'I'  is  made  to  vanish  at 

2  2 

the  walls, the  operator  V  has  an  eigenvalue  (-1/A  ),  where  A  is  the  diffusion  length. 

2 

This  result  leads  to  the  breakdown  condition  v  =  D/A  . 

In  antenna  breakdown  theory ,  it  is  artificial  to  appeal  to  boundary  conditions  on  the 
vanishing  of  'i  in  Equation  (4-2),  since  the  fields  are  not  confined  within  a  vessel  with 
electron-absorbing  walls.  In  this  case  the  coefficient  v/D  is  not  constant,  since  it  depends 
on  the  electric  field  intensity.  No  eigenvalue  exists,  except  by  analogy  with  the  constant 
coefficient  case  in  terms  of  the  value  of  the  Lap*  ician  at  a  given  point.  The  proper  boundary 
conditions  are  that  'it  he  continuous  and  have  a  continuous  gradient  and  that  it  vanish  at 
great  distances  and  at  the  surface  of  the  antenna. 

A  useful  connection  between  the  uniform  field  theory  and  antenna  breakdown  theory  may  be 

214  215  216 

found  in  the  results  of  Herlin  and  Brown  *  and  MacDonald  .  These  authors  studied 
the  effects  of  nonuniform  fields  on  the  diffusion  lengths  of  cavities ,  using  a  power  law 
dependence  of  v/D  on  E.  The  results  show  that  an  effective  diffusion  length  can  be  assigned 
at  the  peak  field  point,  depending  on  the  field  Intensity  distribution.  The  effective  diffusion 
length  for  cylindrical  cavities  depends  on  the  geometry,  especia’*  *  at  small  cavity  radius. 
The  effective  diffusion  length  in  spherical  cavities  can  be  a  fact.  i  the  order  of  one-half 
the  geometrical  diffusion  length.  Assuming  that  the  field  distribution  is  the  controlling 
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factor,  in  light  of  the  fact  that  the  ionization  frequency  is  a  very  strong  function  of  field 
intensity,  then  the  width  of  the  high  field  ~^kion  of  an  antenna  hould  be  related  to  its 
diffusion  length.  By  analogy  with  the  geometrical  diffusion  length,  the  antenna  diffusion 
length  should  be  on  the  order  of  1/"  times  t  je  width  of  the  high  field  region. 


A  simple  illustrative  case  for  microwave  antennas  is  the  waveguide  opening  into  a  ground 
plane.  Assuming  that  the  waveguide  is  operated  in  the  dominant  mode  near  cutoff,  the 
wavelength  A  is  related  to  the  guide  dimensions  a  and  b  by 


1/A2  =  (l/2a)2  +  (l/2b)2 


(4-3) 


But  the  field  distribution  in  the  vicinity  of  the  aperture  will  be  dependent  on  a  and  b  such 
that  the  following  is  approximately  true  for  the  geometrical  diffusion  length. 


1/A2  -  (n/a)2  +  (n/b) 


(4-4) 


Thus ,  for  a  slot  fed  by  an  open  ended  waveguide 


A  s:  2  n  A 


(4-5) 


Whether  the  above  approximate  value  of  the  ratio  A/A  is  valid  for  antennas  operated  over 
the  entire  radar  spectrum  is  open  to  some  question.  At  the  long  wavelength  end  of  the 
spectrum,  there  is  considerable  payoff  in  vehicle  design  to  miniaturize  antennas  to  be 
mounted  on  reentry  vehicles.  The  antenna  dimensions  may  be  forced  to  be  small  compared 
to  the  wavelength,  by  utilizing  special  techniques  such  as  end-fire  design.  As  an  example 
of  such  an  antenna  on  which  much  data  have  been  obtained,  consider  the  conical  equiangular 
spiral  or  log  spiral  antenna.  This  antenna  has  been  analyzed  theoretically  by  integral 
equation  techniques217'218  and  the  radiation  patterns  and  near  fields  have  been  measured 
by  Dyson219’220.  This  antenna  consists  of  conducting  arms  spiraling  around  the  surface  of 


4-4 


o 


O 


o 


a  cone.  The  logarithm  of  the  distance  of  the  arms  from  the  apex,  measured  along  the  cone 
surface,  is  proportional  to  the  angle  of  rotation  of  the  arm  about  the  axis.  Figure  4-1  shows 
results  of  near  field  measurements  for  a  typical  log  spiral  antenna  of  cone  half-angle 
10  degrees  and  80  degrees  pitch.  The  graph  shows  electric  field  amplitude  at  a  distance  of 
0,03  A  from  the  cone  surface  versus  distance  from  the  apex.  In  this  figure  it  is  clear  that 
the  antenna  has  an  active  region  of  finite  extent,  which  for  the  example  shown  is  0.23^  wide 
at  the  3  db  points.  If  this  active  region  is  assumed  equivalent  to  a  slot  of  the  same  width, 
then 


*  -  14  A  (4-6) 

Based  on  the  results  of  this  section  it  is  probably  reasonable  to  assume  that  the  ration  tyA 
takes  values  from  2n  to  about  14  for  most  types  of  antennas.  The  larger  values  are  appli¬ 
cable  to  complicated  antennas,  such  as  periodic  structures,  and  the  smaller  values  apply 
to  simple  slot  antennas. 

An  interesting  example  of  a  fairly  rigorous  calculation  of  near-field  and  diffusion  length  of 
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a  spheroidal  monopole  is  provided  by  Cullen  and  Dobson  ,  who  use  the  power  law  approxi¬ 
mation  to  U/D. 

V/D  =  <IV D)q  (E/e/  (4-7) 

Their  result  is 


A  -  R/a(0  +  l)  (4-8) 

where  R  is  the  r  adius  of  the  antenna  tip  and  a  is  the  root  of  the  Bessel  function. 


2  (.6+1) 


(a) 


0 


(4-9) 
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SECTION  5 

UNEAR  PLASMA  EFFECTS 


An  antenna  radiating  through  the  reentry  induced  body  plasma  will  be  subject  to  certain 
effects  regardless  of  how  low  the  power  level  input  to  the  antenna  is  made.  These  effects 
which  tend  not  to  depend  on  the  magnitude  of  the  local  electromagnetic  field  strength  are 
termed  linear  plasma  effects.  Such  effects  may  take  various  forms.  Absorption  and  re¬ 
flection  of  electromagnetic  waves  by  the  plasma  results  in  attenuation  of  the  transmitted 
signals.  Antenna  pattern  distortion  causes  the  apparent  attenuation  to  depend  on  the  direction 
and  polarization  of  the  receiver.  The  energy  reflected  by  the  plasma  couples  into  the 
antenna  and  causes  its  input  impedance  to  change  from  the  value  it  would  have  in  free  space. 
The  change  in  input  impedance  can  cause  a  mismatch  to  the  device  which  feeds  the  antenna 
and  a  reduction  in  the  power  delivered  to  it.  Severe  mismatch  may  even  create  nonlinear 
circuit  effects  in  the  transmitter. 

The  theoretical  analysis  of  linear  plasma  effects  is  made  tractable  by  simplified  description 
of  the  plasma  constitutive  parameters,  which  work  well  in  the  types  of  plasmas  created  by 
reentry  phenomena.  The  propagation  of  transverse  electromagnetic  waves  through  such 
plasmas  can  be  analyzed  fairly  rigorously  with  comparative  ease,  for  simple  plasma 
geometries.  Such  analyses  are  ob  dously  correct  for  the  radiation,  or  far  zone,  fields  of 
an  antenna;  but  the  presence  of  the  reentry  plasma  in  the  antenna  near-field  complicates  the 
problem.  Thus,  although  the  wave  propagation  theory  is  capable  at  times  of  giving  fairly 
correct  results ,  the  near- field  problem  of  antenna-plasma  interaction  should  actually  be 
treated  in  any  mathematical  analysis.  For  example,  reactive  near  fields  can  be  made 
dissipative  by  the  plasma. 

5.1  SURVEY  OF  THEORY 
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A  recent  book  by  Shkarofsky,  Johnston,  and  Bachynski  gives  a  fairly  thorough  treatment 
of  the  theory  of  plasma  interaction  with  electromagnetic  fields.  Basic  to  all  linear  plasma 
effects  theory  is  the  cold  plasma  approximation,  in  which  all  effects  of  thermal  motion  of 
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the  electrons  are  ignored.  This  approximation  is  well  justified,  since  the  thermal 
velocity  is  much  less  than  the  phase  velocity  of  electromagnetic  waves. 


The  electromagnetic  field  interaction  is  determined  by  the  plasma  conductivity  a, 
which  has  the  following  form  in  sufficiently  weak  magnetic  fields. 


I<2  ' 

"  e0  i>  ^c  ‘  iu,) 


(5-1) 


where  is  the  permittivity  of  free  space,  u»p  is  the  radian  plasma  frequency, 
is  the  momentum  transfer  collision  frequency*,  and  is  the  radian  frequency  of  the 
wave  which  is  assumed  to  have  exp  (-iajt)  time  dependence. 


ufc  =  ^ne/(eom)  =  5.646xl04^rT  (5-2) 

where  n  is  tie  electron  density,  e  is  the  electron  charge,  m  is  the  electron  mass,  and 

-3 

the  numerical  constant  applies  when  n  is  expressed  in  cm  .  The  electric  field  vector  E 
in  a  plasma  satisfies  the  wave  equation. 


7*7xE  =  k2(l-iff/coWj  E  (5-3) 

where  k  =  id/c,  and  c  is  the  speed  of  light.  Thus  propagation  in  a  given  plasma  can 
be  analyzed  in  terms  of  the  parameters  electron  density,  collision  frequency,  and  signal 
frequency. 


*  The  use  of  collision  frequency  as  a  constant  in  the  expression  for  the  electrical 
conductivity  is  not  rigorously  possible,  but  an  effective  constant  collision  frequency 
is  reasonably  accurate  in  most  cases.  See  Section  7. 


The  fact  that  the  reentry  plasma  around  a  slender  cone  vehicle  is  confined  to  a  relatively 

thin  boundary  layer  can  be  used  to  advantage  in  solving  the  propagation  problem*  The 

boundary  conditions  may  be  simplified  to  the  requirements  that  the  tangential  electric 

field  be  continuous  and  the  discontinuity  in  the  tangential  magnetic  field  be  equal  to 
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the  surface  current  density  in  the  sheath.  For  example,  Poeverlein  has  solved 
the  problem  of  a  plane  wave  incident  on  a  plane  thin-plasma  sheath.  His  results  are 
a  function  of  the  surface  density  N  of  electrons,  and  their  form  depends  on  the  wave 
polarization.  For  polarization  in  the  plane  of  incidence  (E  and  the  plasma  normal  are 
coplanar  with  the  direction  of  propagation),  the  volta&o  reflection  and  transmission 
coefficients  R  and  T  are 


— T  r  NX  cos  0 
o  _ 

■'c 


(5-4) 


T  = 


1 

r^NX  cos0 

1  +  V  /<*> " 1 
c 


(5-5) 


2  2  -13 

where  the  classical  electron  radius  r  *  e  /(4s-  e  me  )  =  2. 82  x  10  cm,  A  Is  the  wavelength , 

o  o 

and  0  Is  the  angle  of  Incidence  with  respect  to  the  normal.  For  perpendicular  polarization. 


-T  r  NX 

o _ 

(v  /lO“i)  cos  0 
c 


(5-6) 


T  = 


_ 1 _ 

r  NX 

1  +  -2 - 

/Ct>  -i)  COB  0 
c 


(5-7) 


Even  in  the  simplication  of  the  thin  plasma  approximation,  which  Is  valid  when  the  surface 
electron  density  divided  by  the  peak  volume  electron  density  is  less  than  1/k,  the  geometry 
of  the  reentry  plasma  presents  difficulties  in  terms  of  orthogonal  coordinate  systems.  But 
it  is  necessary  to  bring  into  play  the  geometry  of  either  the  plasma  or  the  antenna  in  order 
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to  go  beyond  the  simple  attenuation  and  reflection  effect.  The  above  results  are  very 

unrealistic  with  regard  to  predicting  antenna  pattern  distortion,  and  they  are  only  qualitatively 

useful  in  the  estimation  of  antenna  impedance  changes.  The  major  advantage  of  these 

equations  is  in  terms  of  estimating  the  relative  importance  of  a  given  electron  density  level 

or  in  terms  of  defining  the  level  of  electron  density  at  which  plasma  effects  become  important. 

Thus  the  surface  density  of  electrons  at  which  plasma  effects  are  significant  is  N  =  (rQX)  . 

Of  course  In  a  thick  plasma  the  criterion  is  w  =  w. 

P 

A  higher  level  of  sophistication  involving  the  effects  of  the  plasma  geometry  on  antenna  pattern 
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distortion  is  provided  by  the  work  of  Banos,  et  al.  .  These  authors  solved  the  problem  of 

a  thin  conical  plasma  sheath  surrounding  a  point  dipole,  by  transformation  of  the  normal  mode 

solutions  to  the  form  of  a  residue  series.  The  roots  of  the  dispersion  relation  which  is 

obtained  indicate  the  contributions  of  surface  waves  to  the  radiated  fields  and  the  pattern 

distortion.  These  calculations  also  have  been  generalized  to  include  the  effects  of  collisions, 
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by  Jordan  and  Pinski  .  The  results  show  that  the  effects  of  the  conical  plasma  sheath  are 
significant  at  ne»  ly  the  same  values  of  surface  plasma  density  as  predicted  by  the  above  plane 
geometry  theory,  but  that  in  the  range  when  the  plasma  has  an  intermediate  effect  on 
attenuation  it  may  have  a  large  effect  on  antenna  radiation  pattern.  When  either  the  electron 
density  or  the  collision  frequency  becomes  large,  the  pattern  distortion  is  again  minimum. 

At  very  large  collision  frequencies  the  patterns  show  very  little  attenuation  or  distortion,  in 
qualitative  agreement  with  the  plane  geometry  theory. 

It  is  physically  feasible  for  the  effect  of  the  plasma  on  antenna  pattern  distortion  to  be  calculated 
fairly  accurately  without  including  near  field  effects,  since  the  far  field  pattern  can  be 
significantly  affected  by  surface  waves  in  the  plasma  at  relatively  large  distances  from  the 
antenna.  The  radiation  from  these  surface  wave  modes  can  dominate  the  far  field,  while 
the  source  of  primary  waves  is  localized  to  the  antenna's  effective  aperture  and  may  not 
appreciably  change  the  radiation  field  geometry  in  some  cases.  In  fact  the  effects  of  plasma 
geometry  and  various  other  obstacles  remotely  located  can  have  predominant  effects.  For 
this  reason  also,  the  pattern  distortion  is  not  a  good  diagnostic  indication  of  the  plasma  density 
over  the  antenna. 
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Near-field  effects  must  be  taken  into  account  in  the  theory  of  plasma  effects  on  antenna 
impedance.  In  this  area  of  the  problem,  the  major  subject  of  study  has  been  the  slot  antenna. 
The  infinite  slotted  cylinder  has  been  studied  by  Smith  and  Golden226,  Swift227,  and  Rusch228. 
and  it  represents  a  convenient  model  for  mathematical  analysis  because  of  the  separability  of 
the  wave  equation  in  cylindrical  coordinates.  Much  work  has  been  done  recently  on  finite 
slots  in  ground  planes  covered  by  plasma  layers.  Villeneuve229  has  computed  the  effective 
input  impedance  of  a  rectangular  waveguide  terminated  in  an  infinite  ground  p  lane,  where  the 
ground  plane  is  covered  by  a  dielectric  slab  of  finite  thickness.  For  a  plasma  slab,  his  results 
show  the  impedance  or  its  reciprocal,  admittance,  to  be  relatively  insensitive  to  the  plasma 
thickness,  when  the  thickness  is  greater  than  0.  625  times  the  larger  waveguide  or  aperture 
dimension.  This  analysis  does  mt  include  the  overdense  plasma  case,  /here  co  /o>  >  1. 

230  231 

Croswell,  et  al,  have  considered  the  admittance  of  a  rectangular  waveguide  terminated 

in  an  infinite  ground  plane  covered  with  an  inhomogeneous  plasma,  both  theoretically  and 
experimentally.  The  experimental  data  is  that  of  Taylor,  obtained  at  Stanford  Research 
Institute.  The  derivation  of  the  admittance  follows  essentially  the  same  method  of  Villeneuve 
but  is  extended  to  electron  densities  greater  than  the  critical  density,  i.  e. ,  w  /w>l.  The 
experimental  data  was  obtained  using  both  shock-heated  and  rf  generated  plasmas  flowing 
past  the  aperture  in  a  ground  plane.  Various  boundary  layer  thicknesses  or  electron  density 
gradients  were  assumed  for  the  calculations  and  the  results  show  that  the  susceptance  is 
strongly  dependent  on  the  electron  density  gradients  at  the  aperture.  Reflection  coefficient 
measurements  for  the  rf -gene rated  plasma  agreed  quite  well  with  the  computed  values  up  to 
about  critical  density,  but  showed  a  wide  disagreement  above  critical.  This  was  due  to  the 
conductance  not  approaching  zero  above  critical  electron  density,  as  predicted. 

The  admittance  of  a  rectangular  aperture  radiating  into  plasma  slab  has  also  been  extensively 
considered  by  Galejs232  and  Galejs  and  Mentzoni233.  This  analysis  differs  from  those  of 
Villeneuve  and  Croswell  in  that  the  waveguide  radiates  into  a  larger  rectangular  waveguide 
Instead  of  free  space.  The  computed  admittance  was  compared  with  measurements  using 
the  afterglow  of  a  capacitor  discharge  as  the  plasma,  and  relatively  good  agreement 
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resulted.  It  is  apparent  that  antenna  admittance  (or  impedance)  is  more  readily  related 
to  the  properties  of  the  plasma  in  the  immediate  vicinity  of  the  antenna  than  are  the  parameters 
of  the  radiation  field.  This  fact  is  the  reason  for  the  relative  success  of  plasma  diagnostics 
by  impedance  measurement. 

An  additional  effect,  that  of  mutual  coupling  between  antennas,  has  been  investigated  by 
232  ,  233 

Galejs  and  Galejs  and  Mentzoni  .  Generally  speaking  the  presence  of  a  plasma  sheath 
tends  to  decrease  the  mutual  coupling  between  slot  antennas  in  a  ground  plane. 

5.2  ILLUSTRATIVE  CALCULATIONS 

The  discussion  above  indicates  that  the  major  impact  of  antenna  near-field/linear  plasma 
interactions  is  obtained  in  terms  of  the  electrical  performance  of  the  antenna  as  a  circuit 
element.  The  antenna/transmitter  system  design  thus  has  a  significant  effect  on  the 
importance  of  the  near-field  effects.  Assuming  that  the  system  is  designed  to  account  for 
such  effects,  by  the  use  of  broadband  components  and  sufficient  antenna  isolation,  then  the 
radiated  field  characteristics  are  of  primary  concern.  Parametric  representation  of  near- 
field  effects  on  antenna  impedanco  is  difficult  to  accomplish,  at  any  rate,  because  of  its 
dependence  on  a  multitude  of  parameters. 

Radiation  pattern  distortion  is  also  difficult  to  treat  parametrically,  but  it  is  generally 

comparable  in  importance  to  the  attenuation  effect  as  an  indication  of  electron  density  levels 

which  have  an  appreciable  effect  on  received  signals.  Thus  the  simple  attenuation  calculation, 

which  includes  losses  due  to  both  absorption  and  reflection,  provides  a  fairly  useful 

illustrative  calculation.  Figures  5-1  and  5-2  show  the  total  attenuation  (A  - 10  log  1 1/T  j 2) 

given  by  Equation  (5-5),  for  normal  incidence.  At  small  values  of  the  collision  frequency, 

the  attenuation  depends  on  the  parameter  NX  ,  while  at  large  values  of  collision  frequency, 

the  parameter  NX  enters,  where  X  =  2jt  c/v  .  These  curves  may  be  used  to  represent 
c  c  c 

either  polarization  at  arbitrary  incidence  angle  by  properly  normalizing  to  cos  0  or  sec  0  ; 
however  the  application  of  this  theory  to  antenna  pattern  distortion  is  not  realistic. 
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The  collision  frequency  v  generally  increases  with  decreasing  altitude,  at  a  rate  dependent 

c 

on  the  flow  field  properties.  A  typical  situation  for  a  9  degree  half-angle  cone  reentering  with 

a  velocity  of  22, 500  ft/sec  is  shown  in  Figure  5-3.  These  results  were  obtained  using 
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Musal's  tables  for  equilibrium  air  collision  frequency,  and  they  show  that  the  collision 
frequency  traverses  the  microwave  region  in  the  same  range  of  altitudes  as  that  in  which  the 
electron  density  becomes  significant  on  a  slender  cone  vehicle. 

These  results  show  that  the  linear  plasma  effects  on  signal  transmission  from  a  slender  cone 
vehicle  become  independent  of  signal  frequency  and  tend  to  decrease  (p  ovided  the  electron 
density  does  not  Increase)  as  the  altitude  decreases,  once  the  altitude  is  less  tha*:  that  at 
which  the  collision  frequency  equals  the  signal  frequency. 
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SECTION  6 

ANTENNA  BREAKDOWN 


The  phenomenon  of  antenna  voltage  breakdown  becomes  an  important  consideration  whenever 
relatively  high  power  levels  are  to  be  transmitted,  especially  when  the  atmosphere  surrounding 
the  antenna  is  reduced  to  pressures  on  the  order  of  one  Torr.  If  the  power  input  to  an  antenna 
is  increased  from  low  levels,  the  output  power  initially  Increases  proportionately,  but  eventu¬ 
ally  a  sudden  decrease  in  output  power  will  result  from  breakdown  of  the  air  surrounding  the 
antenna.  The  strong  electric  fields  near  the  antenna  create  ionization  in  the  air  at  a  greater 
rate  than  the  loss  processes  can  remove  the  electrons  from  this  high  field  region.  This  field- 
induced  plasma  reflects  and  absorbs  the  incident  electromagnetic  energy  so  that  very  little 
power  is  transmitted.  The  value  of  input  power  at  which  breakdown  takes  place  is  called  the 
breakdown  power,  and  it  is  obtained  reproducibly  for  a  given  antenna  at  a  given  pressure 
under  the  proper  conditions,  including  the  presence  of  initiating  electrons. 

Breakdown  power  generally  varies  with  pressure  in  such  a  way  as  to  undergo  a  minimum  at 
a  particular  pressure.  At  low  pressure,  the  electrons  diffuse  rapidly  away  from  the  high 
field  region,  causing  the  breakdown  power  to  increase.  At  high  pressure,  the  increased  rate 
of  collisions  of  electrons  with  neutral  particles  reduces  the  energy  gained  by  the  electrons 
from  the  field,  also  causing  the  breakdown  power  to  increase.  (The  high  pressure  region  in 
air  is  often  referred  to  as  the  attachment  controlled  region,  since  attachment  to  neutral 
molecules  is  the  dominant  electron  loss  mechanism. )  Both  the  minimum  breakdown  power 
and  the  pressure  at  the  minimum  tend  to  increase  as  the  signal  frequency  is  increased. 

The  subject  of  antenna  breakdown  can  be  related  to  the  subject  of  high  frequency  gas  dis¬ 
charges  or  microwave  cavity  breakdown,  for  general  theoretical  description  of  the  processes 
involved.  Hie  advantages  of  this  approach  derive  from  the  extensive  study  of  the  latter  field 
which  has  been  carried  out  during  the  nast  twenty  years.  The  terms  and  equations  used  in 
microwave  breakdown  theory  will  be  employed  in  the  discussion  of  antenja  breakdown. 


The  matter  of  principal  concern  In  this  report  Is  the  effect  of  the  reentry  environment.  Indeed 
there  would  be  no  need  to  discuss  the  subject  of  antenna  breakdown  were  it  not  for  the  effects 
of  reentry.  One  can  most  reliably  determine  the  breakdown  characteristics  of  an  a.itenna 
experimentally,  except  for  the  effects  of  a  given  reentry  environment.  The  important 
question  then  is  how  to  predict  the  change  in  the  breakdown  characteristics  of  an  antenna  in 
the  reentry  environment.  Before  providing  the  solution  to  this  problem,  we  review  the 
classical  theory. 

6.1  SURVEY  OF  CLASSICAL  BREAKDOWN  THEORY 

The  subject  of  microwave  breakdown  in  gases  has  been  thoroughly  surveyed  by  MacDonald  in 
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a  recent  book  .  This  book  presents  practically  all  of  the  basic  theory  of  high  frequency 
gas  discharge'..  Moreover,  it  indicates  how  breakdown  field  strengths  can  be  calculated 
reliably.  The  book  is  thus  useful  in  providing  the  basic  physical  theories  on  the  subject  as 
well  as  the  practical  tools  needed  for  applications.  In  addition,  a  whole  chapter  is  devoted  to 
experimental  methods. 

It  is  instructive  to  begin  the  discussion  of  breakdown  theory  with  a  qualitative  description  of 
the  breakdown  process.  The  impressed  electromagnetic  field  exerts  force  on  any  free  electrons 
which  may  be  present  in  the  gas.  Through  the  process  of  collisions  of  electrons  with  neutral 
gas  particles,  the  field  does  work  on  the  electrons,  raising  their  average  temperature  well 
above  that  of  the  gas.  Ionization  of  neutral  particles  by  impact  of  high-energy  electrons 
competes  with  electron  attachment  to  neutral  particles  to  give  a  net  rate  of  electron  population 
gain.  This  nen  ionization  rate  competes  with  losses  by  diffusion  until  the  point  is  reached 
where  the  ionization  builds  up  at  a  very  rapid  rate  and  the  gas  becomes  conducting.  The 
various  collision  processes  are  very  important  in  determining  breakdown,  in  that  they  control 
the  energy  distribution  as  well  as  the  number  density  of  electrons. 

Generally  speaking  there  are  two  methods  of  analysis  which  can  be  used  to  predict  breakdown, 
the  kinetic  theory  method  and  the  phenomenological  theory. 
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'Hie  kinetic  theory  approach  to  breakdown  prediction  treats  the  electron  distribution  function 

F  (v,  r,  t),  which  represents  the  number  of  electrons  with  velocities  close  to  v,  in  a  small 

volume  of  space  at  the  point  r,  and  at  the  time  t.  This  distribution  function  obeys  the 
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Boltzmann  equation  where  v  is  the  gradient  in  configuration  space,  %  is  the  r,radient  in 
in  velocity  space,  a  is  the  acceleration,  and  C  represents  the  effects  of  collisions. 


S  F  _ 
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(6-1) 


The  collision  teixn  C  may  be  written  as  an  integral  involving  the  distribution  function  and  the 
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collision  cross  section  ,  and  hence  knowledge  of  collision  cross  section  as  a  function  of 
electron  velocity  is  sufficient  in  principle  to  solve  the  Boltzmann  equation.  However,  the 
difficulties  involved  are  so  great  fhat  accurate  results  are  possible  in  only  a  few  cases. 

Hydrogen  and  helium  have  practically  constant  collision  frequencies  above  an  energy  of  about 
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3  or  4  electron  volts,  and  this  fact  has  been  used  by  MacDonald  and  Brown  *  and  Reder 
240 

and  Brown  to  calculate  breakdown  fields  for  these  two  gases  and  for  admixtures  of  mercury. 
The  results  agree  very  well  with  experiment,  expecially  for  the  case  of  helium  containing  a 
small  amount  of  mercury.  The  presence  of  mercury  serves  to  convert  the  metastable 
excitation  level  of  helium  to  ionization  of  mercury  by  an  exchange  reaction.  This  has  the 
effect  that  the  gas  has  no  excitation  levels  below  the  effective  ionization  level,  greatly  simpli¬ 
fying  the  analysis.  The  theoretical  results  for  neon  and  neon-argon  mixtures,  for  which  the 
collision  frequency  is  approximately  proportional  to  the  square  root  of  the  electron  energy, 
do  not  agree  with  experimental  data  quite  as  well  as  for  constant  collision  frequency  gases. 

Of  course  the  use  of  simple  collision  frequency  relationships  for  predicting  atmospheric 
breakdown  by  kinetic  theory  does  not  work.  * 

*  Lenander  and  Epstein241  have  reported  recently  on  the  successful  Drediction  of  breakdown 
fields  of  air,  using  a  computer  piogram  of  Carelton  and  Megill,242  with  detailed  collisional 
data. 
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The  phenomenological  method  of  predicting  breakdown  has  been  widely  used  by  various 
authors  in  classical  breakdown  theory.  Essentially  this  method  is  a  collection  of  extrapolation 
formulas  based  on  the  best  fit  of  reasonable  phenomenological  ideas  to  existing  basic  data. 

This  method  starts  with  the  electron  continuity  equation,  which  can  be  derived  from  the 
Boltzmann  equation.  The  various  terms  in  this  equation  are  then  treated  as  simple  functions 
of  the  proper  variables  of  the  problem,  and  experimental  data  are  used  to  determine  the 
artibrary  constants  in  these  equations. 

The  electron  continuity  equation  is  written  as 


7-  f  + 


dn  dn 
at  dt 


(6-2) 


where  n  is  the  electron  density,  T  is  the  electron  flux,  and  dn/dt  is  the  net  effect  of  all 
sources  and  sinks.  Electron  diffusion  leads  to  the  electron  flux 

f  =  -  7  (Dn)  (6-3) 

where  r  is  the  diffusion  coefficient.  When  the  net  ionisaticn  rate  per  electron  is  represented 
by  the  symbol  v,  the  continuity  equation  may  be  written  as 


dn 

at 
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This  equation  is  generally  solved  by  assuming  separability,,  using 
72  Pn)  =  -Dn/A2 

where  A  is  the  characteristic  diffusion  length.  Thus 
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(6-6) 


-^-=  P  -D/A2)n 

The  Townsend  criterion,  carried  over  from  the  theory  of  direct  current  discharges,  has  been 
shown  by  Berlin  and  Brown  to  apply  to  microwave  discharges.  According  to  this  criterion, 
the  breakdown  point  is  reached  when  the  right-hand  side  of  Equation  6-6  is  zero,  since  a 
slightly  higher  electric  field  would  cause  the  runaway  production  of  electrons. 

When  either  the  electric  field  or  the  electron  mobility  is  not  uniform  in  space,  as  would  be  the 
case  in  antenna  breakdown,  the  diffusion  equation  must  be  solved  in  the  form  of  Equation  6-4  , 
with  the  left-hand  side  set  equal  to  zero.  However,  it  is  possible  to  use  an  effective  diffusion 
length  which  can  be  derived  from  the  solution  of  the  diffusion  equation.  For  example,  Herlin 

pi  c 

and  Brown  derived  the  solution  for  a  -mode  cylindrical  cavity  and  MacDonald  and 

Brown216  derived  the  results  for  a  spherical  cavity.  (See  Section  4  for  a  discussion  of  methods 
of  estimating  the  diffusion  length  for  antennas. ) 
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It  can  be  shown  ,  through  certain  considerations  in  the  kinetic  theory,  that  the  important 
parameter  with  respect  to  electric  field  E  and  the  microwave  radian  frequency  a>  is 

E2  =  v  2E2/(u  2  +  J)  (6-7) 

6  C  C 

where  E  is  the  effective  field  and  v  is  the  collision  frequency  for  momentum  transfer.  The 
e  c 

effective  field  concept  is  useful  in  the  phenomenological  theory  of  breakdown,  as  in  the 
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analysis  of  Brown  ,  Rose  and  Brown  ,  and  Gould  and  Roberts  .  The  Townsend 

criterion  is  written  as 

v =  v  +  D/A2  (6-8) 

i  a 
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where  and  v  are  the  ionization  and  attachment  frequencies,  respectively.  The  first 
Townsend  coefficient  a  is  used  in  place  of  p  according  to  the  equation 


1/  =  0)1  E 

i  e 


where  p  is  the  dc  mobility.  The  ratio  of  diffusion  coefficient  to  mobility  is  invoked  to  give 
the  relation 


or.  _  £  2  u 

p  p  3  (pA)  (E  A) 


(6-10) 


where  p  is  pressure,  u  is  average  electron  energy,  and  0  is  the  number  of  attachments  per 

centimeter  of  path  length.  Data  exist  for  a/p,  0/p,  and  u  as  functions  of  E/p  in  air244, 24?. 

Hence,  Equation  6-10  may  be  used  to  give  E  /p  as  a  function  of  p  A ,  provided  E  can  be 

e  e 

calculated. 


The  difficulty  connected  with  the  phenemenological  theory  is  that  the  effective  field  is  not  a 
sound  concept,  since  the  collision  frequency  is  not  constant.  Brown244  obviates  this  difficulty 
by  calculating  an  effective  constant  collision  frequency  from  the  measured  ac  mobility,  using 


Mac  =  (e/m)  vc/{Uc  +  w 2) 


(6-11) 


where  e/m  is  the  ratio  of  charge  to  mass  for  the  electron.  He  gets  a  value  of  v  equal  to 
9  c 

4,3  x  10  p.  The  results  of  breakdown  calculations  then  agree  with  experimental  data  fairly 

well  over  a  limited  range  of  pressures. 


It  appears  that  the  most  reliable  and  generally  useful  method  of  calculating  breakdown  is  that 
given  by  MacDonald.  This  method  uses  the  effective  field  concept  in  a  way  which  is  much  less 
sensitive  to  the  specific  form  of  collision  frequency  assumed  than  is  the  above  method.  The 


diffusion  coefficient  can  be  calculated  from  the  variation  of  collision  frequency  with  electron 
energy,  if  the  electron  energy  distribution  is  assumed;  however,  it  can  be  shown  that  the  form 
of  the  distribution  function  has  little  effect  on  the  result.  Hence  the  approximation  of 
MacDonald,  Gaskell,  and  Gittem'an248  works  well. 

Dp  =  (29  +  0.9Ee/p)  lo4  cm2  -Torr/sec  (6-12) 


This  equation  was  derived  from  assuming  a  linear  dependence  of  Dp  on  average  electron 
energy  in  terms  of  the  latter's  experimentally  determined  dependence  on  E/p.  Equation 
6-12  is  equivalent  to 


where  X  is  the  wavelength  and 
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where  EX  is  in  volts  and  pX  is  in  Torr/cm. 

The  net  ionization  rate  *  may  be  related  to  the  electric  field  and  the  pressure  by  using 
experimental  breakdown  data.  The  result  is  shown  in  Figure  8-18  of  Reference  235  in  terms 
of  VX  versus  pX  for  fixed  values  of  EX.  Figure  8-19  in  Reference  235  shows  a  plot  of 
Equation  6-14  .  The  solution  to  a  given  problem  in  continuous  wave  breakdown  in  air  is 
found  by  the  points  of  intersection  of  these  two  sets  of  curves,  as  shown  in  Figure  8-20  in 
Reference  235,  since  at  breakdown,  vX  is  equal  to  DX/A2 .  Figure  6-1  shows  curves  of  EX 
versus  pX  for  three  different  values  of  X/A ,  derived  from  MacDonald's  curves. 
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The  above  method  of  calculating  continuous  wave  breakdown  can  be  generalized  to  include 
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pulse  (or  finite  time)  effects,  following  the  approach  of  Gould  and  Roberts  .  For  example, 
for  pulsed  breakdown  at  low  pulse  repetition  rate,  the  net  ionization  frequency  is  determined 

ty 


•'b  T  =  lu  W  (6_15) 

where  T  is  the  pulse  length,  n.  is  the  electron  density  required  for  breakdown,  and  n  is  the 

“  13  °2 

initial  electron  density.  The  value  of  n^  is  usually  taken  as  the  critical  value,  10  /A  , 

(where  X  is  in  cm  and  n^  is  cm  ^  and  nQ  is  a  fairly  small  number.  For  a  given  pulse  length, 

Equation  6-15  determines  v,  ,  which  is  actually  related  to  v  of  the  CW  analysis  by 

b 

V  X  =  DA/A2  +  v.  X  (6-16) 

b 

The  procedure  is  to  start  with  EX  as  a  function  of  pX  from  continuous  wave  theory  and  to 
2 

determine  DX/A  from  Equation  6-13.  Then  v  X  is  determined  from  Equation  6-16.  This 
gives  EX  from  MacDonald's  curves,  and  the  process  can  be  iterated  to  determine  EX  more 
accurately.  In  the  case  of  rapid  pulse  repetition  rates,  the  decay  time  of  the  afterglow 
between  pulses  must  be  taken  into  account.  Then  the  breakdown  criterion  i  j  that  the  buildup 
of  ionization  during  the  pulse  is  just  equal  to  the  decay  between  pulses.  A  slightly  greater 
net  ionization  rate  during  the  pulse  would  give  a  runaway  condition.  As  the  pulse  repetition 
rate  is  increased,  the  breakdown  field  approaches  the  continuous  wave  case. 

It  remains  for  us  to  show  the  limits  of  validity  of  the  above  phenomenological  theory,  referred 
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to  as  the  diffusion  theory  of  microwave  breakdown.  Brown  and  MacDonald  have  studied 
this  question  for  hydrogen,  and  their  approach  can  be  applied  to  air  with  relatively  insignificant 
changes  in  numerical  results.  Briefly,  the  principal  physical  limitations  on  the  validity  of  the 
theory  have  to  do  with  the  fact  that  the  electrons  must  remain  in  the  high  field  region  for  more 
than  a  single  collision  or  a  single  cycle  of  the  alternating  field.  When  the  mean  free  path 
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equals  the  diffusion  length,  (p  A^o.  03  cm  Torr),  we  have  the  mean  free-path  limit.  At  lower 
pressures  the  diffusion  theory  is  inapplicable.  When  the  amplitude  of  the  osclllatoiy 
motion  of  an  electron  in  the  field  exceeds  the  diffusion  length,  we  have  the  oscillation  amplitude 

Q 

limit  (X/  A  is  on  the  order  of  10  or  greater),  beyond  which  the  theory  is  also  inapplicable. 

The  uniform  field  limit  (X  =  2  tt  A )  need  not  be  considered  as  a  boundary  of  the  region  of 
validity,  provided  the  field  distribution  is  accounted  for  in  the  diffusion  length. 

6.2  SURVEY  OF  REENTRY  EFFECTS 

The  subject  of  this  section  is  the  effects  of  reentry  on  antenna  breakdown,  w  ith  emphasis  on 
the  reentry  environment  associated  with  slender  cone  bodies  at  hypersonic  speeds.  The 
discussion  will  depend  on  ideas  and  data  taken  from  written  reports  and  verbal  communications 
with  a  number  of  researchers  in  the  field,  in  the  manner  of  a  survey  of  present  knowledge. 

For  conical  vehicle  reentry,  three  distinct  flow-field  regions  may  be  important  in  antenna 

breakdown:  the  ambient  free  stream,  the  in  viscid  shock  layer,  and  the  boundary  layer.  With 

the  exception  of  relatively  high  altitudes  where  these  regions  are  not  necessarily  distinct,  the 

reentry  environmental  effects  on  breakdown  will  be  different  for  each  of  the  three  regions.  In 

fact,  most  workers  assume  that  a  convenient  separation  can  be  made  such  that  the  effects  of 

each  region  can  be  treated  independently  .  At  any  rate,  the  different  basic  processes  which 

control  breakdown  can  be  identified  for  each  of  the  three  flowfield  regions,  for  example,  as 
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d  ne  by  Reilly  and  numerous  other  authors. 

The  free  stream  ahead  of  the  shock  layer  is  characterized  essentially  only  by  the  fact  that  the 
air  is  moving  relative  to  the  transmitting  antenna.  The  effects  of  altitude  on  the  state  or 
composition  of  the  ambient  air  probably  are  not  significant  in  breakdown;  and  even  if  they  were, 
it  would  be  a  relatively  simple  task  to  account  for  them.  Hence,  the  only  important  effect  in 
the  free  stream  is  convection. 


The  inviscid  flow  region  is  characterized  by  a  state  of  elevated  temperature  and  increased 
molecular  number  density  relative  to  the  free  stream.  On  a  slender  conical  vehicle,  ionization 
is  probably  insignificant  in  this  region,  but  the  detailed  thermodynamic  state  and  chemical 
composition  of  the  shock-heated  air  may  be  quite  important.  At  any  rate,  the  effects  of  the 
inviscid  flow  region  are  felt  in  terms  of  convection,  increased  air  density,  and  elevated 
temperature. 

The  boundary  layer  is  characterized  by  even  higher  temperatures  produced  by  viscous 
dissipation.  The  presence  of  significant  levels  of  thermal  ionization  Is  thus  probable. 

Particle  number  density  and  convection  are  not  so  easy  to  generalize  because  of  the  strong 
gradients  of  these  quantities  which  exist  here.  In  addition,  the  presence  of  chemical 
contamination  of  the  air  by  products  of  ablation  will  be  important.  Thus,  convection,  gas 
density,  temperature,  ionization,  and  chemical  additives,  all  in  the  context  of  strong 
gradients,  must  be  considered  in  the  boundary  layer. 

In  summary,  there  can  be  distinguished  six  phenomena  of  hypersonic  reentry  flow  fields 
which  can  affect  antenna  breakdown  in  one  or  more  of  the  regions  of  the  flow  field.  The 
following  six  subsections  treat  each  of  these  in  varying  detail,  depending  on  the  extent  of 
present  knowledge. 


6.2.1  GAS  DENSITY 

The  classical  theory  of  microwave  breakdown  is  developed  in  terms  of  proper  variables  such 
as  i>/p  and  Dp,  where  p  is  the  pressure  in  Torr.  These  quantities  generally  normalize  in 
this  way  because  of  the  dependence  of  collision  frequencies  on  number  density  of  gas  particles. 
In  voltage  breakdown  experiments  under  standard  laboratory  conditions,  the  practice  is  to 
express  the  state  of  the  gas  being  tested  in  terms  of  pressure,  since  this  is  a  directly 
measurable  quantity.  However,  pressures  in  the  reentry  induced  shock  layer  are  often 
orders  of  magnitude  higher  than  in  the  free  stream.  Also,  because  of  increased  gas  tempera¬ 
ture,  shock-layer  density  is  not  generally  the  same  as  in  the  free  stream.  Whitmer  and 
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MacDonald  have  noted  that  the  physically  important  parameter  for  breakdown  is  the  gas 
density  rather  than  pressure,  since  it  is  the  number  of  gas  particles  per  unit  volume  which 
determines  the  rates  for  all  the  processes  which  control  breakdown.  Thus  the  use  of  either 
the  free- stream  pressure  or  the  shock-layer  pressure  in  predicting  breakdown  at  a  given 
altitude  would  give  erroneous  results.  The  distribution  of  gas  density  in  the  high  field 
intensity  region  is  the  correct  determining  factor. 

In  order  to  give  a  quantity  which  is  commensurable  with  the  classically  used  parameter  of 
pressure,  most  authors  define  a  reduced  pressure  as  p*.  Reilly  defines  the  reduced  pressure 
in  Torr  in  terms  of  the  gas  number  density  N  in  particles  per  cubic  centimeter. 

p*  =  2.9x10  17  N  (6-17) 
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Light  and  Taylor  define  the  reduced  pressure  in  an  equivalent  way  in  terms  of  temperature, 
T,  in  °K: 


p*  *  273  p/T  (f  '  ) 

These  corrections  have  the  effect  of  shifting  the  curve  of  breakdown  field  versus  altitude 

about  30, 000  feet  higher  in  altitude  at  velocities  on  the  order  of  22, 000  ft/sec,  according  to 
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Whitmer  end  MacDonald 
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Gas  density  effect?  have  been  observed  in  reentry  flight  experiments.  Bisbing  showed  that 

the  observed  period  of  breakdown  of  a  slot  antenna  on  a  conical  vehicle  agreed  with  laboratory 

data  in  cold  ah'  when  the  air  density  corresponding  to  the  inviscid  flow  region  was  used  as 
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the  independent  variable.  Nanevicz  also  showed  that  when  the  VHF  breakdown  data  from 
Nlke-Cajun  flights  is  corrected  for  local  density,  it  agrees  roughly  with  laboratory  measure¬ 
ments  in  cold  air. 

The  gas  density  affects  all  the  terms  of  the  continuity  equation  in  their  proper  variables 
form  in  terms  of  the  reduced  pressure  p*  in  place  of  the  pressure  p. 
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6.2.2  IONIZATION 

A.,  appreciable  electron  deneity  probably  exists  In  the  flow  Held  of  a  reentry  body  at  all 
altitudes  of  Interest  for  antenna  breakdown.  This  Is  true  even  at  altitin.»s  too  high  lor  thermal 
Ionization  of  the  air.  because  the  ionosphere  provides  ambient  electrons.  This  initial 
ionization  has  various  effects  on  breakdown.  One  effect  of  some  practical  importance  is  that 
the  breakdown  threshold  will  be  obtained  fairly  often  for  a  given  ret  of  reentry  conditions. 
There  will  always  be  sufficient  ionization  to  initiate  breakdown,  which  requires  at  least  one 

free  electron  before  it  can  proceed. 

The  most  important  effect  of  flow-field  ionization  is  to  change  diffusion  from  free  ambi- 
poiar,  or  at  least  to  create  transitional  diffusion.  That  is,  with  the  buildup  of  a  siz^e 
electron  space  charge,  the  electron  and  molecular  diffusion  processes  become  inextricably 
linked.  The  result  is  to  slow  down  the  diffusion  of  electrons  to  'he  molecular  rate  and  ^ 
thereby  correspondingly  reduce  diffusion  as  an  electron  loss  mechanism.  Allis  and  Rose 
have  indicated  that  the  transition  from  free  to  ambipolar  diffusion  begins  where  the  Debye 
length  becomes  comparable  to  the  diffusion  length.  The  magnitude  of  the  transitional  dif¬ 
fusion  coefficient  D  can  be  calculated  only  by  solving  the  boundary  value  problem  consti¬ 
tuting  the  separate  electron  and  ion  continuity  equations,  with  space  charge,  and  Poisson  s 
equation.  The  result  copends  on  the  geometry  of  the  boundaries  as  well  as  on  the  ionization 

level. 

Allis  and  Rose255  have  perform  'd  calculations  for  h,  ’rogen  and  found  that  the  transition 

from  free  to  ambipolar  diffusion  takes  place  over  five  orders  of  magnitude  of  electron 

density.  Moreover  their  results  indicate  Uat  the  constant  ratio  approximation  is  nearly 

correct  if  the  electron  density  used  in  this  approximation  is  decreased  from  its  actual 

value  by  a  constant  factor.  (The  constant  ratio  approximation  consists  of  assuming  that  the 

ratio  of  electron  and  ion  densities  is  constant  throughout  the  discharge. )  From  Figure  4 

of  their  paper  it  appears  that  this  factor  should  be  about  20,  but  Allis  ’  indicates  a  factor 
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of  five  in  another  calculation.  A  factor  of  ten  may  be  reasonable.  Allis  has  derived 
an  expression  for  the  transitional  diffusion  coefficient  which  follows. 
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where  A  is  the  diffusion  length,  is  the  Debye  length,  p  is  mobility,  n  is  particle  density 

(+  refers  to  ions  and  -  to  electrons),  D  is  the  electron  free  diffusion  coefficient,  and  D  is 

a 

the  ambipolar  diffusion  coefficient. 
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where  f  q  Is  the  permittivity  of  free  space  and  e  is  the  electronic  charge.  In  practice,  the 
ratio  p+  n+/(p_n_)  is  generally  small  compared  with  unity,  except  possibly  when  the  Debye 
length  is  very  large,  in  which  case  the  ratio  is  not  important.  Thus,  a  good  approximation 
which  depends  only  on  the  electron  density  is  the  following  equation: 


,2  .2 

lD  +A 


D  j£  +  D  A2 
a  D 


(6-22) 


fa  thj  limits  of  low  and  high  electron  density, the  value  of  D^  approaches  D  and  D^  respectively, 
fa  order  to  allow  the  aoove  constant  ratio  approximation  formula  to  reproduce  approximately 
the  actual  transition  situation,  apply  the  factor  of  ten  recommended  above. 
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Thus  the  transition  from  free  to  ambipolar  diffusion  in  this  approximation  occurs  when  the 
Debye  length  is  about  one  third  of  the  diffusion  length. 

Two  simplified  approximations  to  the  ambipolar  diffusion  coefficient  are  possible,  considering 
the  fact  that  n  is  much  greater  than  n+.  In  an  active  discharge,  D_  p+  is  much  greater  than 
D+  p  ,  so  that  the  ambipolar  diffusion  coefficient  is  given  approximately  by  the  following 
equation: 


D^  =  (p+  /ju  )  D  (6-24) 

But  in  the  case  of  an  isothermal  plasma,  D  is  equal  to  D+  p  . 

D  =(2n  /n  )  D  =  2D  (6-25) 

It  is  important  to  note  that  the  latter  result  applies  only  to  the  afterglow  plasma,  for  example, 
the  decaying  plasma  which  exists  between  pulses  in  multiple  pluse  breakdown.  Also,  there  is 
not  actually  a  factor  of  two  difference  between  these  two  results,  because  D  /p  is  much 
greater  in  the  active  discharge  than  in  the  isothermal  plasma,  while  n+  is  practically  constant. 

It  is  apparent  that  one  needs  values  for  the  parameters  D  and  p+/p  to  apply  these  equations, 
and  also  one  needs  a  value  of  D  /p  in  order  to  calculate  the  Debye  length  for  Equation  6-23. 
The  latter  \  roblem  is  accurately  enough  handled  for  the  purposes  of  transitional  diffusion  in 
the  constant  ratio  approximation  by  casting  Equation  6-21  in  the  following  form: 

=  €  k  T  /ne^  (6-26) 

D  o  e 

where  k  is  Boltzmann's  constant  and  T  is  the  electron  temperature.  In  the  isothermal 

e 

afterglow  plasma  T  is  equal  to  the  gas  temperature,  and  in  an  active  discharge  in  air 

G 

kT  is  on  the  order  of  a  few  electron  volts.  At  breakdown  in  air,  a  commonly  used  value 
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of  D  is  given  by  Kelly  and  Margenau 

6  2  **1 

Dp  =  1.6x10  cm  Torr  sec  (6-27) 
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But  MacDonald  prefers  to  use  a  value  which  depends  on  the  microwave  field  intensity  E. 

10  *  D  p  *  29  +  0. 9  E  /p  cm2  Torr  sec  *  (6-28) 

e 

Various  values  of  p_/p+  (or  D  /D  )  for  breakdown  in  air  have  appeared  in  the  literature, 

25^  251 

including  Kelly  and  Margenau's  value  of  100  and  Whitmer  and  MacDonald's  value  of  40. 

But  these  are  not  actually  inconsistent,  since  the  latter  authors  applied  Equation  6-25  to  the 

data  on  mobilities  for  nonisothermal  plasmas.  Thus  the  value  of  40  is  actually  80  and  it 

seems  reasonable  to  assume  that  for  breakdown  of  air, 

D  p 

•—  =jL-  e=  100  (6-29) 

a  + 

Flow-field  ionization  can  also  affect  pulsed-signal  breakdown  in  terms  of  the  initial  electron  den 
sity  in  die  pulsed-signal  breakdown  criterion.  The  relative  importance  of  this  effect  is  some¬ 
what  difficult  to  define,  since  it  depends  on  a  somewhat  arbitrary  parameter,  the  final  electron 

density.  The  final  electron  density  level  usually  used  in  the  pulsed  breakdown  criterion  is 

*-3  -g 

the  critical  electron  density,  which  in  unit;;  of  cm  is  approximately  equal  to  10  times  the 
square  of  die  microwave  frequency  (for  example  see  References  250  and  251).  Most  authors 
generally  state  that  when  die  initial  ionization  from  die  flow  field  is  already  equal  to  or 
greater  than  the  critical  value,  then  the  breakdown  field  is  essentially  zero,  since  the  gas 
may  already  be  considered  to  be  broken  down  and  there  can  be  no  penetration  of  the  micro- 
wave  field  through  the  plasma.  On  die  other  hand,  the  CW  breakdown  criterion  normally 
used  does  not  depend  on  the  electron  density  level,  but  only  on  the  rate  of  Increase  of 
electron  density  from  the  microwave  field.  Thus,  no  such  critical  electron  density  effect  as 
described  above  appears.  But  the  two  physical  situations  cannot  be  this  different  merely 
because  one  is  for  pulsed  signals  and  die  other  is  for  CW.  The  difficulty  is  obviously  a 
result  of  the  two  arbitrary  criteria  used  in  defining  breakdown. 


In  the  case  of  the  reentry  induced  flow  field  of  a  slender  cone  body,  it  is  possible  to  have 

greater  than  critical  electron  density  in  a  thin  layer,  such  that  total  signal  blackout  does  not 

occur.  Thus  there  will  be  a  significant  transmission  of  power  through  the  plasma  which  will 

increase  linearly  with  input  power  until  breakdown  is  approached,  when  the  transmitted  power 

will  attenuate  nonlinearly  due  to  the  thickening  of  the  plasma  layer  by  the  breakdown  process. 

This  indicates  that  the  only  physically  meaningful  way  to  calculate  breakdown  in  the  presence 

of  strong  flow-field  ionization  is  to  calculate  signal  transmission  as  a  function  of  input 

power  level,  taking  into  account  the  increased  level  of  ionization  induced  by  the  field,  as  has 
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been  suggested  by  Epstein 

A  careful  consideration  of  the  form  of  the  breakdown  criterion  is  not  only  necessary  for 
pulsed  signal  breakdown  in  the  presence  of  ionization,  but  also  for  breakdown  with  convection, 
as  discussed  in  greatev  detail  below.  A  meaningful  prediction  of  breakdown  in  the  presence 
of  flow -field  ionization  with  convection  requires  one  to  calculate  the  electron  density  distribu¬ 
tion  over  the  antenna  and  the  power  transmitted  through  this  plasma  as  functions  of  the  input 
power.  The  .Importance  of  this  consideration  lies  mainly  in  its  effect  on  the  breakdown 
criterion  and  not  so  much  in  the  phenomenon  of  alteration  of  the  antenna  near-field  and 
impedance  by  the  plasma,  although  the  latter  may  be  important  considerations  if  the  electron 
density  is  initially  very  high. 

6.2.3  HIGH  GAS  TEMPERATURE 

It  has  been  recognized  for  some  time  that  the  flow-field  gas  temperatures  produced  by  hyper¬ 
sonic  reentry  would  hate  a  significant  effoct  on  antenna  breakdown.  This  effect  is  known  to 
be  separate  from  and  additional  to  the  effects  of  ionization  and  gas  density,  which  both 
accompany  high  temperature,  but  it  is  the  least  understood  effect  theoretically. 

Direct  current  breakdown  at  high  temperature  has  been  investigated  experimentally  at  a 

relatively  early  date^^.  It  was  found  that  up  to  temperatures  of  1400°K  no  deviation  from 
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cold  gas  results  could  be  observed  .  The  investigation  of  this  phenomenon  by  the  use  of 
a  shock  tube  in  order  to  produce  higher  temperatures  was  carried  out  by  Sharbaugh  et  al. 


for  nitrogen  261 .  These  investigators  found  that  significant  departures  from  cold  gas  break¬ 
down  occur  at  temperatures  of  2000°K  and  that  at  4000°K  the  breakdown  potential  is  one-half 
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the  cold  gas  value.  Their  analysis  of  this  effect  was  based  on  a  theory  of  space-charge 
alteration  of  the  field  by  ionization  in  the  shock  tube,  although  certain  anomalies  were  present, 
which  they  ascribed  to  ionization  relaxation  in  the  shock  tube.  The  same  investigators 
carried  out  further  experiments  in  both  nitrogen  and  air,  in  which  the  electron  density  was 
also  measured263.  These  data  showed  an  apparent  correlation  between  the  reduction  in 
breakdown  voltage  and  the  electron  density. 

Although  dc  breakdown  is  affected  by  surface  phenomena  and  space-charge  effects,  it  may  be 

that  the  above  results  provide  an  indication  of  effects  of  temperature  on  basic  processes  of 

gaseous  breakdown.  In  fact,  it  Is  apparent  that  they  are  in  qualitative  agreement  with  more 

264,  265,  252  „  ,  264  .  _  .  .  .252 

recent  shock  tube  experiments  on  microwave  breakdown  .  Taylor  and  Light 

have  found  significant  effects  at  temperatures  above  3000°K  in  the  constituents  of  air.  Their 

data  are  plotted  in  Figure  6-2  In  terms  of  the  deduced  values  of  net  ionization  frequency  as 

200 

a  function  of  the  effective  field.  The  cold  air  data  of  Scharfman  ,  as  well  as  data  of 
various  otter  experimenters  reduced  by  Scharfman  to  this  form,  are  plotted  in  this  same 
figure  for  compariosn.  It  is  apparent  that  for  the  shock  tube  conditions  there  is  a  significant 
increase  in  ionization  frequency  over  the  cold  air  values,  corresponding  to  a  lowering  of  the 
breakdown  field  by  a  factor  of  the  order  of  two. 

Some  attempts267  have  been  made  to  predict  theoretically  high  temperature  air  effects,  but 
tto  results  of  these  fail  to  agree  with  the  experimental  data  of  Figure  6-2.  Both  experimenters 
apparently  have  eliminated  shock-induced  ionization  effects  in  their  experiments.  Also  the 
effects  of  attachment  cannot  be  important,  since  the  cold  air  attachment  frequency  is  orders 
of  magnitude  lower  than  the  values  of  ionization  frequency  obtained  in  the  experiments.  Thus 
the  observed  effects  are  true  gas  temperature  effects,  except  for  the  possible  air  composition 
effects  of  shock-induced  dissociation.  *  It  is  apparent  that  excitation  of  air  species  at  high 

*  In  an  attempt  to  explain  the  data,  Light  suggested  a  simple  model,  which  depends  on  the 
dissociation  of  air  molecules  to  a  significant  extent.  This  model  is  based  on  the  effect  of 
vfb rationally  excited  molecules  in  changing  the  energy-loss  rates  of  electrons  in  inelastic 
collisions. 
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Figure  6-2.  High  Gas  Temperature  Effects  on  Ionization  Frequency 
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temperature  is  an  important  factor  in  changing  the  electron  energy  distribution  (and  hence  the 
ionization  frequency)  through  the  effect  on  the  inelastic  collision  frequencies  of  electrons 
with  neutral  particles,  but  the  theoretical  prediction  of  this  effect  is  impossible  without 
solving  the  Boltzmann  equation  for  breakdown  of  high  temperature  air. 


As  an  alternative  to  the  difficult  task  of  rigorous  theoretical  prediction,  a  simple  curve  fit  of 
the  breakdown  data  has  been  obtained  in  the  form  of  a  power  law,  as  follows: 


u/p* 
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(E  /I00p*)2 
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,  -(3000/T)5 

1-e 


2.67 

4 

-6.4  x  10 
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This  function  has  been  plotted  in  Figure  6-3  in  comparison  with  Light's  high  temperature  data 
and  the  cold  air  data.  It  can  be  seen  to  be  a  good  fit,  considering  the  probable  accuracy  of 
the  data  of  Light,  who  used  theoretical  calculations  of  thermochemical  relaxation  behind  the 
shock  to  deduce  the  air  temperatures  in  his  experiments.  (Note  that  there  is  only  200  K 
difference  between  each  set  of  points. ) 
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It  should  be  pointed  out  at  this  point  that  the  cold  air  data  of  Scharfman  and  Morita  used  in 

these  figures  represents  a  departure  from  the  theory  of  MacDonald  at  low  pressure  (see 

Paragraph  6.1),  which  does  not  permit  p/p  versus  E  /p  to  collapse  uto  a  single  curve 

e  244 

independent  of  pA.  The  Scharfman  method  is  similar  in  this  regard  to  the  Brown  theory, 

which  preceded  the  MacDonald  theory.  MacDonald's  theory  would  appear  to  be  the  best  for 

cold  air,  since  it  has  been  well  tested  against  cavity  breakdown  data.  But  it  is  fallacious  to 

268  .  24 

assume  that  it  gives  the  true  values  of  ionization  frequency  in  air.  Epstein  and  Lenander 
recently  have  concluded  from  numerical  solutions  of  the  Boltzmann  equation  that  the  use  of  the 
effective  field  concept  and  the  neglect  of  nonuniform  field  effects  may  lead  to  inaccurate 
inferences  of  ionization  frequency  at  low  pressure.  The  Scharfman  data  were  obtained  with 
short  pulse  techniques,  as  were  all  of  the  high  temperature  data,  and  the  results  were  not 
dependent  on  the  calculation  of  diffusion  losses,  so  that  they  possibly  could  be  more  reliable. 
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Figure  6-3.  Curve  Fit  to  Ionization  Frequency 


Unfortunately  the  ionization  frequency  has  never  been  measured  directly  except  for  hydrogen, 
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by  Cottingham  and  Buchsbaum  ,  using  an  optical  method. 

In  conclusion,  the  format  of  Scharfman,  shown  in  the  figures  of  this  section,  is  recommended 
in  predictions  of  gas  temperature  effects.  It  has  the  advantage  of  simplicity  in  terms  of  an 
excellent  algebraic  curve  fit  w'thin  the  accuracy  of  the  data.  The  high  temperature  data  could 
be  converted  to  the  MacDonald  format  (X  =  3.21  cm  in  the  shock  tube  experiment);  but  the  data 
would  become  much  too  meager  because  of  the  additional  parameter  involved.  The  data  are 
too  meager  even  in  the  economical  format  used. 

6.2.4  CONVECTION 

The  period  of  time  in  which  a  given  volume  element  of  air  comes  under  the  influence  of  the 
field  of  an  antenna  mounted  on  a  reentry  vehicle  is  on  the  order  of  other  characteristic  times 
for  breakdown.  Hence,  the  collection  effect  can  be  important.  However,  many  alternative 

♦  i 

theoretical  analyses  of  this  effect  have  appeared,  and  the  most  important  aspect  of  the 
theoretical  prediction  problem  is  to  resolve  which  theory  applies. 

Three  theoretical  formulations  of  the  effect  of  convection  on  breakdown  are  presently  known, 

and  in  fact  four  are  possible  logically.  These  theories  derive  from  the  combination  of  two 

alternative  views  of  the  electron  flux  in  a  moving  plasma  and  two  views  of  the  effects  of  the 

electron  flux  on  breakdown.  The  contribution  of  convection  to  the  electron  flux  in  a  plasma 

can  be  assumed  to  be  equal  to  die  product  of  the  electron  density  and  the  flow  velocity  (free 

convection),  or  one  can  take  into  account  space-charge  effects  by  which  there  is  an  additional 

"force”  on  the  electrons  caused  by  the  convection  of  the  ions.  The  convective  flux  can  be 

coupled  to  the  diffusive  flux  in  the  continuity  equation,  or  the  convective  flux  can  be  treated 

as  an  electron  loss  term  whereby  the  electrons  are  present  in  the  microwave  field  for  a 

finite  time.  The  possible  combinations  of  these  assumptions  give  four  theories,  of  which 

only  the  combination  of  space  charge  enhanced  flux  in  the  finite  time  theory  has  not  been 
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reported.  Romig  and  later  I^ante  used  free  convection  coupled  to  diffusion,  Kelly  and 
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Margenau  used  free  convection  in  the  finite  time  theory,  and  Cottingham  used  space- 
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charge  enhanced  convection  coupled  to  diffusion. 


Each  oi  the  four  convection  theories  adds  a  term  to  the  breakdown  criterion  for  the  ionization 
frequency. 


(6-31) 


where  v  is  the  ionization  frequency  for  breakdown  in  the  absence  of  convection  and  V  is 
o  v 

the  added  effect  of  convection.  Table  6-1  gives  the  form  of  the  contribution  V  for  each  of 

v 

the  four  theories. 


Table  6-1.  Contribution  of  Convection  to  the  Breakdown  Criterion  for  v 


Coupled  to 

Diffusion 

Finite 

Time 

Free 

Convection 

2 

V 

~4D 

(v/L^in  (n^/n  ) 

Space  Charge 
Enhanced 

v2D 

4D2 

a 

(vD/LDjAn  (n^) 

In  this  table  v  is  the  flow  velocity,  L  is  the  length  of  the  microwave  field  region  parallel  to 
the  flow,  D&  is  the  amblpolar  diffusion  coefficient,  and  D  is  the  diffusion  coefficient  (which 
may  te  free,  ambipolar,  or  transitional,  depending  on  the  plasma  density). 


The  derivation  of  these  theories  begins  with  the  electron  continuity  equation  in  the  form 


^  •  r  +  — 

at 


V 

n 


(6-32) 


where  the  electron  flux  T  is  given  by  either  of  two  alternative  formulas.  In  the  free 
convection  theory  it  is  assumed  that 
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F  =  nv  -  V  (Dn) 


(6-33) 


Thus,  neglecting  V  •  (v/D) 

72  (Dn)  -  v/D  •  7  (Dn)  +  v  =  77  (6-34) 

not 

In  the  finite  time  theory,  the  coupling  between  diffusion  and  convection  is  neglected.  Then 
the  eigenvalues  of  the  two  differential  operators  can  be  assigned  independently. 

7 2  (Dn)  =  -Dn//?  (6-35) 

v  •  7  (Dn)/Dn  =  a  (6-36) 

where  a  is  a  separation  constant.  Integration  over  a  region  of  length  L  gives 


n  =  noexp  (a  L/v)  (6-37) 

Defining  breakdown  in  terms  of  electron  density  ^  and  substitution  back  into  Equation  6-34 
gives  the  breakdown  criterion. 

v  *  D/A2  +  v/L  in  ~  +  l/n  ~  (6-38) 

not  ' 

o 


The  other  alternative  in  the  free  convection  theory  is  to  couple  the  diffusion  and  convection 
together.  A  change  of  variables  then  transforms  Equation  6-34  to  canonical  form.  Let 


(6-39) 
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where  v  is  distance  along  the  streamline.  The  equation  for  0  is,  again  ignoring  (7  •  v/2D), 


2  %  _ 

7  0+  i//D  -  l/Dn  ^  -  (v/2D)  0  =  o 


(6-40) 


Now  if  the  eignevalue  of  v  iu  -1/A  ,  then  the  breakdown  criterion  is 


v  =  D/A2  +  \2/db  +  i/n  Is 

o  t 


(6-41) 


The  other  theory  of  convective  flux  recognizes  the  possibility  of  space-charge  coupling  by 
treating  ions  (+)  and  electrons  (-)  separately. 


F_  =  n_v  -v  (D_n_)  -n  En 


(6-42) 


T+  =  n+v-v  <D+n+)+p+En+ 


(5-43) 


where  n  is  the  mobility  and  E  is  the  space-charge  field.  In  the  steady  state  the  electron  and 
ion  fluxes  are  equal,  giving 


r  =  r  = 


K  +  **+>  nJ*+v  -  n+v  (D_n_)  -  p_n_7  <D+n+) 


H  n  +  u  n 

-  -  +  4- 


(6-44) 


Two  limiting  forms  of  this  result  may  be  applied.  At  low  electron  densities  D  n  =  D  n 

,  .  +  +  - 
giving 


K+"+  1 

=  \  7p-^p_  J  D-d-v  ~  *  (D.n.) 


=  ( D_/D  )  n  v  -  7  (D  n  ) 


(6-45) 
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Also,  at  high  electron  density  n+  =  n  ,  giving 

r  *  n  v  -( — 77 - J  V  p  n  )- ( — — - V  (Dj»J 

-  \  f  -  -  \  M_  +  M+  /  +  + 


(6-46) 


Under  the  reasonable  assumption  that  the  ratio  of  mobilities  is  independent  of  position,  this 
gives 


r  »  n  v  -  V  P^n  )  (6-47) 

Now  since  D  =  D_  at  low  electron  density  and  D  *  D  at  high  electron  density,  we  may  infer 
*  a. 

the  following  approximation  during  transitional  ^.Jfusion  p  -  Dg). 

r  =  (v/D  )d  n  -  V  P  n)  (6-48) 

A'  8  S 

The  velocity  is  thus  enhanced  by  the  ratio  D  /D  ,  giving  the  results  shown  in  the  second 

S  ft 

row  of  Table  6-1. 

It  is  presently  possible  to  resolve  the  discrepancies  presented  by  the  above  four  theories, 
to  the  extent  that  they  can  all  be  shown  to  be  special  cases  of  the  general  situation.  First  we 
point  out  that  whan  the  initial  electron  density  is  sufficiently  high,  the  two  rows  of  Table  6-1 
automatically  become  identical,  because  D  =  D  .  Thus  there  is  no  need  to  resolve  the 

A 

discrepancy  between  the  free  convection  and  space-charge  enhanced  convection  theories  except 
at  low  initial  electron  density. 

The  resolution  of  the  disparity  between  columns  of  the  table  has  been  established  by  W.  C. 
Taylor  of  Stanford  Research  Institute  (private  communication),  who  has  pointed  out  that 
the  finite-time  theory  is  actually  a  special  case  of  the  diffusion-coupled  theory.  This  is 
demonstrated  by  solving  Equation  6-40  for  the  electron  density  at  the  downstream  edge  of  an 
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aperture  antenna  in  a  gas  flow  of  uniform  velocity.  We  depart  from  Taylor's  analysis  and 
Fante*  s  solution  to  pose  the  following  one-dimensional  problem  : 


-(v/D  )~  +vft>  n  =  0 
,  2  V  v’  dx  s 

dx 


(6-49) 


where  D  is  either  D  or  D  ,  depending  on  whether  the  free-convection  theory  or  the  space- 
v  s  a 

charge  enhanced  convection  theory  is  under  discussion.  Outside  the  region  o  <  x  <  L,  the 
ionization  frequency  vanishes  and  the  value  of  is  taken  much  smaller  than  its  value  inside 
this  region.  The  values  of  n  and  its  derivative  are  matched  at  the  boundaries.  The  solution 
for  the  electron  density  at  x  ^  L  is 
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A  =  v/v. 
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*  This  particular  problem  is  posed  to  obviate  the  following  difficulties: 

1)  Fante  solved  a  problem  in  two  dimensions  with  uniform  flow  and  with  a  constant  background 
ionization  rate,  as  an  inhomogeneous  differential  equation;  but  unfortunately  the  particular 
solution  becomes  infinite  when  the  classical  breakdown  criterion  for  parallel  plate  breakdown 
in  the  transverse  direction  is  satisfied  (y/D  =  7r2  A^2),  a  fact  apparently  unnoticed  by  Fante;  and 

2)  Taylor  included  diffusion  perpendicular  to  the  flow,  which  gives  results  which  do  not  repro¬ 
duce  a  constant  background  electron  density  when  v  vanishes.  No  sacrifice  of  generality  is 
made  in  terms  of  the  relative  importance  of  convection  when  the  problem  is  posed  in  simplified 
form. 
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where  is  the  breakdown  value  from  the  simple  diffusion  coupled  theory  and  nQ  is  the 
electron  density  at  negative  x. 


(*/L)2  +  (v/2Dy)2 


(6-51) 


Breakdown  defined  as  infinite  n/n  is  determined  by  the  solution  of  the  transcendental 

o 

equation 


(bV^i  ),  x/2  <  B  ylTi 


<  7T 


(6-52) 


In  addition  when  V  is  much  less  than  V^t  the  above  equations  have  the  following  simple 
approximate  form. 

n/nQ  =  exp  fpIX^/vD^  ,  v « (6-53) 


If  breakdown  in  defined  in  terms  of  a  finite  electron  density  n^,  at  which  the  exponential 
approximation  (6-50)  still  applies,  then  the  breakdown  criterion  is 

V  =  ^vDg/LDvjin  (t^/n^,  V  « (6-54) 

The  exact  solution,  Equation  6-50,  and  the  exponential  approximation,  Equation  6-53,  are 

shown  for  certain  specific  cases  in  Figures  6-4  through  6-6,  for  a  velocity  of  20  kft/sec, 

6  2  “1 

assuming  D  jp  =1.6  x  10  cm  /Torr/sec  and  D  /D  =100.  Three  values  of  pL  are 

assumed,  1,  10,  and  100  cm/Torr.  The  first  case  considered  is  the  free-convection  theory 

assuming  free  diffusion  applies.  Figure  6-4  shows  that  the  simple  exponential  approximation 

does  not  hold,  and  in  addition *the  actual  breakdown  occurs  at  much  lower  values  than  v,m 

b 

Thus  in  the  free  convection  theory  the  breakdown  ionisation  frequency  is  given  approximately 

by  the  smaller  of  the  two  criteria  given  by  Equations  6-52  and  6-54,  with  D  =  D  . 

s  v 
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Figure  6-5  shows  the  results  assuming  amblpolar  diffusion  applies  (for  either  theory  of  con¬ 
vective  flux).  This  time  the  exponential  approximation  predicts  breakdown  very  precisely  at 
high  pressure,  and  the  simple  criterion  v  =  begins  to  apply  at  low  pressure.  Therefore 
the  correct  breakdown  criterion  is  less  than  either  V.  or  Equation  6-54  with  D  =  D  ,  and  a 

D  S  V 

fair  approximation  is  to  take  the  smaller  of  the  two. 

Figure  6-6  is  for  the  Cottingharo  theory,  in  which  space-charge  enhanced  convection  is 

assumed  and  transitional  diffusion  is  considered.  These  curves  were  calculated  applying 

Equation  6-23  to  n/n  ,  assuming  T  /  (n  A2)  =  2  x  104  cm/°K.  Also,  v  in  this  figure  has 
o  e  o  o 

the  meaning  of  the  ambipolar  value,  the  same  as  in  the  previous  figure.  Two  different  forms 
of  exponential  approximation  are  shown,  that  for  the  ordinary  L/v  theory,  and  that  for  the 
Cottingham  velocity  with  free  diffusion.  The  transition  from  free  to  ambipolar  diffusion 
shows  up  at  n  =  10^n  .  All  three  curves  are  asymptotic  to  v  =  l^.  The  interpretation  of  these 
results  is  confusing  in  light  of  the  multi-valued  electron  density  as  a  function  of  the  ionization 
frequency.  But  if  the  criterion  of  a  finite  electron  density  is  used,  then  the  breakdown 
ionization  frequency  is  single  valued.  The  initial  condition  corresponds  to  having  approximately 
one  electron  in  a  volume  of  gas  equal  to  tiie  extent  of  the  field  region,  so  that  breakdown 
probably  corresponds  to  the  portion  of  the  curves  at  which  the  transition  is  complete.  Thus 
breakdown  is  determined  approximately  by  the  smaller  of  or  Equation  6-54  with  D^=Dv=D^. 

The  conclusions  from  Figures  6-5  and  6-6  are  essentially  equivalent,  and  it  is  possible  to 

show  that  the  results  of  Figures  6-4  and  6-5  are  consistent  if  transitional  diffusion  is 

assumed  and  if  the  difference  in  values  of  V,  is  accounted  for.  Hence  ambipolar  diffusion 

b 

always  accompanies  breakdown  with  convection,  since  the  curves  of  the  Figure  6-4  never 
reach  large  electron  densities  before  v  =  it.  Also  the  contribution  of  convection  to  the  break¬ 
down  criterion  is  always  smaller  than  either  of  the  following,  regardless  of  the  initial 
electron  density. 

Vy  ^v/L)  ,  Vv  *  v2/4Da  (6"55) 
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6.2.5  GRADIENTS  OF  GAS  PROPERTIES 

As  previously  noted,  the  boundary  layer  surrounding  a  slender  cone  vehicle  in  hypersonic 
flight  has  strong  gradients  of  density,  temperature,  velocity,  and  electron  density.  In  addi¬ 
tion,  the  entire  flow  field  is  highly  nonunifcrm,  with  large  changes  in  properties  in  passing 
from  one  region  to  another.  Thus  the  importance  of  these  nonuniformities  must  be  evaluated 
before  the  breakdown  analysis  can  be  simplified. 

271 

Fante  has  analyzed  the  situation  of  nonuniform  media  in  terms  of  the  effect  of  an  ambi- 
polar-dif  fusion  layer  and  an  outer  free -diffusion  layer  covering  an  aperture  antenna.  The 
results  of  this  calculation  show  that  even  a  small  ambipolar-diffusion  layer  can  drastically 
reduce  the  breakdown  field  strength.  In  general,  it  can  be  seen  from  the  effects  of  reentry 
that  several  factors  contribute  to  spatial  nonuniformity  of  different  terms  in  the  continuity 
equation.  The  flow  velocity  is  obviously  nonuniform  in  space  as  determined  by  the  flow- 
field  properties.  The  diffusion  coefficient  will  have  gradients  due  to  molecular  density 
gradients  and  ionization  gradients  (transition  to  ambipolar  diffusion  in  high  electron-density 
regions).  If  the  effective  field  is  included  in  the  formula  for  the  diffusion  coefficient, 
gradients  will  also  be  due  to  field  nonuniformity.  The  frequencies  of  ionization  and  attach¬ 
ment  have  spatial  gradients  contributed  by  variations  in  particle  number  density,  field 
intensity,  and  also  air  temperature.  Profiles  of  chemical  species  which  are  distributed  in 
the  flow  field  would  conceivably  contribute  to  gradients  in  both  the  diffusion  coefficient 
(molecular  weight  effect)  and  the  ionization  and  attachment  frequencies.  Finally,  the  ioniza¬ 
tion  gradients  would  affect  pulse  breakdown  in  terms  of  the  spatial  variation  of  initial  electron 
density,  although  this  effect  is  probably  not  very  important. 

These  effects  can  be  accounted  for  in  principle  by  solving  the  continuity  equation  as  a  differen¬ 
tial  equation  with  nonconstant  coefficient,  just  as  in  the  problem  on  nonuniform  antenna  fields. 
A  fair  approximation  is  to  treat  the  breakdown  on  a  slender  cone  vehicle  separately  in  terms 
of  three  different  regions,  boundary  layer,  invlscid  layer,  and  free  stream.  The  contribution 
of  the  first  two  regions  to  the  effect  should  be  to  give  an  effective  diffusion  length  of  the  order 
of  the  thickness  of  the  particular  region  in  question. 
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The  presence  of  high-electron  density  In  a  thin  layer  Is  not  necessarily  a  sufficient  condition 
for  breakdown,  but  rather  the  adjacent  flow-field  region  must  also  ionize  to  cause  the  plasma 
to  thicken,  leading  to  breakdown.  This  aspect  of  the  problem  is  related  to  the  general  question 
of  breakdown  criteria  discussed  above. 

6.2.6  EFFECTS  OF  GAS  CHEMICAL  COMPOSITION 

Two  sources  of  chemical  composition  effects  in  the  flow  fields  of  slender  reentry  vehicles  are 
dissociation  of  the  air  and  ablation  of  gases  from  the  heat  protection  shield.  These  effects 
can  be  accounted  for  in  principle  by  reviewing  the  basic  data  on  each  of  the  species  involved. 
However,  much  of  the  breakdown  data  have  been  obtained  in  terms  of  materials  as  candidates 
for  delib*  irate  injection  into  the  flow  as  alleviants.  Therefore  these  data  are  reviewed  in 
Section  8. 

6.3  ILLUSTRATIVE  CALCULATIONS 

The  breakdown  power  of  an  antenna  in  the  reentry  environment  can  be  estimated  from  the 
results  developed  above.  In  fact,  it  is  possible  to  write  these  results  in  terms  of  a  closed 
form  algebraic  equation,  utilizing  the  following  equation  for  the  power  output  of  an  antenna  273. 

P  =  1.33xl0"3E2A  (6.56) 

where  P  is  in  watts,  E  is  the  maximum  rms  field  intensity  in  volts/cm,  and  A  is  the  effective 
2 

aperture  in  cm  ,  defined  in  terms  of  the  near  field  and  the  input  power. 


Equations  for  the  effects  of  gas  density,  ionization,  temperature,  and  convection  are  discussed 
above  and  can  be  used  here.  The  most  uncertain  of  these  equations  is  the  high  gas  temperature  * 
Equation  6-30,  which  is  a  curve  fit  to  some  very  meager  data  and  which  is  based  on  a  simplified 
view  of  cold  air  breakdown.  Before  proceeding  »o  the  application  of  this  equation  to  illustrative 
calculations,  let  us  indicate  the  extent  of  the  discrepancy  between  this  simplified  theory  and 
that  of  MacDonald  .  Figure  6-7  shows  this  comparison  for  three  values  of  X/A .  The 
limitations  and  difficulties  of  the  data  and  theories  are  discussed  more  fully  above,  but 
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1. 


Figure  6-7  can  be  used,  if  desired,  to  correct  the  values  of  breakdown  power  calculated 
below  to  the  values  corresponding  to  the  MacDonald  theory.  Of  course  the  method  used  here 
is  conservative,  since  it  will  always  predict  lower  br<  ak>  own  power  than  MacDonald's  theory. 


The  resulting  equation  for  CW  breakdown  power  can  be  written  in  proper  variable  form  by 
normalizing  P. 


P/Pq  =  (A/A  ) 


,  -(3000/T)' 
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=  (A/A  )  f  ( T) 
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=  1. 12  x  1013/A2 
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In  these  equations  P  is  in  watts,  T  and  T  are  in  °K,  p*  and  p  are  in  Torr,  and  all  other 
o  e 

quantities  are  in  cgs  units.  This  normalization  removes  the  temperature  effect  as  a  simple 

mult4>licative  factor  f  (T)  and  expresses  the  remaining  reentry  effects  in  terms  of  proper 

variables.  For  fixed  values  of  A/A  and  A/L,  P  is  a  function  of  the  proper  variables  p*  A, 

2  o 

nQ  A  ,  and  v,  if  the  electron  temperature  is  held  constant. 

P  is  plotted  parametrically  as  a  function  of  these  proper  variables  in  Figures  6-8  through 

6  4 

6-10  using  the  following  constants:  D  p*  =  1. 6x  10  ,  T  =  2  x  10  ,  D  /D  =  100. 

-  e  -  a 
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SECTION  7 

NONUNEAR  INTERACTIONS 


Electromagnetic  fields  of  sufficient  intensity  can  induce  a  change  in  the  state  of  a  partially 
ionized  plasma  due  to  field  energy  coupling  into  the  electron  gas.  As  a  result,  the  electron 
distribution  function  is  altered  and  becomes  dependent  on  the  electric  field  intensity. 
Consequently,  transport  coefficients  such  as  thermal  and  electrical  conductivities,  which 
are  deduced  from  the  distribution  function,  become  dependent  on  the  local  field  intensities. 
While  various  forms  of  nonlinearities  may  be  exhibited  from  the  interaction  of  electromagnetic 
radiation  with  a  plasma,  the  particular  class  covered  in  this  section  is  restricted  to  the 
case  where  the  incident  field  deposits  energy  in  .he  plasma  and  modifies  the  electron  tempera¬ 
ture,  electron  density,  and  effective  collision  f'equency.  The  present  interest  is  in  the 
propagation  of  electromagnetic  radiation  in  a  plasma,  and  the  conductivity  coefficient  is 
of  paramount  importance. 

7.1  SURVEY 
274 

Papa  has  performed  a  comprehensive  survey  of  nonlinear  electromagnetic  wave 
propagation  in  plasmas.  His  conclusions  indicate  that  there  are  many  kinds  of  nonlinear 
interactions  possible,  resulting  from  the  fact  that  the  physical  mechanisms  involved 
have  such  a  large  number  of  characteristic  times,  especially  in  terms  of  electronic 
interactions.  The  various  possible  combinations  of  relative  magnitudes  of  these  charac¬ 
teristic  times  lead  to  the  definition  of  myriad  possible  situations.  For  example,  if  the 
electron-electron  collision  time  is  much  less  than  the  times  for  the  relaxation  of  electron 
density  and  temperature  and  the  diffusion  time,  then  the  isotropic  part  of  the  distribution 
function  is  Maxwellian.  But  If  the  electron  temperature  relaxation  time  is  much  shorter 

than  all  other  characteristic  times  except  the  period  of  the  electromagnetic  field,  then 
275 

the  Margenau  distribution  applies. 

When  the  power  level  of  a  signal  in  a  plasma  is  increased,  the  first  nonlinear  effects  to 
be  observed  are  collision  frequency  effects,  especially  in  plasma  formed  in  air.  The 
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dependence  of  collision  frequency  on  electron  velocity  changes  the  effective  collision 
frequency  as  the  electron  energy  gained  from  the  field  is  increased.  Strictly  speaking, 
the  plasma-electromagnetic  wave  interaction  is  always  nonlinear,  since  the  electron 
energy  depends  on  the  field  strength  under  all  conditions.  In  addition,  the  plasma 
constitutive  parameters  depend  on  the  signal  frequency,  and  they  are  not  additive  in 
the  case  of  mixtures  of  different  species. 

Various  approaches  to  the  definition  of  effective  parameters  for  a  plasma  have  been  used. 

The  problem  is  to  define  two  parameters,  which  are  related  to  the  complex  conductivity 
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of  the  plasma.  Molmund  used  a  complex  collision  frequency  to  obviate  the  difficulty 

277 

of  non-constant  collision  frequency,  while  Shkarofsky  used  auxiliary  real  functions. 

278 

The  method  of  Whitmer  and  Herrmann*  ,  who  defined  effective  values  of  collision 
frequency  and  plasma  frequency  for  direct  use  in  the  conductivity  formula,  is  the  most 
convenient  to  use.  The  latter  authors  give  useful  numerical  results  for  actual  gases  in 
graphical  form. 
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Bakshi  et  al.  have  shown  that,  all  the  above  approaches  are  equivalent,  and  they  have 
derived  a  general  procedure  for  calculating  effective  parameters.  The  effective  small 
signal  conductivity,  defined  as  the  ratio  of  the  current  density  to  the  electric  field,  in 
the  limit  of  vanishing  field,  is  defined  in  terms  of  response  functions. 


a(co)  =  c0  2  [c  (u> )  +  is  (u>)] 

(7-1) 

00 

C  (u>)  =  f  R  (y'  cos  uj  dy 

(7-2) 

OB 

S  (u>)  *  f  R  (y)  sin  wy  dy 

(7-3) 

where  R(y)  is  the  average  response  of  an  electron  in  the  plasma  to  a  unit  amplitude 
step  function  ox  electric  Held  starting  at  y  ■  0.  These  authors  also  give  graphical 
results  for  idealized  forms  of  dependence  of  collision  frequency  on  electron  speed. 
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Haskell  et  al.  have  considered  the  above  theory  of  effective  parameters  for 

strong  fields,  i.  e.  non-Maxwellian  velocity  distributions.  As  in  the  case  of  vanishing 

fields  considered  above,  each  plasma  parameter  undergoes  a  transition  from  a  low 

frequency  value  to  a  high  frequency  value  in  the  neighborhood  of  the  average  collision 

frequency.  However,  in  this  case  the  results  depend  on  the  field  strength  at  low 

frequency.  But  even  though  the  low-to-high  frequency  transition  occurs  when  the 

collision  frequency  is  independent  of  velocity,  it  is  independent  of  field  strength 

in  this  case.  Thus,  for  a  constant  electron  density,  all  nonlinear  effects  must  come 

from  non-constant  collision  frequency  effects.  For  collision  frequency  as  an  increasing 

function  of  velocity,  both  the  effective  plasma  frequency  and  the  effective  collision 

frequency  increase  with  increasing  field  strength  at  low  frequency  (or  low  altitude). 

2 

Unfortunately,  the  ratio  /v c>  which  is  important  in  thin  plasma  for  slender  reentry 

vehicles  at  low  altitudes,  tends  to  increase  with  field  strength. 

The  nonlinear  collision  frequency  effect  is  absent  at  high  altitude,  so  that  nonlinear 

production  of  electrons,  or  breakdown,  is  the  only  important  phenomenon.  Nonlinear 

breakdown  can  take  a  variety  of  forms.  If  very  little  ambient  plasma  is  present, 

very  little  effect  will  be  observed  until  the  breakdown  field  is  reached.  The  ionization 

frequency  increases  too  rapidly  with  field  strength  for  appreciable  prebreakdown 

nonlinear  effects  to  occur.  Subsequent  increase  in  power  level  beyond  breakdown 

causes  the  discharge  plasma  to  thicken,  so  that  output  power  is  generally  limited  to 

a  low  leve*.  When  the  initial  electron  density  is  high  (near  critical)  the  process  of 

breakdown  has  the  apparent  behavior  of  a  nonlinear  effect,  siuce  the  transmitted 

power  will  begin  to  level  off  as  input  power  is  increased  and  the  plasma  layer  thickens. 
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These  effects  have  been  observed  experimentally  by  Chown  et  al.  ' 

7.2  BROADBAND  SIGNAL  EFFECTS 

The  most  important  aspect  of  nonlinear  phenomena  is  their  effects  on  broadband  signals. 
In  terms  of  prebreakdown  nonlinear  collision  frequency  effects,  it  appears  that  the  most 
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important  effects  will  occur  in  a  band  of  frequencies  near  the  average  collision  frequency. 
But  the  plasma  dispersion  tends  to  decrease  as  the  field  increases  ,  at  field  strengths 
below  breakdown,  from  the  results  of  Haskell.  Thus  the  only  important  broadband 
signal  effect  is  breakdown.  The  problem  of  multi-frequency  propagation  has  been 

OQO 

considered  by  Mott  in  conjunction  with  antenna  breakdown  in  air.  The  breakdown 

depends  not  only  on  frequency  but  also  on  the  envelope  of  the  electric  field.  The  envelope 

may  be  determined  by  phasor  addition  of  the  individual  components.  A  minimum  level 

is  obtained  for  n  equal  amplitude  signals  by  in-phase  addition.  The  resulting  envelope 

peaks  determi.  e  breakdown  by  acting  as  a  continuous  wave  breakdown  field  n  times 

as  large  as  the  single  field.  The  resulting  breakdown  power  is  l/n  times  that  for 

one  frequency.  Additional  work  on  multi-frequency  signals  has  been  carried  out  by 
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Fante  and  Mullin  .  Here  an  attempt  is  made  to  determine  the  effective  electric 
field  strength  for  a  composite  signal  such  as  noise.  For  a  stationary  random  signal 
E  (t),  the  effective  field  is  given  by 

.y-Wl 

O  ui  +  v 

c 

where  v  is  the  velocity-independent  collision  frequency,  and  S  (f)  is  the  spectral 
c 

density.  The  transmission  of  power  by  the  noise  signal  i6  compared  to  that  of  a  single 
frequency  component.  At  high  collision  frequencies  there  is  no  difference  between  the 
grata.  However,  when  the  collision  frequency  is  less  than  the  average  signal  frequency, 
the  single  frequency  can  transmit  more  power  than  the  noise  signal. 
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SECTION  8 


ALLEVIATION  TECHNIQUES 

The  alleviation  of  reentry  plasma  effects  on  signal  transmission  and  of  reentry  flow  field 
effects  on  antenna  breakdown  is  generally  a  matter  of  considerable  concern  where  these 
effects  are  apt  to  be  serious.  The  solution  of  the  problem  of  alleviation  is  dependent  on 
gaining  an  understanding  of  the  phenomena  involved  in  terms  of  the  causative  factors  and 
the  mechanisms  of  interaction.  Plasma  effects  on  signal  transmission  depend  on  the  flow- 
field  ionization  and  collision  frequency  as  sole  causative  factors;  whereas  the  effects  of 
reentry  on  breakdown  are  caused  by  a  variety  of  phenomena.  Both  types  of  effect  depend 
on  basically  the  same  mechanism  of  interaction  of  the  antenna  field  with  the  medium.  The 
discussion  of  alleviation  is  best  organized  according  to  the  discussion  of  alleviation  of  the 
causative  environment  or  the  interaction,  since  overlap  may  occur  between  attenuation  and 
breakdown.  The  environment  control  is  discussed  under  the  heading  of  plasma  alleviation 
and  the  interaction  is  treated  under  electromagnetic  techniques. 

8.1  PLASMA  ALLEVIATION 

It  is  desirable  from  the  standpoint  of  both  plasma  attenuation  and  breakdown  to  minimize  the 
flow  field  ionization  in  the  vicinity  of  the  antenna.  Aerodynamic  shaping  and  material  in¬ 
jection  are  considered  in  this  vein.  In  terms  of  body  shape  effects,  the  nose  bluntness  of  a 
slender  cone  vehicle  is  the  most  important  contributor  to  ionization.  Thus,  the  use  of  a 

sharp  nose  tends  to  minimize  shock-induced  plasma.  Various  other  aspects  of  body  shaping 

285 

arc  also  considered  in  an  Aerospace  Corporation  report,  including  the  communications 
fin  and  flow  diverters.  But  these  techniques  have  no  application  to  slender  cone  vehicles, 
especially  when  the  nose  radius  can  be  made  sufficiently  small.  Various  ideas  for  liquid 
or  gas  coolants  have  been  suggested.  Some  of  these  techniques  may  possibly  alleviate  the 
high  gas  temperature  effect  on  breakdown,  as  well  as  plasma  effects  on  both  attenuation 
and  breakdown. 
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Chemical  additives  may  take  the  form  of  either  ablation  products  from  the  heat  shield  or 
active  mass  injection. 

286 

Aisenberg  has  conducted  analytical  and  experimental  studies  on  alleviation  of  plasma 
sheath  and  antenna  breakdown  effects.  They  have  been  evaluating  materials  for  the  ability 
to  reduce  electron  concentration  and  to  increase  breakdown  power  level  In  thermal  plasmas. 
The  primary  alleviation  mechanism  studied  is  attachment  of  electrons  by  electronegative 
gases.  The  plasma  simulation  facility  used  in  these  studies  is  an  argon  arc  jet.  Conditions 
simulated  include  atmospheric  pressure  down  to  13  Torr  and  temperatures  of  3000  to  5000°K. 
Diagnostics  are  conducted  by  means  of  K-band  transmission  through  the  test  section. 
Additives  tested  include  SF^,  BF^,  N^,  Freon  116,  air,  and  argon. 

Aisenberg  concluded  that  materials  which  are  good  attachers  in  glow  discharges  often  do  not 
work  in  thermal  plasmas  because  of  dissociation.  A  figure  of  merit  for  attachers  would 
appear  to  depend  on  the  molucular  weight,  the  electron  attachment  cross  section  (including 
its  energy  dependence),  and  the  resistance  to  thermal  dissociation.  Highly  electronegative 
materials  have  higher  attachment  cross  sections  in  the  molecular  form  than  in  the  atomic 
form.  Survival  lifetimes  of  molecular  additives  in  thermal  plasmas  are  thus  important  in 
determining  their  effectiveness  as  plasma  alleviative  substances. 

The  ideal  additive  material  for  alleviation  of  both  ionization  and  breakdown  is  pictured  as 
one  made  up  of  "giant  molecules,  "  which  would  actually  be  extremely  small  solid  particles 
or  liquid  drops.  These  particles  would  have  maximum  survival  times  as  well  as  providing 
surface  recombination  as  a  mechanism  for  electron  removal.  Refractory  solids  would  have 
to  bo  avoided  because  of  thermionic  emission  and  Saha- Langmuir  ionization. 

MIT  Aerophysics  Laboratory  has  been  conducting  experiments  in  a  wind  tunnel  under  the 
direction  of  E.  E.  Covert  .  In  this  experiment  a  supersonic  arc  jet  is  blown  over  the 
nose  of  a  cone  model  in  a  Mach  4  wind  tunnel.  Breakdown  of  an  X-band  slot  on  the  cone  was 
de.ermlned  under  conditions  of  injection  of  SFg.  They  have  observed  that  the  electron 
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density  was  reduced  by  an  order  of  magnitude,  but  that  the  breakdown  was  unaffected  unless 

the  additive  was  injected  immediately  in  front  of  the  anteana.  Good  results  were  also 

obtained  with  Teflon  products.  Thompson  has  measured  the  effects  of  CO,  CF^,  and 

SF  on  electron  density  and  antenna  breakdown  in  an  arc  channel.  The  first  two  of  these 
6 

materials  are  products  of  aDlation  for  carbon  and  Teflon  heat  shields,  respectively;  and 

the  last  is  a  popular  electrophilic  quenchant  candidate.  Twenty  percent  by  volume  of  each 

of  these  materials  was  added  to  the  flow,  and  experiments  were  conducted  at  temperatures 

between  3000  and  5000°K.  Carbon  monoxide  was  found  to  have  no  effect  on  either  electron 

density  or  breakdown.  CF,  reduced  the  electron  density  by  a  factor  of  2  to  3  and  increased 

4 

the  breakdown  power  by  a  factor  of  3.  SF  reduced  the  electron  density  by  more  than  a 
factor  of  3  and  increased  the  breakdown  power  by  a  factor  of  2.  The  authors  attribute  the 
breakdown  alleviation  by  the  fluorides  to  the  additional  electron  attachment  provided  by 
fluorine  and  fluoride  radicals. 

An  additional  aspect  of  plasma  alleviation  is  the  problem  of  collision  frequency  control, 
which  can  have  a  payoff  in  both  attenuation  and  breakdown.  If  the  plasma  is  thin  compared 
to  the  wavelength  of  the  signal,  as  it  tends  to  be  for  slender  cones,  then  it  is  beneficial  to 
make  the  collision  frequency  a  maximum.  Thus  one  should  locate  the  antenna  on  a  forward 
facing  part  of  the  body,  preferably  on  the  cone  frustum.  (Optimum  body  station  on  the  cone 

I 

from  the  standpoint  of  plasma  density  is  a  complicated  question  involving  nose  bluntness 
and  ablation  effects  in  a  tradeoff,  which  is  beyond  the  scope  of  this  study. )  This  approach 
tends  also  to  minimize  the  plasma  thickness  over  the  antenna.  The  effectiveness  of  this 
technique  in  alleviating  breakdown  is  dependent  on  the  desired  altitude  of  operation  of  the 
system.  It  is  a  decided  disadvantage  in  a  high  altitude  ECM  system,  since  the  increase  in 
collision  frequency  tends  to  push  breakdown  to  higher  altitudes  but  it  is  probably  helpful  in 
increasing  the  overall  power  handling  capability  if  the  breakdown  minimum  can  be  moved  to 
an  altitude  where  ionization  and  high  temperature  effects  are  weak. 
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8.2  ELECTROMAGNETIC  TECHNIQJES 

The  most  basic  technique  for  control  of  the  electromagnetic  interaction  is  the  antenna 
design.  The  use  of  a  non-charring  dielectric  window  or  covering  the  entire  antenna  area 
with  ncn-charring  low  loss  ablation  material,  such  as  Teflon,  reduces  the  near  field 
intensity  in  the  flow  field.  This  can  alleviate  linear  plasma  effects  by  reducing  the  contri¬ 
bution  of  the  plasma  to  the  antenna  input  impedance.  Breakdown  is  alleviated  by  reducing 

the  strength  of  the  near  field  and  by  increasing  the  rate  of  diffusion  loss,  as  noted  by 

288  289  281 

Scharfman  and  Morita  *  .  Chown  et  al.  point  out  that  the  alleviation  of  break¬ 

down  by  antennas  designed  to  maximize  the  breakdown  field  strength  does  not  necessarily 
increase  the  power  handling  capability,  especially  if  the  antenna  Q  is  increased  in  order  to 
do  so.  High-Q  antennas  are  characterized  by  the  fact  that  there  is  a  large  amount  of  stored 
energy  near  the  antenna  in  comparison  with  the  energy  radiated  per  cycle.  This  means 
that  the  field  will  be  quite  high  in  the  vicinity  of  the  antenna.  The  high-Q  antenna  is  unde¬ 
sirable  for  several  reasons— one  is  that  the  high  fields  imply  low  power-handling  capability, 
another  is  tha*  the  impedance  of  a  high-Q  antenna  will  be  more  sensitive  to  the  presence  of 
the  plasma  environment.  Reactive  fields,  upon  immersion  in  a  plasma,  become  dissipative, 
providing  another  means  of  absorbing  the  RF  power.  Also  themultipactor  discharge  pheno¬ 
menon  which  is  possibly  detrimental  to  the  vehicle  or  antenna  in  terms  of  heating  if  not  powe 
transmission,  is  enhanced  in  high-Q  antennas.  Multipacting  occurs  when  the  electron  mean 
free  path  exceeds  the  limits  of  validity  of  the  diffusion  theory  of  breakdown  (see  Paragraph 
6. 1).  It  tends  to  flatten  out  the  curve  of  breakdown  field  versus  pressure  in  such  a  way  that 
the  multipactor  breakdown  field  remains  approximately  equal  to  the  normal  breakdown  field 
at  the  mean  free  path  limit. 

288  289 

Scharfman  and  Morita  *  also  have  measured  the  effect  of  a  dc  bias  voltage  on  antenna 
breakdown,  to  determine  if  it  would  perform  the  function  of  sweeping  the  electrons  out  of  the 
high  field  region.  However,  this  technique  is  much  less  effective  for  well  designed  antennas 
which  have  low  near-field  intensity  than  for  those  which  already  have  a  low  power-handling 
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capability.  The  technique  does  not  work  wh  in  the  antenna  Is  covered  by  dielectric  material 
since  the  electrons  attracted  to  the  dielectric  surface  cancel  the  biasing  field  In  the  plasma. 

Static  magnetic  fields  can,  in  principle,  be  used  to  alleviate  plasma  attenuation.  Attenuation 

for  right-hand  circularly  polarized  waves  propagating  normal  to  a  plasma  layer  will  be 

reduced  in  the  presence  of  a  magnetic  field  in  the  direction  of  propagation.  This  is  because 

the  Lorentz  force  checks  the  ability  of  the  electron  generated  field  to  cancel  the  imposed 
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field.  Rothman  and  Morita  consider  the  problem  analytically  fcr  rectangular  coordinates. 
The  conditions  for  enhancement  are  that  the  electron  cyclotron  frequency  be  greater  than 
the  signal  frequency,  and  that  the  square  of  the  collision  frequency  be  much  less  than  the 
square  of  the  difference  between  the  cyclotron  frequency  and  the  imposed  angular  frequency. 
Rothman  and  Morita  carried  out  experimental  laboratory  work  using  an  RF  discharge  whlcn 
demonstrated  qualitative  enhancement  and  polarization  in  the  presence  of  a  magnetic  field. 

_  „  291,  232 

oamaddar  considers  magnetic  window  effects  in  special  cases  of  nonrectangular 

geometries.  In  particular,  he  considers  cylindrical  and  conical  geometries  which  approximate 
vehicle  shapes.  His  method  utilized  calculating  the  dielectric  tensor  as  a  function  of  the 
applied  magnetic  field,  and  the  plasma  and  collision  frequencies.  To  gain  effective  control 
over  attenuation,  mutual  orientation  of  the  antenna  and  imposed  magnetic  field  must  be 
chosen  so  that  the  field  components  are  independent  of  components  of  the  dielectric  tensor 
parallel  to  the  imposed  magnetic  field.  This  approach  has  been  analyzed  285  in  terms  of 
weight  requirements  and  found  to  be  feasible  only  if  a  greater  system  loss  figure,  than  with 
the  other  allevlants  (20  db  rather  than  10  db),  could  be  tolerated.  The  use  of  superconducting 
magnet  has  some  promise  for  short  flights,  but  cooling  equipment  becomes  a  limiting 
factor, 

293 

Lax  et  al  have  studied  breakdown  in  transverse  magnetic  fields,  which  tend  to  worsen 
the  problem  through  two  mechanisms.  The  magnetic  field  enhances  the  transfer  of  energy 
to  the  electrons  and  it  decreases  the  diffusion  losses  perpendicular  to  the  field. 
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The  breakdown  Held  is  reduced  by  very  significant  amounts  in  the  neighborhood  of  the 
cyclotron  resonance  (2.  8  MHz/Gauss),  although  it  is  slightly  increased  at  much  higher 
values  of  magnetic  field.  These  conclusions  are  based  on  studies  of  well  behaved  noble 
gases  and  there  are  only  qualitative  indications  that  they  would  still  hold  for  air,  in 
which  static  magnetic  fields  can  create  serious  nonlinear  behavior.  But  static  fields  are 
probably  not  feasible  alleviation  techniques  for  slender  body  reentry  effects  at  any  rate. 


SECTION  9 


RECOMMENDATIONS 


The  present  state  of  knowledge  of  reentry  effects  on  signal  transmission  is  suffici3ntly  incom¬ 
plete  to  warrant  further  work  on  this  subject.  In  order  to  arrive  at  research  recommenda¬ 
tions  which  are  best  suitea  to  the  systems  problems  to  which  the  tecanolody  is  applied,  it  is 
necessary  to  evaluate  the  results  of  this  study  from  the  standpoint  of  certain  criteria. 
Accordingly,  three  Ck.teria  are  applied  to  assessing  the  importance  of  a  given  effect,  as 
follows : 

a.  The  effect  must  be  severe. 

b.  It  must  be  relatively  difficult  to  predict. 

c.  It  must  he  difficult  to  control  or  alleviate. 

9.1  EVALUATION 

Linear  plasma  effects  are  wholly  dependent  on  the  presence  of  flow-field  ionization.  At  low 
altitudes,  collision  frequency  effects  become  important  in  such  a  way  on  slender  vehicles  as 
to  make  increased  collision  frequency  desirable.  While  these  effects  are  potentially  very 
severe,  it  is  also  possible  to  alleviate  them  by  body  shaping  and  mate  ials  selection.  The 
predictability  of  these  effects  is  relatively  poor,  but  it  is  possible  to  estimate  the  electron 
density  levels  bexow  which  they  are  alleviated. 

Nonlinear  plasma  effects  are  unimportant.  These  effects  deviate  from  linear  effects  as  a 
function  of  the  collision  frequency.  Thus  they  are  present  only  at  low  altitudes,  when  the 
collision  frequency  exceeds  the  signal  radian  frequency,  and  they  are  not  very  severe, 
especially  on  a  slender  cone  vehicle.  The  theoretical  prediction  of  these  effects  is  very 
difficult,  but  they  may  be  alleviated  in  th-i  same  way  as  linear  plasma  effects. 

Antenna  breakdown  effects  are  the  most  Important.  They  score  high  on  all  three  evaluation 
criteria  of  severity,  unpredictability,  and  uncontrollability.  The  six  phenomena  which  have 
been  identified  as  important  influences  on  breakdown  have  varying  degrees  of  predictability 
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and  controllabil'ty,  not  only  in  terms  of  the  breakdown  mechanism  but  also  as  a  function  of 
the  flow-field  prediction  and  control.  For  example,  little  is  presently  known  about  the  high 
gas  temperature  effect  on  the  breakdown  mechanism,  aod  available  techniques  for  describing 
the  thermal  nonequilibrium  effect  have  not  been  applied  to  the  flow  field  properties  of  a  sler.de 
cone  at  high  altitudes. 

Table  9-1  summarizes  the  evaluation  of  reentry  effects  on  antenna  breakdown.  The  items 
marked  "1"  are  the  most  important,  and  the  phenomena  are  listed  in  the  order  of  importance 
based  on  these  criteria.  Temperature  is  most  important  because  so  little  is  known  about  it. 
Gradients  effects  can  be  predicted  by  solving  simple  boundary  value  problems  and  using  the 
resulting  effective  diffusion  length,  but  the  breakdown  prediction  should  be  performed  in  the 
context  of  the  actual  How-field  profiles.  The  convection  effect  requires  more  careful  consi¬ 
deration  of  the  mathematical  criterion  for  breakdown. 


Table  9-1.  Evaluation  of  Flow- Field  Effects  on  Breakdown 


Phenomenon 

Seventy 

Predictability 

Controllability 

Temperature 

1 

1 

1 

Gradients 

2 

2 

1 

Convoction 

2 

2 

1 

Ionization 

1 

3 

2 

Density 

1 

3 

2 

Chemical  Composition 

3 

2 

2 

1  *  most  severe,  least  predictable,  or  least  controllable 

Ionization  effects  on  breakdown  correspond  to  lower  electron  densities  than  for  linear  plasma 
effects,  so  that  they  are  not  so  easily  controlled.  Density  effects  are  well  understood  and 
somewhat  controllable  by  antenna  placement  on  the  vehicle.  Chemical  composition  effects 
are  severe  only  with  heavy  mass  injection  and  then  only  in  the  boundary  layer. 
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In  addition  to  the  consideration  of  the  importance  of  a  given  flow-field  effect,  it  is  necessary 
to  consider  the  Importance  of  the  flow-field  prediction  and  control  problems.  Breakdown 
generally  occurs  at  much  higher  altitudes  than  linear  plasma  effects  on  slender  vehicles,  and 
this  fact  leads  to  difficulties  In  the  flow  field  prediction,  especially  if  the  vibrational  and  elec¬ 
tronic  excited  states  of  the  gas  are  important  la  antenna  breakdown.  In  addition,  the  best  way 
to  account  for  flow-field  gradients  effects  is  to  couple  the  breakdown  theory  to  the  flow-field 
computation. 

9.2  RESEARCH  RECOMMENDATIONS 

Recommendations  for  needed  research  to  improve  the  predictability  and  controllability  of 
reentry  effects  on  breakdown  can  be  derived  from  the  above  evaluation.  Three  topics  are 
recommended  for  further  study;  gas  temperature,  convection,  and  gradients  effects  on  break 
down. 

9. 2. 1  HIGH  GAS  TEMPERATURE 

More  research  is  needed  in  order  to  extend  and  validate  the  data  already  obtained,  for  which 
the  range  of  conditions  of  pressure  and  temperature  is  too  narrow  to  permit  a  good  empirical 
formula  to  be  derived.  The  formula  derived  in  this  study  represent*'  a  very  great  extrapolation, 
which  does  not  correspond  fully  in  the  cold  air  limit  to  the  best  available  theory.  Data  at 
lower  pressures  and  higher  temperatures  are  required  to  fill  this  gap.  The  following  tasks 
should  be  performed  in  order  to  obtain  a  reliable  phenomenological  model: 

a.  Measure  microwave  breakdown  in  air  in  a  shock  tube  at  various  values  of  pressure, 
shock  velocity,  pulse  duration,  and  station  behind  the  shock. 

b.  Measure  the  air  properties  for  the  above  conditions,  including  gas  temperature, 
pressure,  density,  and  electron  density. 

c.  Calculate  the  gas  properties  in  the  above  experiment,  using  the  best  techniques 
available  and  including  coupled  vlbratlon-dlssociatlon-vlbratlon  (CVDV)  effects 
and  allowing  for  arbitrary  preferred  dissociation  from  the  upper  vibrational 
energy  levels. 
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d.  Compare  the  calculated  and  measured  gas  properties  and  adjust  the  proba¬ 
bility  for  preferred  dissociation  to  match  the  experin  ental  data. 

e.  Reduce  the  breakdown  data  to  give  curves  of  normaiitod  ionization  frequency 
versus  normalized  electric  field  as  a  function  of  the  gas  properties,  including 
vibrational  and  electronic  temperatures,  if  necessary'. 

f.  Derive  a  phenomenological  theory  for  the  gas  temperature  effect  on  ionization 
frequency,  using  the  above  data  to  determine  the  values  of  arbitrary  parameters. 

9.2.2  CONVECTION 

A  set  of  experiments  on  microwave  breakdown  in  flowing  gases  at  a  wide  rt  nge  of  conditions 
of  pressure,  velocity,  and  ionization  level  is  needed  to  verify  the  theories  and  show  their 
regions  of  validity. 

A  further  study  of  convection  theory  is  needed  to  resolve  the  discrepancies  among 
the  various  theories,  and  to  interpret  the  results  of  experiments.  The  basic  problem 
is  to  solve  the  coupling  among  diffusion,  convection,  and  space  charge.  The  analysis 
which  was  done  in  this  study  neglected  to  do  this,  using  Instead  simplif  ed  formulas  for 
the  transitional  parameters  without  explicitly  accounting  for  space  charge.  Results  would 
be  compared  with  experiments. 

9. 2. 3  GRADIENTS  AND  BREAKDOWN  CRITERIA 

The  effects  of  flow-field  gradients  on  breakdown  need  to  be  studied  by  solving  the  breakdown 
equations  for  typical  flow-field  profiles  of  gas  velocity,  density,  temperature,  and  ionization, 
using  realistic  antenna-field  distributions.  This  work  should  be  done  in  conjunction  with 
a  study  of  breakdown  criteria  for  slender  body  flow-field  effects.  In  other  words,  the 
antenna  Impedance  and  radiated  power  3hould  be  calculated  as  a  function  of  input  power  in 
conjunction  with  the  breakdown  equations.  To  do  this  properly  may  mean  coupling  the 
flow-field,  the  breakdown,  and  the  propagation  equations  together  into  a  single  computer 
program. 
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9.2.4  FLIGHT  EXPERIMENTS 

Flight  experiments  are  needed  to  provide  realistic  results  on  reentry  effects  on  antenna 
breakdown.  Two  types  of  experiments  can  be  fruitful.  Systems  tests  would  be  useful  In 
testing  the  breakdown  behavior  of  a  given  particular  antenna  system,  and  basic  data 
experiments  would  be  designed  to  show  the  Importance  of  one  or  another  of  the  six  basic 
reentry  phenomena.  Recommendations  for  flight  experiments  are  given  In  detail  In 
Volume  2  of  this  report. 
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SECTION  10 


CONCLUSION 


Ti  °  most  important  effect  of  slender  cone  vehicle  reentry  on  signal  transmission  Is 
the  perturbation  of  the  breakdown  characteristics  of  the  antenna.  The  breakdown 
phenomenon  is  important  in  terms  of  the  severity,  unpredictability,  and  uncontrol¬ 
lability  of  reentry  effects.  Not  all  of  the  six  presently  identifiable  phenomena  of 
the  reentry  environment  are  equally  important,  but  the  fact  that  there  are  so  many 
factors  contributing  to  the  effects  adds  to  their  unpredictability  and  the  lack  of  ability 
to  control  them.  Linear  plasma  effects  are  controllable  and  relatively  predictable 
on  a  slender  cone  vehicle,  so  that  their  importance  is  less  than  reentry  antenna 
breakdown  effects.  Nonlinear  effects  are  unimportant,  because  they  are  not  very 
severe,  although  they  are  the  least  predictable. 

Signal  transmission  systems  for  reentry  vehicles  can  be  designed  to  anticipate  the 
reentry  effects,  to  a  certain  extent.  Slender  cone  vehicles  which  reenter  the 
atmosphere  at  ICBM  velocity  can  be  designed  to  eliminate  linear  (and  nonlinear) 
plasma  effects  under  the  proper  design  constraints.  The  designer  must  have  the 
freedom  to  select  the  vehicle  shape,  materials,  and  antenna  location.  A  good 
broadband  antenna  design  should  be  used  also.  Of  course  the  plasma  effects  cannot 
be  alleviated  by  these  means  below  the  altitude  of  reentry  to  which  the  vehicle  shape 
and  materials  survive,  unless  the  vehicle  velocity  is  sufficiently  low  when  the  shape 
or  materials  advantage  is  lost. 

No  such  situation  exists  in  the  case  of  reentry  effects  on  antenna  .breakdown.  Here 
the  alleviation  of  flow  field  plasma,  as  shown  above,  has  some  beneficial  effect,  but 
not  enough  is  known  about  various  other  environmental  influences  on  breakdown 
to  permit  their  reliable  alleviation,  especially  in  the  high  altitude  portion  of  reentry. 
More  research  in  flow-field  effects  and  antenna  breakdown  phenomenology  is  needed. 
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